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Finally, this manual Is accompanied by a disk containing the following 


INSTALl.BAT 

PCCARES.INI 

PCCARES. FOR 

PCCARES.EXE 

PCCARESB.EXE 

PCCARES.DAT 

TEMPLET. INP 


Pc'cARF? 1nStd11 batch file 

PC-CARES Initialization file 

PC-CARES FORTRAN source code 

PC-CARES DOS execution file 

Pr'fARF^ D0S> os/2 bound execution file 

pr padcc ? ample problem input file 
PC-CARES input templet file 


UP by ,”fl ng addl tlona?^cop?es°us 1 ng^^ OS/ 2 ’coPY^ ^ S !) 0Uld be backed 

verston^f H&SS F " £^s 

the Lahey FORTRAN compeer), fffs E ? ec S2 el E*! MiCr ° Soft W<* BASIC o? 
get a copy of the newer version for which thu ^° U contact Y°ur source and 
BASIC version does not calcu ate the "T al was bitten. The older 

Lahey FORTRAN verslo? . 1 * on?^ run on coefflcie "‘ and the 

Further, neither will execute under CK /9 S L wi * h math ^-processors, 
advantage to upgrade. b tect mode, so It will be to your 


You may run PC-CARES off a floppy disk or inctaii 
It Is recommended that you put the PC-CARES rndl ?! I U on your hard d1sk - 
when you Install It. For your convenient V iZl i"Jj S 0 r n separate directory 
Included on the distribution disk which is uspH L f n’ ^STALL.BAT, has been 
mand prompt with the default drive let to \L U fo V^ s: At the D0S com- 
s 1 mp 1 y type INSTALL C: where C Is Ihe Srfvt / V °l the distr ’bution disk 

The batch file will create thed 1 rectorv d r!\prrJJj? nat 2 r of your hard d1sk - 
the distribution disk In this directory 7 ??? p1ace the conten ts of 

users ^INSTALL^^im Wo OS™* FwOS/P <S 

rename INSTALL.BAT to INSTALLED and exec^e^ SZX'&H c^nV^t. 


2.0 PROGRAM CAPABILITY AND DESCRIPTION 


Version S.O^ompner^fJj^hJ^^ Microsoft FORTRAN 

tlon supported by this compiler for the JS M /isS^° nS Si dynam1c arra ^ al loca- 
te current version of this program Jni makfut o? r a 0S/ K° Perat,ng systems - 
is installed on your compute?. V 1 t u ^ ^-processor if one 


user to choose both the^umber 1 of C constan? U temDe °t the f compner allows the 
maximum number of test specimens D e? tem?^ P ^ re fracture sets and the 
these parameters In the k! 5 wo?d d e 11 * at runt1me by setting 

This Initialization file also allows vou to °?»w f1 ,e Called PC CARES. INI . 

temperature'fra^^ ^ ? s No?e "^It 

-> .. . MM s; r-JM yxts.s.rs*s jssar* - 


ho»e,e? de 5 ufL t ^h”“l™M )r ??:?"!_ s ! 2e '5 640 kB. 


This limit mav hp iriwor 


PRINT spooler, or other such utilities. Remember if you want to check out the 
amount of available memory under DOS, the CHKDSK command gives this information 
along with the report on the disk medium. For more information on CHKDSK con- 
sult your DOS reference manual. You may have to remove all the memory resident 
programs from memory. If you execute memory resident programs from your 
AUTOEXEC.BAT file, you may want to delete or comment out those commands before 
you reboot to achieve as much free memory as you can. If you are using a LIMS 
PC at NASA Lewis, it is recommended that you consult the LIMS manual or Compu- 
ter Services before altering AUTOEXEC.BAT. If you should attempt to allocate 
more memory than is available, the program will simply halt and display an 
error message, so there is no harm in experimenting with the memory allocation. 

Under OS/2 versions 1.1 and 1.2 you are limited only by the virtual memory 
capabilities of this operating system and the architecture of the 80286 chic 
which permits processes (programs) up to 16 mB of virtual address space. (Note 
that if you have a 80386 system with OS/2 versions 1.1 or 1.2, you are still 
limited to 16 mB as these versions operate the 386 chip in 286 mode.) In prac- 
tice, you are limited only by the sum of the amount of physical RAM and the 
amount of hard disk space available on your system. If you are running other 
processes concurrently, this amount will be reduced accordingly. 

The primary function of PC-CARES is statistical analysis of the data 
obtained from the fracture of simple, uniaxial tensile or flexural specimens 
and estimation of the Weibull and Batdorf material parameters from this data. 


The weakest-link mechanism is expressed with the classical Weibull two- 
parameter formulation, which, for volume flaw reliability, is 



and for surface flaw reliability is 



( 2 . 1 ) 


( 2 . 2 ) 


where P s is the survival probability and a is the applied uniaxial tensile 
stress. Here V is the volume of stressed material, and A is the area. The 
subscripts V and S denote parameters that are a function of material volume 
and surface area, respectively. The scale parameter o 0 has dimensions of 
stress x (volume)^ m V or stress x (area)^ m S. The scale parameter corresponds 
to the stress level at which 63.2 percent of specimens with unit volume or area 
would fracture. The shape parameter (or Weibull modulus), denoted by m, is a 
dimensionless quantity and measures the degree of strength dispersion of the 
flaw di stri bution . 


In general, the parameters are obtained from the fracture stresses of many 
specimens (30 or more are recommended) whose geometry and loading configura- 
tions are held constant. Solutions for the three-point modul us-of-rupture 
(MOR) bending bar, four-point MOR bending bar (ref. 4), and the pure tensile 
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specimen (ref. 5) maintained at a constant specified temperature have been 
incorporated into the PC-CARES program. Since the material parameters are a 
function of temperature, different constant-temperature data sets can be input 
and the corresponding parameter estimates will be calculated. The amount of 
specimens per each constant temperature data set and the total number of these 
sets is only limited by the computer memory available to the program. In 
addition, each specimen can be identified by its mode of failure, either 
volume flaw, surface flaw, or some other mode so that parameter estimates for 
competing failure modes can be obtained. The statistical accuracy of the 
parameter estimates compared with the true material parameters depends on the 
number of specimens tested, assuming that the true distribution is a Weibull 
distribution . 

Figure 2.1 shows the flowchart for the calculation of the statistical 
strength parameters of the two-parameter Weibull distribution for volume-flaw 
and surface-flaw-induced fracture, with complete (single mode) or censored 
(multiple mode) samples, and the calculation of other statistical quantities. 
Following the input of specimen geometry, fracture stresses, and respective 
flaw origins, PC-CARES will first identify any potential bad data (outliers). 
The outlier test developed by Stefansky (ref. 6) and subsequently used by Neal, 
Vangel , and Todt (ref. 7) is incorporated into the program. Although the tech- 
nique is based on the normal distribution and, therefore, its application to 
the Weibull distribution is not rigorous, it serves as a guideline to the user. 
Data detected as outliers are flagged with a warning message, and any further 
action is left to the discretion of the user. 

Weibull parameter estimates are obtained for the specimen surface and/or 
volume as requested by the user, taking into account the fracture origin data 
also supplied by the user. Biased estimates of the Weibull shape parameter 
and characteristic strength are obtained from either least-squares analysis or 
the maximum likelihood method for complete samples and/or censored samples. 
PC-CARES uses the Weibull log-likelihood equations given in Nelson (ref. 8) 
and the rank increment adjustment method described by Johnson (ref. 9), for 
complete and censored statistics. 

Because the estimates of parameters are obtained from a finite amount of 
data, they contain an inherent uncertainty that can be character i zed by bounds 
in which the true parameters are likely to lie. Methods have been developed 
to evaluate confidence limits that quantify this range with a level of proba- 
bility as a function of sample size. For the maximum likelihood method witn 
a complete sample, unbiasing factors for the shape parameter m, and 5 and 
95 percentile confidence limits for m and the characteristic strength oq , 
are provided (ref. 10). The characteristic strength, or characteristic 
modulus of rupture, is similar to the Weibull scale parameter except that it 
includes the effect of the total specimen volume or area. For a censored 
sample, an asymptotic approximation of the 90-percent confidence limits is 
calculated. No unbiasing of parameters or estimation of confidence limits is 
given when the least-squares option is requested. 

The ability of the parameter estimates to reasonably fit the empirical 
data is measured with the Kolmogorov-Smi rnov (K-S) and Anderson-Darling (A-D) 
goodness-of-f It tests. These tests are extensively discussed by D'Agostino and 
Stephens (ref. 11). The tests quantify discrepancies between the experimental 
data and the estimated Weibull distribution by a significance level associated 
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with the hypothesis that the data were generated from the proposed distribu- 
tion. The A-D test is more sensitive than the K-S test to discrepancies at low 
and high probabilities of failure. The Kanofsky-Srini vason 90-percent confi- 
dence band values (ref. 12) about the Weibull line are given as an additional 
test of the goodness-of-f i t of the data to the Weibull distribution. 

After the shape and characteristic strength parameters are estimated and 
analyzed, PC-CARES calculates the other material parameters. The biased esti- 
mate of the shape parameter m and the estimated characteristic strength a@ 
are used along with the specimen geometry to calculate the Weibul 1 _s c a 1 e 
parameter <j 0 . The Batdorf normalized crack density coefficient kg, which is 
explained in the appendix section THEORY, is computed from the selected frac- 
ture criterion, crack geometry, and the biased estimate of the shape parameter. 
Figure 2.2 shows the fracture criteria and flaw geometries available to the 
user which must be specified in order to calculate the Batdorf crack density 
coefficient. If the user selects to calculate the Batdorf crack density coef- 
ficient by setting it to the value that is the solution for the normal stress 
fracture criterion, then the user need not specify the fracture criterion and 
the crack geometry. 

The simple P I A fracture theory does not require a specific crack geometry, 
and for uniaxial stress states, it reduces to equation (2.1) or (2.2), which- 
ever is appropriate. The Weibull normal stress averaging method is also inde- 
pendent of crack geometry, since it only considers impending mode I (opening 
mode) crack growth, and neglects mode II (sliding mode) and mode III (tearing 
mode) effects. Batdorf's fracture theory can be used with several different 
mixed-mode fracture criteria and crack geometries. The combination of a par- 
ticular flaw shape and fracture criterion results in an effective stress equa- 
tion involving far-field principal stresses in terms of normal and shear 
stresses acting on the crack plane. The coplanar crack extension criterion for 
shear-sensitive materials available in PC-CARES is the total strain energy 
release rate theory. Out-of-plane crack extension criteria are approximated by 
a simple semi-empirical equation (refs 13 and 14). This equation involves a 
parameter that can be varied to model the maximum tangential stress theory, the 
minimum strain energy density criterion, the maximum strain energy release rate 
theory, or experimental results . For comparison, Griffith's maximum tensile 
stress analysis for volume flaws is also included. The highlighted boxes in 
figure 2.2 show the recommended fracture criteria and flaw shapes. 
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MAIN PROGRAM 



FIGURE 7.1. - FLOWCHART TOR CALCULATION OF MATERIAL STATISTICAL STRENGTH PARAMETERS. 


6 













X'! 

/ t 
/ / 

<" / 




FRACTURE CRITERIA 


FIGURE 2.2. - AVAILABLE FAILURE CRITERIA AND CRACK SHAPES. 
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3.0 PC-CARES INPUT INFORMATION 

PC-CARES uses two files during Its execution, the Initialization file 
(PCCARES.INI) and the Input file (an example Is Included on the PC-CARES disk 
called PCCARES.DAT) . Using the initialization file with PC-CARES is optional. 
For preparation of both the. initialization file and the PC-CARES Input file, it 
Is assumed that the user has access to a DOS or OS/2 full screen text editor 
with block editing capabilities. 

Input to PC-CARES for both the initialization file and the input file Is 
keyword driven. The keywords can be present In any order within each input 
section, but they must start In the first column of the file. Examples of both 
the Initialization file and the input file are given in figures 3.1 to 3.7. 

Note that underneath each keyword a location is given that specifies where the 
data value or values are input. An explanation of each keyword is provided to 
the right, and a list of available choices Is given, if applicable. If inte- 
ger Input is required, then the input field is between two asterisks (*), and 
entries must be right justified. Real number input is read In an F10.4 format, 
and asterisks are not present to define the field width. A maximum of 20 lines 
between keywords Is allowed before an error message is generated, and there- 
fore, the user can Insert short notes as desired. 


3.1 PC-CARES Initialization File Description 

Upon execution, PC-CARES first attempts to locate the file called 
PCCARES.INI In the directory from which PC-CARES was executed. If the file Is 
not present a warning message Is issued to the screen and the default para- 
meters for filenames and array allocation sizes are used. Note that if you 
intend to make use of the initialization file, it must be located in the same 
directory as the PC-CARES execution file, otherwise the defaults are used. 

The Initialization file, PCCARES.INI, consists of two categories of input: 
(1) File Control Input and (2) Dynamic Memory Control Input. The File Control 
Input consists of two path/filenames (up to 20 characters) which PC-CARES uses 
as the Input and output filenames, so that the user can specify different 
directory and/or filenames for the input and output files to/from the PC-CARES 
program. The Dynamic Memory Control Input can contain two indices: one to 

control the number of constant temperature data sets and the other to control 
the maximum number of fracture stress data per constant temperature data set. 
PC-CARES continues to read the Initialization file until either It reads the 
$ENDI keyword or end of file is encountered. Note that keywords that are not 
found assume default values. 

The user is advised to keep an unaltered copy of the PCCARES.INI file as 
a backup. Details on specific Input preparation are described In 3.2.1 File 
Control Input and 3.2.2 Dynamic Memory Control Input sections of this manual. 
Each keyword Is discussed briefly In these sections, and the format field for 
the Input Is denoted In parentheses next to the keyword as well as the default 
val ues . 
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3.2 PC-CARES Initialization File Preparation 
3.2.1 File Control Input 

The File Control Input optionally controls which files PC-CARES uses for 
input and output, overriding the default values. Each file/path name may be up 
to 20 characters in length. For example see figure 3.1. In this example, if 
PC-CARES was located in the C:\PCCARES directory the program accesses the input 
file PPCARES.DAT in the C:\PCCARES\DATA directory and produces the output file 
PCCARES.OUT in the C:\PCCARES\OUTPUT directory. Please note that the output 
directory must be created prior to program execution. The keywords pertaining 
to the File Control Input section of the PC-CARES Initialization file are 
defined as follows (the input format field is in parentheses): 

INFILE (A20) ; Default : INFILE = 'CARES.DAT' 

INFILE controls where PC-CARES reads its input. 

OUTFILE (A20); Default : OUTFILE = 'CARES. OUT' 

OUTFILE controls where PC-CARES sends its output. 


3.2.2 Dynamic Memory Control Input 

The Dynamic Memory Control Input optionally controls the maximum number of 
constant temperature test sets and the maximum number of fracture stresses per 
temperature set, and hence the amount of array space the program allocates. 

This dynamic array allocation makes the compiled program smaller and gives the 
user the flexibility to determine how much memory the program will use during 
execution. Setting IMAXT = 20 and IMAXF = 200 is about the maximum allocation 
for the PC-CARES running under DOS with its maximum 640K program size. Fig- 
ure 3.1 shows an example of the Dynamic Memory Control Input section of the 
PC-CARES Initialization file. The keywords IMAXF and IMAXT are defined as fol- 
lows (the Input format field is in parentheses): 

IMAXF (2X.I5) ; Default : IMAXF = 150 

IMAXF controls the maximum number of fracture stress data per constant 
temperature set. 

IMAXT (2X.I5); Default : IMAXT = 10 

IMAXT controls the maximum number of constant temperature fracture stress 

sets. 


3.3 PC-CARES Input File Description 

For the PC-CARES input file, two categories of Input are required for 
execution: (1) Master Control Input and (2) Material Control Input (which 

Includes temperature-dependent material data). The Master Control Input Is a 
set of control indices which directs the overall program execution. It speci- 
fies the number of brittle material statistical characterizations, the number 
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of Guassian quadrature points to be used to perform numerical integration, and 
whether the fracture stresses should be reprinted in the output. The Material 
Control Input consists of control indices and either the data required to esti- 
mate the statistical material parameters or direct input of the Weibull statis- 
tical parameter values themselves, for various temperatures to calculate the 
Batdorf crack density coefficient only. This input category includes the 
choices of fracture criteria and flaw shapes shown in figure 2.2. The choice 
of fracture criteria and flaw shape is required to calculate the Batdorf crack 
density coefficient since it is a function of these conditions as well as the 
fracture stress data. 


3.4 PC-CARES Input File Preparation 

To control the execution of the PC-CARES program, the user must prepare an 
input file consisting of the Master Control Input and the Material Control 
Input. On the disk provided with the program is a file called TEMPLET. INP that 
can be used to construct an input file for a particular problem. The Master 
Control Input always comes at the beginning of a file. 

After reading the initialization file (if one is present) PC-CARES contin- 
ues execution by searching for the keywords associated with the Master Control 
Input in the PC-CARES input file. The end of the Master Control Input occurs 
when the $ENDX keyword is encountered. Following the Master Control Input, 
PC-CARES searches for keywords specific to the Material Control Input. The 
$ENDM and $ENDT keywords signal the end of two different sections of the Mate- 
rial Control Input. Keywords not found between SEND intervals may assume 
default values. Because PC-CARES has a multiple material capability, each 
section of input for a particular material is separated by a SENDT card. The 
TEMPLET. INP file has only two materials characterized. Modifying the file for 
more materials involves block copying sections of the original file, appending 
them to the end of the file, and modifying the copied input values accordingly. 
The user is advised to keep an unaltered copy of the TEMPLET. INP file as a 
backup. Details on specific input preparation are described in the Master Con- 
trol Input and Material Control Input sections of this manual. Each keywo-d is 
discussed briefly in these sections, and the format field for the input is 
denoted in parentheses next to the keyword as well as the default value where 
appl i cable . 


3.4.1 Master Control Input 

The Master Control Input section from the TEMPLET. INP file is reproduced 
in figure 3.2. If parameter keywords in the input file are omitted, they 
assume their default values. 

NGP (4X,I2); Default : NGP = 15 

NGP controls the number of Gaussian integration points that are used in 
the reliability calculations. An entry of 30 will give better accuracy but at 
the penalty of larger CPU requirements than an entry of 15. There are also 
other options for NGP in PC-CARES, but users should not specify less than 
15 Gaussian points. 
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NMATS, NMATV (4X,I2); Defaults : NMATS = NMATV = 0 

The keyword NMATS represents the number of materials for which surface 
flaw analysis is performed. NMATV represents the number of materials for wnich 
volume flaw analysis is performed. A component consisting of one material may 
have one set of statistical material parameters to characterize the surface and 
another set for the volume, for which NMATS = 1 and NMATV = 1. Statistical 
material parameters are a function of processing, microstructure, and environ- 
ment. The PC-CARES program is capable of analyzing a single material with 
multiple statistical material characterizations or many materials with multiple 
statistical material characterizations. For example, if a single material has 
two different surface finishes, then NMATS = 2 is used because two different 
sets of statistical material parameters are required. 

TITLE (72A1 ) 

The input associated with the TITLE keyword is reproduced in the program 
output for problem identification. 

$ENDX 


$ENDX signifies the end of the MASTER CONTROL INPUT. 


3.4.2 Material Control Input 

A sample of the Material Control Input section from the TEMPLET. INP file 
is reproduced in figures 3.3 and 3.4. The figures are an example of the input 
required for PC-CARES to estimate the volume flaw statistical material parame- 
ters from experimental fracture data. In figures 3.5 and 3.6 an example of 
the input needed to estimate the surface flaw statistical material parameters 
from experimental fracture data is shown. Note that the Material Control 
Input actually consists of two different data partitions. Figures 3.3 and 
3.4, for example, make up a single section of the Material Control Input. In 
figure 3.3, the control indices, material constants, and geometric variables 
necessary to calculate volume flaw statistical parameters are shown. In fig- 
ure 3.4, the temperature-dependent fracture data are given. The temperature- 
dependent fracture data (MOR), or temperature-dependent values of the Wei bull 
shape and scale parameters (PARAM), are always placed immediately following the 
control indices for that material. 

The total number of Material Control Input sections is equal to the sum 
of NMATS + NMATV from the Master Control Input. Note that keywords that are 
not found assume default values. 


3.4.2. 1 Material and Specimen Dependent Data 

The following keywords are the control indices, material indices, and geo- 
metric variables necessary for calculation of volume and surface flaw statisti- 
cal parameters as shown in figures 3.3 and 3.5. 
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C (F10.4); Default : NONE 


If ID2S = 5 or ID2V = 5 (i.e., if She tty ' s mixed-mode fracture criterion 
is selected), then the value of the empirical constant C, denoted by the key- 
word C, must be specified. If this criterion is not selected, this input is 
ignored and can be deleted. 

DL1 , DL2 , DH, DW (F10.4); Default : NONE 

If ID1 = 2 or 5 (i.e., if statistical material parameters are to be deter- 
mined from four-point MOR fracture specimens), then the specimen dimensions 
must be input. DL1 represents the length between the two outer symmetrical 
loads. DL2 is the length between the two inner central loads. DH is the total 
height of the test specimen cross section. DW is the total width of the test 
specimen cross section. By setting DL2 equal to zero, data obtained from 
three-point bend tests can also be used to obtain appropriate Weibull 
parameters . 

ID1 ( 4X ,11); Default : NONE 

ID1 is a control index for specifying the form of the data to be input 
for obtaining the statistical material parameters. Either the Weibull shape 
and scale parameters are directly specified or experimental fracture data are 
input. The fracture data can be either from four-point modul us-of-rupture 
bend bars or from tensile test specimens. If the fracture data are assumed to 
be all from one failure mode (all volume flaws or all surface flaws), then 
IDl = 1 or 2 can be chosen. If ID1 = 1 or 2, then fracture origins are not to 
be input with the specimen fracture stresses, and PC-CARES assumes that the 
fracture origins are consistent with the ID4 input index. If IDl = 4 or 5, 
then fracture origins must be supplied with the fracture data. 

ID2S.ID2V (4X,I1); Default : NONE 

These control indices are for selection of a fracture criterion. ID2S is 
for a surface flaw fracture criterion (see fig. 3.5). ID2V is for a volume 
flaw fracture criterion (see fig. 3.3). If ID4 = 1, then ID2V should be speci- 
fied. If ID4 = 2, then ID2S should be specified. If both ID2S and ID2V are 

specified in the same input section, the entry not consistent with the ID4 
index is ignored. Shetty's mixed-mode fracture criterion is recommended for 
both surface and volume flaw analysis. 

ID3S , ID3V (4X.I1 > ; Default : NONE 

The ID3S and ID3V control indices are for selection of a crack geometry. 
ID3S is for surface flaw geometry (see fig. 3.5). ID3V is for volume flaw 
geometry (see fig. 3.3). If ID4 = 1, then ID3V should be specified. If 

ID4 = 2, then ID3S should be specified. If both ID3S and ID3V are specified in 

the same input section, the entry not consistent with the ID4 index is ignored. 
The penny-shaped crack is recommended for volume flaw analysis and the semicir- 
cular crack is recommended for surface flaw analysis. 
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ID4 (4X.I1); Default : NONE 

ID4 controls the calculation of volume- or surface-based statistical mate- 
rial parameters. From the fracture data supplied, the Weibull shape and scale 
parameters along with the normalized Batdorf crack density coefficient are 
estimated. If the Weibull shape and scale parameters are directly input, then 
the normalized Batdorf crack density coefficient is calculated. 

IKBAT (4X.I1 ) ; Default : IKBAT = 0 

IKBAT selects the method of calculating the normalized Batdorf crack 
density coefficient. If IKBAT = 0, then the crack density coefficient is set 
to the value that is the solution for the normal stress fracture criterion, 
regardless of the fracture criterion and crack geometry selected by the user, 
and therefore the ID2S and ID3S (or ID2V and ID3V) indices do not have to be* 
specified for surface or volume flaw analysis, respectively. If IKBAT = 1, 
then the crack density coefficient is calculated based on the fracture cri- 
terion and crack geometry selected by ID2S and ID3S or by ID2V and ID3V. 

IKBAT = 0 gives more conservative reliability predictions and usually agrees 
more closely with test data than IKBAT = 1 does; it is therefore recommended 
as the best choice unless specific data exist that indicate otherwise. 

MATID (IX, 17) 

MATID is the material identification number that is associated with the 
statistical material parameter data. 

MLORLE (4X, 1 1 ) ; Default : MLORLE = 0 

MLORLE is the control index for the method of estimation of the Weibull 
shape parameter m and characteristic strength <j@ from experimental fracture 
data. MLORLE is ignored if the Weibull shape and scale parameters are directly 
i nput . 

PR (F10.4) ; Default : PR = 0.25 

PR is Poisson's ratio. It is assumed to be temperature independent. 

TITLE (72A1 ) 

The input associated with the TITLE keyword is reproduced in the program 
output for material identification. 

VAGAGE (F10.4); Default : NONE 

VAGAGE is the gage volume or area of a tensile test specimen. If ID4 = 2 
and ID1 = 1 or 4 (i.e., if surface flaw analysis is specified and the statis- 
tical material parameters are to be determined from simple tension tests), 
then the gage surface area of the specimen must be specified. If ID4 = 1 and 
ID1 = 1 or 4 (i.e., if volume flaw analysis is specified and the statistical 
material parameters are to be determined from simple tension tests), then the 
gage volume of the specimen must be specified. 
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$ENDM 


$ENDM signifies the end of a section of the Material Control Input. The 
temperature-dependent specimen fracture data or the Weibull shape and scale 
parameters are assumed to immediately follow. 


3. 4. 2. 2 Temperature-Dependent Fracture or Statistical Material Parameters Data 

Immediately following the $ENDM keyword, which signals the end of the 
material and specimen dependent data, the temperature-dependent experimental 
fracture data or Weibull shape and scale parameters are input. Data for up to 
IMAXT different temperatures can be specified, where IMAXT is the maximum num- 
ber of constant temperature data sets as described in section 3.2.2 Dynamic 
Memory Control Input, of this manual. Figures 3.4 and 3.6 show examples of the 
input for experimental fracture stress data. Figure 3.7 shows an example of 
the input for the Weibull shape and scale parameters. 

MOR (3A1 .3E18.10) or (3E18.10) 

MOR indicates that experimental fracture stresses will be input. Frac- 
ture stresses can be input in any order, with a maximum of IMAXF specimen fail- 
ure stresses input for each temperature. There are two styles of input. If 
ID1 = 1 or 2 ( 1 . e . , if the fracture data are assumed to be a complete sample), 
then fracture stresses only are input and the Input format is 3E18.10 as shown 
In figure 3.6. Referring to figure 3.4, if ID1 = 4 or 5 ( 1 . e . , if the frac- 
ture origins and the fracture stresses are to be input), then the Input format 
Is 3A1 , 3E18.10. The 3A1 represents three fields of single alphanumeric char- 
acters. This field is for fracture origin Input. An "S" Indicates a surface 
flaw origin. A "V" represents a volume flaw origin. A "U" indicates an 
unknown flaw origin. Each fracture stress has a corresponding fracture origin. 
In figure 3.4, each line of fracture data consists of three fracture origins 
followed by their respective failure stresses. Fracture data values should be 
unique, and multiple identical values should not be input (change the values 
slightly). 

NUT (3X, 13) 

NUT is the sample size of the experimental fracture data for the tempera- 
ture Indicated by TDEG. NUT is specified If ID1 does not equal 3 (statistical 
material parameters are not being directly Input). Different numbers of speci- 
mens are permitted at different temperatures. 

PARAM (2E18.10) 

PARAM signals that the Weibull shape and scale parameter will be input for 
the temperature Indicated by TDEG. Referring to figure 3.7, the Weibull shape 
and then the Weibull scale parameters are entered at the indicated space with a 
format of 2E18.10. The Weibull shape parameter Is dimensionless. The Weibull 

scale parameter has units of stress x (volume) ' m V for volume flaw analysis 
and units of stress x (area) 1/ms for surface flaw analysis. Note that ID1 = 3 
must be specified. 
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TDEG (F10.4) 


TOEG is the input keyword for the temperature of the fracture data or of 
the statistical material parameters that immediately follow. Temperature can 
be specified in any units desired and is used for identification purposes only 
in the PC-CARES code. 

$ENDT 

$ENDT signals the end of the temperature-dependent data. 

Another section of the MATERIAL CONTROL INPUT follows, if required. 
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— PC-CARBS INITIALIZATION FILE — 


FILE CONTROL INPUT 


INFILE : PATH AND FILE NAME FOR CARES MASTER CONTROL INPUT FILE (LUA) 

(LENGTH OF PAT11/FILE NAME MUST BE < 20 CHARACTERS) 

DATA\PCCARES . DAT 

(DEFAULT: ’PCCARES.DAT’) 

OUTFILE : PATH AND FILE NAME FOR CARES GENERAL OUTPUT FILE (LUB) 

— (LENGTH OF PATH/FILE NAME MUST BE < 20 CHARACTERS) 

DATA\PCCARES . OUT 

(DEFAULT: ’PCCARES.OUT’) 


DYNAMIC MEMORY CONTROL INPUT 


IMAXF : MAXIMUM NUMBER OF FRACTURE STRESSES THAT CAN BE INPUT PER RUN 

(IMAXF > 0. NOTE - THIS PARAMETER CONTROLS THE AMOUNT OF 

*00150* ARRAY SPACE DYNAMICALLY ALLOCATED AT RUN-TIME.) 

(DEFAULT: 150) 

IMAXT : MAXIMUM NUMBER OF TEMPERATURE SETS THAT CAN BE INPUT PER RUN 

(IMAXT > 0. NOTE - THIS PARAMETER CONTROLS TnE AMOUNT OF 

*00010* ARRAY SPACE DYNAMICALLY ALLOCATED AT RUN-TIME.) 

(DEFAULT: 10) 


$ENDI : END OF INITIALIZATION FILE 


FIGURF 3.1. - IT CARTS INI1 IAI l/AI ION I II I . 
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— PC-CARES TEMPLET INPUT FILE — 


************* 






MASTER CONTROL INPUT 


TITLE : PROBLEM TITLE (ECHOED IN PC-CARES OUTPUT) 

SAMPLE INPUT FOR PC-CARES USERS MANUAL 


: NO. OF MATERIALS FOR SURFACE FLAW ANALYSIS 

(NMATS+NMATV < 101) 

*01* (DEFAULT: NMATS = 0) 


lTV : NO. OF MATERIALS FOR VOLUME FLAW ANALYSIS 

(NMATS+NMATV < 101) 

*01* (DEFAULT: NMATV = 0) 


IPRINT : CONTROL INDEX FOR STRESS OUTPUT 

(DEFAULT: IPRINT = 0) 

*1* 0 : DO NOT PRINT FRACTURE DATA 

! . PRINT FRACTURE DATA 

NOP : NO. OF GAUSSIAN QUADRATURE POINTS (15 OR 30) 

(DEFAULT: NGP =15) 


************* 

$ENDX 


END OF MASTER CONTROL INPUT 


******** 




FIGURE 3.2. - PC-CARES MASTER 


CONTROL INPUT. 
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MATERIAL CONTROL INPUT 

TITLE : MATERIAL TITLE (ECHOED IN PC-CARES OUTPUT) 

SAMPLE INPUT FOR PC-CARES USERS MANUAL 


MATID : MATERIAL IDENTIFICATION NUMBER 

(NO DEFAULT) 

*0000300* 


IDI 


* 5 * 


CONTROL INDEX FOR EXPERIMENTAL DATA 
(NO DEFAULT) 

1 : UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

2 : FOUR-POINT BEND TEST DATA 

3 : DIRECT INPUT OF THE WEIBULL PARAMETERS, M AND SP 

(SHAPE PARAMETER AND SCALE PARAMETER) 

4 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

5 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

FOUR-POINT BEND TEST DATA 


: CONTROL INDEX FOR VOLUME OR SURFACE FLAW ANALYSIS 

(NO DEFAULT) 

*1* 1 : VOLUME 

2 : SURFACE 


ID2V 


*5* 


ID3V 


* 2 * 


IKBAT 


* 0 * 


CONTROL INDEX FOR VOLUME FRACTURE CRITERION 
(NO DEFAULT) 

1 : NORMAL STRESS FRACTURE CRITERION 

(SHEAR-INSENSITIVE CRACK) 

2 : MAXIMUM TENSILE STRESS CRITERION 

3 : COPLANAR STRAIN ENERGY RELEASE RATE CRITERION 

(G SUB T) 

4 : WEIBULL PIA MODEL 

5 : SHETTY’S SEMI-EMPIRICAL CRITERION 

CONTROL INDEX FOR SHAPE OF VOLUME CRACKS 
(NO DEFAULT) 

1 : GRIFFITH CRACK 

2 : PENNY-SHAPED CRACK 

CONTROL INDEX FOR METHOD OF CALCULATING BATDORF CRACK 
DENSITY COEFFICIENT (K SUB B) FROM TEST DATA 
(DEFAULT: IKBAT = 0) 

0 • SHEAR-INSENSITIVE METHOD (MODE I FRACTURE ASSUMED) 

1 • SHEAR-SENSITIVE METHOD (FRACTURE ASSUMED TO OCCUR 

ACCORDING TO THE FRACTURE CRITERION AND CRACK SHAPE 
SELECTED BY THE ID2 AND ID3 INDICES) 


PR : POISSON’S RATIO 

(DEFAULT: PR = 0.25) 

00000.2500 


C • CONSTANT FOR SHETTY’S SEMI -EMPIRICAL MIXED-MODE FRACTURE 

CRITERION (KI/KIC) + (KII/ (C*KIC) ) **2 = 1 

00000 8000 OBSERVED VALUES RANGE FROM 0.8 TO 2. (REF. D.K. SHETTY) 
1 NOTE: AS C APPROACHES INFINITY, PREDICTED FAILURE 

FIGURE 3.3. - PC-CARES MATERIAL CONTROL INPUT FOR VOLUME FLAW ANALYSIS. 
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PROBABILITIES APPROACH NORMAL STRESS CRITERION VALUES 
(DEFAULT C = 1.0) 

MLORLE : CONTROL INDEX FOR METHOD OF CALCULATING WEIBULL 

PARAMETERS FROM THE EXPERIMENTAL FRACTURE DATA 

*1+ (DEFAULT: MLORLE = 0) 

0 . MAXIMUM LIKELIHOOD 

1 : LEAST- SQUARES LINEAR REGRESSION 

DH : HEIGHT OF THE FOUR POINT BEND BAR 

(NO DEFAULT) 

00000.0710 


DL1 : OUTER LOAD SPAN OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

00001.0000 


DL2 : INNER LOAD SPAN OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

00000 . 5000 


DW : WIDTH OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

00000.1477 


SENDM : END OF TEMPERATURE INDEPENDENT MATERIAL CONTROL*INPUT***** 

FIGURE 3.5. - CONCLUDED. 
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TEMPERATURE DEPENDENT MATERIAL CONTROL INPUT DATA 
FOR THE ABOVE MATERIAL 

mm!!!!!!! r mhm m !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!m}{! m urn ! 
' PLEASE * NOTE * THE * FOLLOWING : 

1. FRACTURE STRESSES FOR A GIVEN TEMPERATURE CAN BE INPUT IN 
ARBITRARY ORDER. 

2. THE DEFAULT MAXIMUM NUMBER OF TEMPERATURE SETS IS 10. 

3. THE DEFAULT MAXIMUM NUMBER OF FRACTURE SPECIMENS PER TEMPERATURE IS 
150. 

4. REGARDLESS OF THE FRACTURE ORIGIN LOCATION, THE FRACTURE STRESS 
INPUT VALUE IS THE EXTREME FIBER STRESS WITHIN THE INNER LOAD SPAN 
OF THE MOR BAR. 

t I I I H T I 1 ! II f I I t II t II I H J ! V ! t I I ; J J J j ; I H ( { M I [I j ! I } ! I j ! ( lj f M ! f 11 f f I M I I 


TDEG : TEMPERATURE OF THIS SET 


00070.0000 


NUT : NUMBER OF FRACTURE SPECIMENS AT THIS TEMPERATURE 


*015* 


MOR 

: S-URFACE, 

v- 

-OLUME, OR U-NKNOWN FLAW AND RESPECTIVE STRESS 

vw 

0.457500E+05 


0.461000E+05 

0 . 481000E+05 

vw 

0.481250E+05 


0. 491250E+05 

0.491880E+05 

vw 

0.495000E+05 


0.496250E+05 

0 . 496500E+05 

vsv 

0.497500E+05 


0 . 498500E+05 

0 . 498900E+05 

SSU 

0 . 506250E+05 


0.516250E+05 

0 . 522500E+05 


j. 


j. 



END OF DATA 

FOR THE ABOVE TEMPERATURE 

TDEG 

: TEMPERATURE 

OF THIS SET 


00500.0000 




NUT 

: NUMBER OF FRACTURE SPECIMENS AT THIS TEMPERATURE 

*015* 

i> 




MOR 

: S-URFACE, 

V- 

-OLUME, OR U-NKNOWN 

FLAW AND RESPECTIVE STRESS 

vw 

0.407500E+05 


0.411000E+05 

0.431000E+05 

vw 

0 .431250E+05 


0. 441250E+05 

0.441880E+05 

vw 

0 . 445000E+05 


0.446250E+05 

0.446500E+05 

vw 

0.447500E+05 


0.448500E+05 

0.448900E+05 

wv 

0.456250E+05 


0 . 466250E+05 

0 .472500E+05 


END OF DATA 

FOR THE ABOVE TEMPERATURE 

TDEG 

: TEMPERATURE 

OF THIS SET 



01000.0000 


FIGURE 3.4. - TEMPERATURE-DEPENDENT FRACTURE DATA FOR VOLUME FLAW ANALYSIS (CENSORED DATA OPTION). 
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NUT : NUMBER GF FRACTURE SPECIMENS AT THIS TEMPERATURE 


*015* 


MGR : S-URFACE, V-OLUME, OR U-NKNOWN FLAW AND RESPECTIVE STRESS 

— * * * * 

UVV 0. 357500E+05 0.361000E+05 0.381000E+05 

VW 0. 381250E+05 0.391250E+05 0.391880E+05 

VVV 0. 395000E+05 0.396250E+05 0.396500E+05 

WS 0.397500E+05 0.398500E+05 0.398900E+05 

VSS 0 . 406250E+05 0.416250E*05 0.422500E+05 

— * * * * 

END OF DATA FOR THE ABOVE TEMPERATURE 

******************** **************************************************** 
$ENDT : END OF DATA FOR THE ABOVE MATERIAL 

************************************************************************ 

FIGURE 3.<l. - CONCLUDED. 


21 



MATERIAL CONTROL INPUT 


TITLE : MATERIAL TITLE (ECHOED IN PC-CARES OUTPUT) 


SAMPLE INPUT FOR PC-CARES USERS MANUAL 


MATID : MATERIAL IDENTIFICATION NUMBER 

(NO DEFAULT) 

*0000300+ 


ID1 : CONTROL INDEX FOR EXPERIMENTAL DATA 

(NO DEFAULT) 

+2* 1 : UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

2 : FOUR-POINT BEND TEST DATA 

3 : DIRECT INPUT OF THE WEIBULL PARAMETERS , M AND SP 

(SHAPE PARAMETER AND SCALE PARAMETER) 

4 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

5 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

FOUR-POINT BEND TEST DATA 

ID4 : CONTROL INDEX FOR VOLUME OR SURFACE FLAW ANALYSIS 

(NO DEFAULT) 

*2* 1 : VOLUME 

2 : SURFACE 

ID2S : CONTROL INDEX FOR SURFACE FRACTURE CRITERION 

(NO DEFAULT) 

*5* 1 : NORMAL STRESS FRACTURE CRITERION 

(SHEAR-INSENSITIVE CRACK) 

3 : COPLANAR STRAIN ENERGY RELEASE RATE CRITERION 

(G SUB T) 

4 : WEIBULL PIA MODEL 

5 : S HETTY ’ S SEMI -EMPIRICAL CRITERION 

ID3S : CONTROL INDEX FOR SHAPE OF SURFACE CRACKS 

(NO DEFAULT) 

*4* 1 : GRIFFITH CRACK 

(ASSOCIATED WITH STRAIN ENERGY RELEASE RATE CRIT.) 

(ASSOCIATED WITH SHETTY’S SEMI-EMPIRICAL CRITERION) 

3 : GRIFFITH NOTCH 

(ASSOCIATED WITH STRAIN ENERGY RELEASE RATE CRIT.) 
(ASSOCIATED WITH SHETTY’S SEMI -EMPIRICAL CRITERION) 

4 : SEMICIRCULAR CRACK 

(ASSOCIATED WITH SHETTY’S SEMI-EMPIRICAL CRITERION) 

IKBAT : CONTROL INDEX FOR METHOD OF CALCULATING BATDORF CRACK 

DENSITY COEFFICIENT (K SUB B) FROM TEST DATA 

*0* (DEFAULT: IKBAT = 0) 

0 : SHEAR-INSENSITIVE METHOD (MODE I FRACTURE ASSUMED) 

1 : SHEAR-SENSITIVE METHOD (FRACTURE ASSUMED TO OCCUR 
ACCORDING TO THE FRACTURE CRITERION AND CRACK SHAPE 
SELECTED BY THE ID2 AND ID3 INDICES) 

PR : POISSON’S RATIO 

(DEFAULT: PR = 0.25) 

00000.2500 


FIGURE 3.5. - PC-CARES MATERIAL CONTROL INPUT FOR SURFACE FLAW ANALYSIS. 
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c : CONSTANT FOR SIIETTY’S SEMI-EMPIRICAL MIXED-MODE FRACTURE 

CRITERION (KI/KIC) + (KII/ (C*KIC)) +*2 = 1 

00000.8000 OBSERVED VALUES RANGE FROM 0.8 TO 2. (REF. D.K. SIIETTY) 

NOTE: AS C APPROACHES INFINITY, PREDICTED FAILURE 

PROBABILITIES APPROACH NORMAL STRESS CRITERION VALUES 
(DEFAULT: C = 1.0) 

MLORLE : CONTROL INDEX FOR METHOD OF CALCULATING WEIBULL 

PARAMETERS FROM THE EXPERIMENTAL FRACTURE DATA 

♦1* (DEFAULT: MLORLE = 0) 

0 : MAXIMUM LIKELIHOOD 

1 : LEAST-SQUARES LINEAR REGRESSION 

DH : HEIGHT OF THE FOUR-POINT BEND BAR 

( N0 DEFAULT) 

00000.0710 


DL1 : OUTER LOAD SPAN OF THE FOUR-POINT BEND BAR 

(No DEFAULT) 

00001 . 0000 


DL2 : INNER LOAD SPAN OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

00000.5000 


DW : WIDTH OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

00000.1477 


JENDM : END OF TEMPERATURE INDEPENDENT MATERIAL CONTROL INPUT 

FIGURE 3.5. - CONCLUDED. 
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TEMPERATURE DEPENDENT MATERIAL CONTROL INPUT DATA 
FOR THE ABOVE MATERIAL 

j j i ; ;n ; n n i j in!!;!!!; j n! j H!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! 
PLEASE NOTE THE FOLLOWING: 

1. FRACTURE STRESSES FOR A GIVEN TEMPERATURE CAN BE INPUT IN 
ARBITRARY ORDER. 

2. THE DEFAULT MAXIMUM NUMBER OF TEMPERATURE SETS IS 10. 

3. THE DEFAULT MAXIMUM NUMBER OF FRACTURE SPECIMENS PER TEMPERATURE IS 
150. 

4 REGARDLESS OF THE FRACTURE ORIGIN LOCATION, THE FRACTURE STRESS 
INPUT VALUE IS THE EXTREME FIBER STRESS WITHIN THE INNER LOAD SPAN 
OF THE MOR BAR. 

!!!!!!!!! I n j 1 1 !! j !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! 

TDEG : TEMPERATURE OF THIS SET 


00070.0000 


NUT : NUMBER OF FRACTURE SPECIMENS AT THIS TEMPERATURE 


*015* 


MOR : FRACTURE STRESSES 


0 . 457500E+05 
0 . 481250E+05 
0.495000E+05 
0.497500E+05 
0 . 506250E+05 


0. 461000E+05 
0 . 491250E+05 
0 . 496250E+05 
0 . 498500E+05 
0.516250E+05 


0.481000E+05 
0.491880E+05 
0 . 496500E+05 
0 . 498900E+05 
0 . 522500E+05 


* 

END OF DATA FOR THE ABOVE TEMPERATURE 

TDEG 

: TEMPERATURE OF THIS SET 

00500 . 0000 

- 

NTJT 

: NUMBER OF FRACTURE SPECIMENS AT THIS TEMPERATURE 

*015* 


MOR 

* 

: FRACTURE STRESSES 


0.407500E+05 
0.431250E+05 
0 . 445000E+05 
0.447500E+05 
0 . 456250E+05 


0.411000E+05 
0.441250E+05 
0.446250E+05 
0 . 448500E+05 
0 . 466250E+05 


0.431000E+05 
0.441880E+05 
0 . 446500E+05 
0.448900E+05 
0.472500E+05 


♦ 


* — — — — — — — “* * 

END OF DATA FOR THE ABOVE TEMPERATURE 


TDEG : TEMPERATURE OF THIS SET 


01000.0000 


FIGURE 3.6. - TEMPERATURE-DEPENDENT FRACTURE DATA FOR SURFACE FLAW ANALYSIS (COMPLETE SAMPLE OPTION). 
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NUT : NUMBER OF FRACTURE SPECIMENS AT THIS TEMPERATURE 


+015* 


MOR : FRACTURE STRESSES 

* * + * 

0. 357500E+05 0.361000E+05 0.381000E+05 

0.381250E+05 0.391250E+05 0.391880E+05 

0. 395000E+05 0.396250E+05 0.396500E+05 

0 . 397500E+05 0.398500E+05 0.398900E+05 

0 . 406250E+05 0.416250E+05 0.422500E+05 

* + + * 

END OF DATA FOR THE ABOVE TEMPERATURE 

***** + ** + ***** + + + * + + +!***** + + ****** + * + + * + * + t*********t***+ + + ++ + ** + ** t + ** t 

$ENDT : END OF DATA FOR THE ABOVE MATERIAL 

****** + ************ + ** + *** + + + + ** + **** + * + ** + ****** ^.fc^******^^*^^^^ 

FIGURE 3.6. - CONCLUDED. 
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TEMPERATURE DEPENDENT MATERIAL CONTROL INPUT DATA 
FOR THE ADOVE MATERIAL 

!!!!!!!!!! | M !!!!!!!!!! !!! !!!!!!!!!!!!!!! M !!!!!!!!!!!!!!!!!!!!!!!!!!!! ! 
PLEASE NOTE THE FOLLOWING: 

1 THE DEFAULT MAXIMUM NUMBER OF TEMPERATURE SETS IS 10. 

!!!!!!!!!!!!!! 1 1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! M 

TDEG : TEMPERATURE OF THIS SET 


00070 . 0000 


PARAM : WE1BULL MODULUS (SHAPE PARAMETER) AND SCALE PARAMETER 

+-WEIBULL MODULUS- ‘-SCALE PARAMETER-* 

0.765000E+01 0.878910E+05 

* * * 

END OF DATA FOR THE ABOVE TEMPERATURE 


END OF DATA FOR THE ABOVE MATERIAL 


FIGURE 3./. - DIRECT INI’UI Of TEMPI RA1 URE -DEPENHEN T STATISTICAL MAI ER I AL PAR/U1EIERS. 


26 


4.0 EXECUTION OF THE PC-CARES PROGRAM 

Prior to the execution of the PC-CARES program the user must prepare the 
Initialization file (optional) and the PC-CARES input file as per the instruc- 
tions of section 3.0 PC-CARES INPUT INFORMATION. Further, if one intends to 
use a special directory other than the default directory for the PC-CARES 
output file, then that directory must be created prior to execution of the 
program. 

As can be noted from the distribution disk list In section 1.0 INTRODUC- 
TION, there are two executable files included, PCCARES.EXE and PCCARESB.EXE. 

The former execution file is the PC-CARES code linked for DOS mode (or the 
OS/2 DOS compatabi 1 i ty box) execution only. The latter code is the PC-CARES 
bound execution file, which means that this file will run in either DOS (real) 
mode or OS/2 (protected) mode. The reason for including the real mode 
execution file when the bound execution file will run in either mode is that 
the real mode execution file is smaller and will allow the user to have more 
array space for fracture test data which may make a significant difference 
given the DOS 640K memory limitation. For users running DOS only, it is recom- 
mended that you use the PCCARES.EXE execution file. 

If you do have OS/2 version 1.1 or higher you can take advantage of the 
larger address space of protect mode and multiasking by executing PCCARESB.EXE 
from the OS/2 command prompt. 

Execution of PC-CARES is straight forward. After setting up the necessary 
input files and setting the default directory to the location of the execution 
files, simply type PCCARES (or PCCARESB) at the DOS (or OS/2) command prompt. 
Immediately the message 'EXECUTING PC-CARES...' should appear on the screen. 
Except for any warning or error messages the program will run to completion 
without any user interaction. When the program finishes the message 

PC-CARES EXECUTION COMPLETE 
Stop - Program Terminated 

will appear on the screen followed by the command prompt. You then may view 
the output of the program by printing the file, using the TYPE command, or 
simply using your screen editor to view it. 

The program will generate two types of error messages: (1) Program Con- 

trol Errors and (2) Data Control Errors. The former will always appear on the 
screen and Immediately halt the program. These errors are generated by the 
inability of the program to either open the necessary files or allocate the 
array space dynamically; these are essentially errors associated with the Ini- 
tialization file. In the latter case the program will appear to terminate 
normally but upon viewing the output file an error message will be found. The 
Data Control Errors result from missing or inconsistent input from either the 
Master Control Input or the Material Control Input of the PC-CARES input file. 
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5.0 PC-CARES OUTPUT INFORMATION 


The first part of the PC-CARES output is an echo of the choices selected 
(or default values) from the Master Control Input. The PRINTA subroutine 
echoes these data. 

Then PC-CARES proceeds to the PRINTB subroutine to echo the user inputs 
for each section of the Material Control Input. The results of the analysis of 
the data from the Material Control Input are output in the PRINTP subroutine. 

If statistical material parameters are directly input, then output pertaining 
to calculated values of the normalized Batdorf crack density coefficient will 
follow. If statistical material parameters are determined from experimental 
fracture data, then the output will identify the method of solution, the con- 
trol Index used for experimental data, the number of specimens in each batch, 
and the temperature of each test. In addition, the output echoes the input 
values of all specimen fracture stresses with proper failure mode identifica- 
tion if the user has set IPRINT = 1 in the Master Control Input. Any data 
value that deviates greatly from the rest of the sample is detected as an out- 
lier, and its corresponding significance level is printed. Three levels of 
significance are available for outliers: 1 , 5, or 10 percent. The lower the 

significance level, the more extreme is the deviation of the data point from 
the rest of the distribution. A 1-percent significance level indicates that 
there is a l-in-100 chance that the data point is actually a member of the same 
population as the other data, assuming a normal distribution. 

Next, the biased and the unbiased value of the shape parameter, the spec- 
imen characteristic strength, the upper and lower bound values at 90-percent 
confidence level for both the shape parameter and the specimen characteristic 
strength, the specimen Weibull mean value, and the corresponding standard devi- 
ation are printed for each specified temperature. For censored statistics, 
these values are generated first for the volume flaw analysis and subsequently 
for the surface flaw analysis. Not all of this information is available for 
all methods of material parameter estimation, and section 2.0 PROGRAM CAPABIL- 
ITY AND DESCRIPTION of this manual should be consulted for further information. 

The Kolmogorov-Smi rnov (K-S) goodness-of-f i t test is done for each data 
point, and the corresponding K-S statistics (D+ and D-) and significance 
level are listed. Similarly, the K-S statistic D for the overall population 
is printed along with the significance level. This overall statistic is the 
absolute maximum of individual specimen data D+ and D- factors. For the 
Anderson-Darling (A-D) goodness-of-f i t test, the A-D statistic A 2 is deter- 
mined for the overall population and its associated significance level is 
printed. The lower the significance level, the worse is the fit of the exper- 
imental data to be proposed distribution. For these tests, a 1-percent level 
of significance indicates that there is a l-in-100 chance that the specimen 
fracture data is from the estimated distribution. 

The next table generated by PC-CARES from the PRINTP subroutine contains 
data to construct Kanofsky-Srinivasan 90-percent confidence bands about the 
Weibull distribution. The table includes fracture stress data, the correspond- 
ing Weibull probability of failure values, the 90-percent upper and lower con- 
fidence band values about the Weibull line, and the median rank value for each 
data point. These statistical quantities are calculated with either tabular 
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values or approximating polynomial functions. Experimental fracture data lying 
outside these bands are an indication of poor fit to the Weibull distribution. 

The last table from the PRINTP subroutine summarizes the material param- 
eters listed as a function of temperature. These include the biased Weibull 
modulus, the Batdorf crack density coefficient, and the material Weibull scale 
parameter (unit volume or unit area characteristic strength). The values given 
correspond to the experimental temperatures input. Information on the selected 
fracture criterion and crack shape is printed for shear-sensitive fracture 
models. Crack shape is not required for the shear-insensitive fracture cri- 
terion or for the PIA model, and it need not be identified for those cases. 


6.0 EXAMPLE PROBLEM 

The following example and discussion of estimating the statistical mate- 
rial parameters was, like the THEORY section in the appendix, obtained from the 
original CARES manual. The CARES input file for this problem is included on 
the distribution disk in the file called PCCARES.DAT to allow the user to run 
the analysis in order to familiarize himself (herself) with PC-CARES operation. 
Following the text below is a copy of both the input file, PCCARES.DAT and the 
PC-CARES generated output. 

Example - Statistical Material Parameter Estimation 

To validate the methods used to estimate statistical material parameters, 
we compared results from the fracture of four-point bend bars broken at NASA 
Lewis and analyzed by CARES with results independently obtained by Bruckner- 
Foit and Munz (ref. 15) for the International Energy Agency (IEA) Annex II, 
Subtask 4 (ref. 16). The IEA Annex II agreement is focused on cooperative 
research and development among the United States, West Germany, and Sweden in 
the areas of structural ceramics. Subtask 4 of the agreement addresses mechan- 
ical property measurement methods with initial research concentrating solely on 
four-point flexure testing. Three different materials were analyzed, namely a 
hot isostatic pressed (HIPed) silicon carbide (SiC) from Elektroschmel zwerke 
Kempten (ESK), West Germany, a HIPed silicon nitride (Si 3 N 4 ) from ASEA CERAMA, 
Sweden; and a sintered silicon nitride from GTE WESGO, USA, although only 
results from the ESK and ASEA materials are discussed herein. 

In November 1986, 400 HIPed SiC flexure bars from West Germany were 
distributed by Oak Ridge National Laboratory (ORNL) (Oakridge, Tennessee) to 
the five participating U S. laboratories, including NASA Lewis. The bars were 
fractured at these laboratories and the fracture stress data sets were 
returned to ORNL as complete data without censoring four different failure 
modes. Shortly thereafter, 400 S i 3 N 4 bars from Sweden were also received by 
ORNL and subsequently distributed to the same U.S. laboratories for fracture 
testing. Again, the fracture stress data sets were returned to ORNL as 
complete samples. The number of specimens of a particular material given to 
each U.S. laboratory was 80. The specimens had cross-sectional dimensions of 
3.5 mm (0.138 in.) in width and 4.5 mm (0.177 in.) in height. The specimens 
were tested in four-point bending with an outer span of 40 mm (1.57 in.) and 
an inner span of 20 mm (0.787 in.). The nominal loading rate was 0.5 mm/min 
(0.020 in. /min), and the testing temperature was approximately 20 °C (68 °F). 
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Details of the statistical analyses of these data sets are given in refer- 
ences 15 and 16. The results of the 80 silicon carbide flexure bars tested at 
NASA Lewis, which are shown in table I, were analyzed with the CARES code to 
calculate the maximum likelihood estimates (MLE's) of the Weibull parameters. 
The Weibull parameter values from CARES, summarized in table II, match the pre- 
dictions from reference 15 reasonably well. The SiC fracture data are plotted 
in figure 6.1 along with the proposed Weibull line and the Kanof sky-Sri ni vasan 
90-percent confidence bands. Since all of the data are within the 90-percent 
bands and the goodness-of-f 1 t significance levels are high, it is concluded 
that the fracture data show good Weibull behavior. 

ASEA CERAMA HIPed Si 3 N 4 bars (ref. 16) from Sweden were also fractured at 
NASA Lewis, and subsequently, the statistical material parameters were estima- 
ted with CARES by using the maximum likelihood method. A comparison of the 
Si 3 N 4 results with those in reference 15 is also shown in table II. Agreement 
between estimates from the two sources is excellent. When the 80 ASEA silicon 
nitride bars were analyzed by the CARES code as a complete sample, the signif- 
icance levels of 54 and 35 percent from the Kolmogorov-Smi rnov and Anderson- 
Darling goodness-of-f i t tests, respectively, were relatively low, indicating a 
questionable fit to the proposed Weibull distribution. The lower significance 
level for the Anderson-Darling test indicated greater deviation occurring in 
the low strength region of the distribution. From the outlier test included 
in the CARES code analysis package, the highest strength fracture stress was 
detected to be an outlier at the 1-percent significance level. Several of the 
lower strengths were flagged as outliers at various significance levels (1, 5, 
or 10 percent). Figure 6.2 shows a Weibull plot of the data. From the figure 
it appears that the data are bimodal with an outlier point at the highest 
strength. 

Because of the observed trends, the data were re-analyzed assuming a 
censored distribution and removing the highest strength outlier point 
(af = 817.2 MPa (1.1 85x1 0^ ps i ) ) as bad data. Although it is possible that 
both failure modes were surface induced, for the sake of this example it is 
assumed that the low-strength failures were predominantly due to volume flaws 
and that the high-strength specimens fractured predominantly because of surface 
flaws. Since results from fractography of the individual specimens to identify 
the various failure modes were not available, the fracture origins had to be 
arbitrarily assigned prior to parameter estimation. Note that identifying 
individual specimen flaw origins is especially important for small sample 
sizes where a plot of the data does not yield clear trends. However, for the 
NASA Lewis Si 3 N 4 data, the sample size was large, and clear trends could be 
observed, although extra care would be required to determine if the trends were 
surface flaw or volume flaw based. From inspection of figure 6.2, we decided 
to assign the lowest nine strengths as due to volume flaws and the remainder as 
due to surface flaws. The cracks were arbitrarily assumed to be Griffith 
cracks, and the total strain energy release rate fracture criterion was used. 
This assumption was used only in the calculation of kgs and k gy. Th e <-$ 
significance level increased from 0.54 to 0.68, and the A-D significance level 
increased from 0.35 to 0.58. This improvement supports the initial assumption 
of bimodal behavior. The value of m (the superscript ~ indicates an esti- 
mated parameter) changed from 13.4 for the complete sample to m$ = 22.8 and 
my = 4.13. The value of ctq changed from 686 MPa (9.950xl0 4 psi) for the 
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complete sample to o@s = 692 MPa (1.004x10^ psi) and oq\j = 1128 MPa 
(1.636x10^ psi) for the surface and volume flaw distributions, respectively. 
Further Improvements in the goodness-of-f i t scores may be gained by correctly 
identifying the location of fracture origins. 

From equation (A. 82) the normalized Batdorf crack density coefficient for 
volume flaws Is (my + 1) = 4.13 + 1 = 5.13, and from equation (A. 72) the scale 
parameter a 0 y is 17.9 MPa <m) 3/4 . 1 3 (3.742xl0 4 psi < 1 n. )3/4. 1 3) . re- 
surface flaws the normalized Batdorf crack density coefficient Is 6.05. 
whereas a n c calculated by using equation (A. 87a) is 461.3 MPa (m)2/22.8 
(9.234xl0 4 psi (In. >2/22.8). 

For the Si 3 N 4 fracture data from NASA Lewis, we have obtained goodness-of- 
fit significance levels as high as 0.78 and 0.88 for the K-S and A-D tests, 
respectively, by assuming a particular bimodal flaw distribution. For this 

case, 13 volume flaws were assumed, and the MLE's were mg = 21.00, my = 6.79, 
oqs = 693 MPa (1.005xl0 5 psi), and o ev = 876 MPa <1.271xl0 5 psi). The 13 vol- 
ume flaws did not correspond to the 13 lowest fracture strengths. On the basis 
of these goodness-of-f i t scores, it is concluded that the data show good 
bimodal Weibull behavior. 

It should be noted from figure 6.2 that the assumed volume flaw distribu- 
tion dominates the failure response at low probabilities of failure. There- 
fore, in component design, it is essential to properly account for competing 
failure modes; otherwise nonconservative design predictions may result. 
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TABLE I. - EXTREME FIBER FRACTURE STRESSES 
OF ESK HIPed SILICON CARBIDE (SiC) BARS 


FI exure 
bar 

Strength, 

MPa 

1 

281 .2 

2 

291 .0 

3 

358.2 

4 

385.4 

5 

389.0 

6 

390.8 

7 

391.8 

8 

402.8 

9 

412.5 

10 

413.3 

11 

413.9 

12 

417.8 

13 

418.2 

14 

426.9 

15 

437.6 

16 

440.0 

17 

441.0 

18 

442.5 

19 

443.8 

20 

444.9 

21 

446.2 

22 

451.5 

23 

452.1 

24 

452.7 

25 

470.4 

26 

474.1 

27 

475.5 

28 

475.5 

29 

479.2 

30 

483.5 

31 

484.8 

32 

486.2 

33 

488.6 

34 

492.5 

35 

493.2 

36 

496.0 

37 

505.7 

38 

511.9 

39 

512.5 

40 

513.8 


Flexure 

bar 

Strength, 

MPa 

41 

516.2 

42 

519.8 

43 

527.6 

44 

530.7 

45 

530.7 

46 

545.7 

47 

548.8 

48 

552.7 

49 

559.6 

50 

562.4 

51 

563.3 

52 

566.1 

53 

566.5 

54 

570.1 

55 

572.8 

56 

575.0 

57 

576.1 

58 

580.0 

59 

582.6 

60 

588.0 

61 

588.6 

62 

591.0 

63 

591 .0 

64 

593.3 

65 

598.7 

66 

599.6 

67 

610.0 

68 

612.7 

69 

619.9 

70 

619.9 

71 

622.2 

72 

622.3 

73 

640.5 

74 

649.0 

75 

657.2 

76 

660.0 

77 

664.3 

78 

673.5 

79 

673.9 

80 

725.3 
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FRACTURE STRESS, MPa 

FIGURE 6.2, - COMPLETE SAMPLE, ASSUMED CENSORED 
SAMPLE WE I BULL DISTRIBUTIONS, AND OUTLIER IN 
ASEA CERAMA HIPped SILICON NITRIDE (S^). 
(FRACTURE STRESS DATA GENERATED AT NASA-LEWIS; 
NOT ALL DATA POINTS SHOWN.) 
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PC-CARES Templet Input File 


*********************** **************** ***************************** 
*************************************************************** ++ .^ + 

MASTER CONTROL INPUT 

TITLE : PROBLEM TITLE (ECHOED IN PC-CARES OUTPUT) 

EXAMPLE PROBLEM : STATISTICAL MATERIAL PARAMETER ESTIMATION 


NMATS : NO. OF MATERIALS FOR SURFACE FLAW ANALYSIS 

(NMATS+NMATV < 101) 

*01* (DEFAULT: NMATS = 0) 

NMATV : NO. OF MATERIALS FOR VOLUME FLAW ANALYSIS 

(NMATS+NMATV < 101) . 

*01* (DEFAULT: NMATV = 0) 

IPRINT : CONTROL INDEX FOR STRESS OUTPUT 

(DEFAULT: IPRINT = 0) 

*1* 0 : DO NOT PRINT FRACTURE DATA 

1 ; print FRACTURE DATA 

NGP : NO. OF GAUSSIAN QUADRATURE POINTS (15 OR 30) 

(DEFAULT: NGP = 15) 

*30* 


******************************** **************************************** 
SENDX : END OF MASTER CONTROL INPUT 

************************************************************************ 

MATERIAL CONTROL INPUT 

TITLE : MATERIAL TITLE (ECHOED IN PC-CARES OUTPUT) 

SI3N4 SPECIMEN DATA FROM ASEA CERAMA FOR VOLUME FLAW ANALYSIS 


MATID : MATERIAL IDENTIFICATION NO. FROM THE FINITE ELEMENT 

MATERIAL PROPERTY CARD (IF POSTPROCESSING IS NOT 

*0000001* BEING PERFORMED THIS ENTRY SHOULD BE SOME UNIQUE NO.) 

(NO DEFAULT) 


ID1 


*5* 


CONTROL INDEX FOR EXPERIMENTAL DATA 
(NO DEFAULT) 

1 : UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

2 : FOUR-POINT BEND TEST DATA 

3 : DIRECT INPUT OF THE WEIBULL PARAMETERS, M AND SP 

(SHAPE PARAMETER AND SCALE PARAMETER) 

4 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

5 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

FOUR-POINT BEND TEST DATA 
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ID4 : CONTROL INDEX FOR VOLUME OR SURFACE FLAW ANALYSIS 

(NO DEFAULT) 

*1* 1 : VOLUME 

2 : SURFACE 

ID2V : CONTROL INDEX FOR VOLUME FRACTURE CRITERION 

(NO DEFAULT) 

*3* 1 : NORMAL STRESS FRACTURE CRITERION 

(SHEAR-INSENSITIVE CRACK) 

2 : MAXIMUM TENSILE STRESS CRITERION 

3 : COPLANAR STRAIN ENERGY RELEASE RATE CRITERION 

(G SUB T) 

4 : WEIBULL PIA MODEL 

5 : SHETTY’S SEMI-EMPIRICAL CRITERION 

ID3V : CONTROL INDEX FOR SHAPE OF VOLUME CRACKS 

(NO DEFAULT) 

*1* 1 : GRIFFITH CRACK 

2 : PENNY-SHAPED CRACK 

IKBAT : CONTROL INDEX FOR METHOD OF CALCULATING BATDORF CRACK 

DENSITY COEFFICIENT (K SUB B) FROM TEST DATA 

*1* (DEFAULT: IKBAT = 0) 

0 : SHEAR-INSENSITIVE METHOD (MODE I FRACTURE ASSUMED) 

1 : SHEAR-SENSITIVE METHOD (FRACTURE ASSUMED TO OCCUR 
ACCORDING TO THE FRACTURE CRITERION AND CRACK SHAPE 
SELECTED BY THE ID2 AND ID3 INDICES) 

PR : POISSON’S RATIO 

(DEFAULT: PR = 0.25) 

00000.2500 


MLORLE : CONTROL INDEX FOR METHOD OF CALCULATING WEIBULL 

PARAMETERS FROM THE EXPERIMENTAL FRACTURE DATA 

*0* (DEFAULT: MLORLE = 0) 

0 : MAXIMUM LIKELIHOOD 

1 : LEAST-SQUARES LINEAR REGRESSION 

DH : HEIGHT OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

0000.00350 


DL1 : OUTER LOAD SPAN OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

0000.04000 


DL2 : INNER LOAD SPAN OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

0000.02000 


DW : WIDTH OF THE FOUR-POINT BEND BAR 

(NO DEFAULT) 

0000.00450 


SENDM : END OF TEMPERATURE INDEPENDENT MATERIAL CONTROL INPUT 
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TEMPERATURE DEPENDENT MATERIAL CONTROL INPUT DATA 
FOR THE ABOVE MATERIAL 

I ! j !!!!!! j !!!! j !!!!!! i t !!!!!!!!!! f !!!!!!!! m !!!!!! j !!!!! j !!! j !!!! j !!!!! ! 

PLEASE NOTE THE FOLLOWING: 

1. FRACTURE STRESSES FOR A GIVEN TEMPERATURE CAN BE INPUT IN 
ARBITRARY ORDER. 

2. THE DEFAULT MAXIMUM NUMBER OF TEMPERATURE SETS IS 10. 

3. THE DEFAULT MAXIMUM NUMBER OF FRACTURE SPECIMENS PER TEMPERATURE IS 
150. 

4. REGARDLESS OF THE FRACTURE ORIGIN LOCATION, THE FRACTURE STRESS 
INPUT VALUE IS THE EXTREME FIBER STRESS WITHIN THE INNER LOAD SPAN 
OF THE MOR BAR. 

1 1 1 1 1 1 f t f i ii ii r 1 1 i 1 1 1 1 1 1 ; ! t j j j ; f i ; ; ; ; i i ; j j t ; i f f ; ; j i r f i t i r 1 1 i i m i i i m it m 


TDEG : TEMPERATURE OF THIS SET 


00070.0000 


NUT : NUMBER OF FRACTURE SPECIMENS AT THIS TEMPERATURE 


*079* 


MOR : S-URFACE, V-OLUME, OR U-NKNOWN FLAW AND RESPECTIVE STRESS 


* 

SW 0 . 6257900000E+03 
SSS 0.691 1000000E+03 
SSS 0 . 6950100000E+03 
SSS 0. 6616800000E+03 
SSS 0 . 6707500000E+03 
SSS 0 . 6401800000E+03 
VSS 0 . 5725000000E+03 
SSS 0 . 6597400000E+03 
SSS 0 . 6621400000E+03 
SSS 0 . 7156300000E+03 
VVS 0 . 5605800000E+03 
SSS 0.6311 300000E+03 
SSS 0 . 6770800000E+03 
SSS 0 . 7156300000E+03 
VSS 0 . 5457700000E+03 
SSS 0 . 6837300000E+03 
SSS 0 . 6490400000E+03 
SSS 0 . 7044700000E+03 
SSS 0 . 6215900000E+03 
SSS 0 . 6515900000E+03 
SSV 0 . 7 100500000E+03 
SSS 0. 6688200000E+03 
SSS 0. 6093200000E+03 
SSS 0. 6210000000E+03 
SVS 0 . 6620500000E+03 
SSS 0 . 7286600000E+03 
V 0 . 5199000000E+03 


0 . 6034700000E+03 
0 . 6831 100000E+03 
0 . 6721000000E+03 
0 . 7270500000E+03 
0 . 7029700000E+03 
0 . 7045300000E+03 
0, 6251400000E+03 
0 . 6501400000E+03 
0. 6997700000E+03 
0 . 6045300000E+03 
0 . 4159200000E+03 
0 . 7254400000E+03 
0 . 7251700000E+03 
0 . 7173200000E+03 
0 . 7323300000E+03 
0 . 6144200000E+03 
0 . 6867500000E+03 
0 . 6570300000E+03 
0.71651 00000E+03 
0 . 7 0288 1 0000E+03 
0 . 6550600000E+03 
0 . 6597400000E+03 
0 . 6711600000E+03 
0 . 687 OOOOOOOE + 03 
0 . 5950400000E+03 
0 . 6223700000E+03 


0 . 5272300000E+03 
0 . 6854900000E+03 
0 . 6569200000E+03 
0 . 7259500000E+03 
0 . 6114900000E+03 
0 . 7254400000E+03 
0 . 6741300000E+03 
0.71 64300000E+03 
0 . 6645400000E+03 
0 . 6664400000E+03 
0 . 6214600000E+03 
0 . 69401 OOOOOE+03 
0.712 6000000E+03 
0 . 7029200000E+03 
0 . 6717300000E+03 
0 . 6636100000E+03 
0 . 7157700000E+03 
0 . 6641400000E+03 
0 . 7028800000E+03 
0 . 7060500000E+03 
0 . 4583600000E+03 
0 . 6426400000E+03 
0 . 6768300000E+03 
0. 6221500000E+03 
0 . 6783500000E+03 
0 . 6798500000E+03 


* + + + 

END OF DATA FOR THE ABOVE TEMPERATURE 


SENDT : END OF DATA FOR THE ABOVE MATERIAL 
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MATERIAL CONTROL INPUT 


TITLE : MATERIAL TITLE (ECHOED IN PC-CARES OUTPUT) 


SI3N4 SPECIMEN DATA FROM ASEA CERAMA FOR SURFACE FLAW ANALYSIS 


MATID : MATERIAL IDENTIFICATION NO. FROM THE FINITE ELEMENT 

MATERLAL PROPERTY CARD (IF POSTPROCESSING IS NOT 

*0000002* BEING PERFORMED THIS ENTRY SHOULD BE SOME UNIQUE NO.) 

(NO DEFAULT) 

ID 1 : CONTROL INDEX FOR EXPERIMENTAL DATA 

(NO DEFAULT) 

* 5 * 1 : UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

2 : FOUR-POINT BEND TEST DATA 

3 : DIRECT INPUT OF THE WEIBULL PARAMETERS, M AND SP 

(SHAPE PARAMETER AND SCALE PARAMETER) 

4 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

5 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

FOUR-POINT BEND TEST DATA 

ID 4 : CONTROL INDEX FOR VOLUME OR SURFACE FLAW ANALYSIS 

(NO DEFAULT) 

*2* 1 • VOLUME 

2 : SURFACE 

ID2S : CONTROL INDEX FOR SURFACE FRACTURE CRITERION 

(NO DEFAULT) 

* 3 * 1 : NORMAL STRESS FRACTURE CRITERION 

(SHEAR-INSENSITIVE CRACK) 

3 : COPLANAR STRAIN ENERGY RELEASE RATE CRITERION 

(G SUB T) 

4 : WEIBULL PIA MODEL 

5 : SHETTY’S SEMI-EMPIRICAL CRITERION 

ID3S : CONTROL INDEX FOR SHAPE OF SURFACE CRACKS 

(NO DEFAULT) 

* 1 * 1 : GRIFFITH CRACK 

(ASSOCIATED WITH STRAIN ENERGY RELEASE RATE CRIT.) 

(ASSOCIATED WITH SHETTY’S SEMI-EMPIRICAL CRITERION) 

3 : GRIFFITH NOTCH 

(ASSOCIATED WITH STRAIN ENERGY RELEASE RATE CRIT.) 
(ASSOCIATED WITH SHETTY’S SEMI-EMPIRICAL CRITERION) 

4 : SEMICIRCULAR CRACK 

(ASSOCIATED WITH SHETTY’S SEMI-EMPIRICAL CRITERION) 

IKBAT : CONTROL INDEX FOR METHOD OF CALCULATING BATDORF CRACK 

DENSITY COEFFICIENT (K SUB B) FROM TEST DATA 

* 1 * (DEFAULT: IKBAT = 0) 

0 : SHEAR-INSENSITIVE METHOD (MODE I FRACTURE ASSUMED) 

1 : SHEAR-SENSITIVE METHOD (FRACTURE ASSUMED TO OCCUR 
ACCORDING TO THE FRACTURE CRITERION AND CRACK SHAPE 
SELECTED BY THE ID2 AND ID3 INDICES) 
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PR 

: POISSON’S RATIO 

00000.2500 

(DEFAULT: PR = 0.25) 

MLORLE 

*0* 

: CONTROL INDEX FOR METHOD OF CALCULATING WEIBULL 
PARAMETERS FROM THE EXPERIMENTAL FRACTURE DATA 
(DEFAULT: MLORLE = 0) 

0 : MAXIMUM LIKELIHOOD 

1 : LEAST-SQUARES LINEAR REGRESSION 

DH 

0000.00350 

: HEIGHT OF THE FOUR-POINT BEND BAR 
(NO DEFAULT) 

DL1 

0000.04000 

: OUTER LOAD SPAN OF THE FOUR-POINT BEND BAR 
(NO DEFAULT) 

DL2 

0000.02000 

: INNER LOAD SPAN OF THE FOUR-POINT BEND BAR 
(NO DEFAULT) 

DW 

0000.00450 

: WIDTH OF THE FOUR-POINT BEND BAR 
(NO DEFAULT) 


************************************************************************ 
SENDM : END OF TEMPERATURE INDEPENDENT MATERIAL CONTROL INPUT 

************************************************************************ 

TEMPERATURE DEPENDENT MATERIAL CONTROL INPUT DATA 
FOR THE ABOVE MATERIAL 

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

' PLEASE NOTE ’ THE * FOLLOWING i 

1. FRACTURE STRESSES FOR A GIVEN TEMPERATURE CAN BE INPUT IN 
ARBITRARY ORDER . 

2. THE DEFAULT MAXIMUM NUMBER OF TEMPERATURE SETS IS 10. 

3. THE DEFAULT MAXIMUM NUMBER OF FRACTURE SPECIMENS PER TEMPERATURE IS 
150. 

4. REGARDLESS OF THE FRACTURE ORIGIN LOCATION, THE FRACTURE STRESS 
INPUT VALUE IS THE EXTREME FIBER STRESS WITHIN THE INNER LOAD SPAN 
OF THE MOR BAR. 

I M ! I ! t ! I I I II I ! I ! It I If I I I U ! It I I I It I I l I M I ! f I II I II II ! I ! It I f I I I I I I II I I If ! 


TDEG : TEMPERATURE OF THIS SET 


00070.0000 


NUT : NUMBER OF FRACTURE SPECIMENS AT THIS TEMPERATURE 


*079* 
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: S-URFACE, V-OLUME, OR U-NKNOWN FLAW AND RESPECTIVE STRESS 


MOR 

* * * * 

SW 0 . 6257900000E+03 0 . 6034700000E+03 0.5272300000E+03 

SSS 0 . 6911000000E+03 0. 6831100000E+03 0 . 6854900000E+03 

SSS 0 . 6950100000E+03 0.6721000000E+03 0.6569200000E+03 

SSS 0. 6616800000E+03 0 . 7270500000E+03 0.7259500000E+03 

SSS 0. 6707500000E+03 0 . 7029700000E+03 0. 6114900000E+03 

SSS 0 . 6401800000E+03 0. 7045300000E+03 0 . 7254400000E+03 

VSS 0 . 5725000000E+03 0. 6251400000E+03 0. 6741300000E+03 

SSS 0 . 6597400000E+03 0 . 6501400000E+03 0. 7164300000E+03 

SSS 0.6621 400000E+02 0 . 6997700000E+03 0 . 6645400000E+03 

SSS 0 . 7156300000E+03 0 . 6045300000E+03 0 . 6664400000E+03 

WS 0 . 5605800000E+03 0 . 4159200000E+03 0. 6214600000E+03 

SSS 0 . 631 1300000E+03 0 . 7254400000E+03 0 . 69401 OOOOOE+03 

SSS 0. 6770800000E+03 0.7251700000E+03 0 . 7126000000E+03 

SSS 0 . 7156300000E+03 0.7173200000E+03 0.7029200000E+03 

VSS 0. 5457700000E+03 0 . 7323300000E+03 0. 6717300000E+03 

SSS 0 . 6837300000E+03 0 . 6144200000E+03 0.6636100000E+03 

SSS 0 . 6490400000E+03 0. 6867500000E+03 0.7157700000E+03 
SSS 0.7044700000E+03 0 . 6570300000E+03 0.6641400000E+03 

SSS 0.6215 900000E+03 0.7165100000E+03 0 . 7028800000E+03 

SSS 0. 6515900000E+03 0 . 70288 10000E+03 0 . 7060500000E+03 

SSV 0 . 7100500000E+03 0 . 6550600000E+03 0 . 4583600000E+03 

SSS 0. 6688200000E+03 0. 6597400000E+03 0 . 6426400000E+03 

SSS 0 . 6093200000E+03 0. 6711600000E+03 0 . 6768300000E+03 

SSS 0 . 6210000000E+03 0 . 6870000000E+03 0. 6221500000E+03 

SVS 0 . 6620500000E+03 0.5950400000E+03 0. 6783500000E+03 

SSS 0 . 7286600000E+03 0.6223700000E+03 0. 6798500000E+03 

V 0 . 5199000000E+03 

* * * * 

END OF DATA FOR THE ABOVE TEMPERATURE 


SENDT : END OF DATA FOR THE ABOVE MATERIAL 
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PC-CARES Output File 


cccccc 
c c 


c 

c 

c 


A 

A A 
A A 
A A 
AAAAAAAAA 


RRRRRRRR 
R R 

R R 

RRRRRRRR 


0000000000000000008000000000000000000000000000000000000000000000000000000 
QQQQOOQQOQQOQQOOOOOQOQQOOOOOOOOQOOOOOOOOOOOOOOOOOOOOOQOOOOOOQOOQQOQOQOQOO 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 

OOQOQQOOOOQQQQOQQOOOOQQOOQOOOOOOOOOOOOQOOOOQQOOQOOOOOOOOOOOOOOOOOQOQOOOQO 

OOOOOQOOQQQQOOQOOOOOOOOOQOOOOOOOQOOOOQOQQOOOOOQOQOOOOOOOOOOOOOOOQOQOOOOQO 


c c 
cccccc 


R 


R 



CERAMICS ANALYSIS AND RELIABILITY EVALUATION 



00 


00 


00 

sssssss 

00 

; s 

00 

\ 

00 

sssssss 

00 

s 

00 

! S 

00 

SSSSSSS 

00 


00 


00 


00 

1UCTURES 

00 


00 


00 


* * 

* ECHO OF MASTER CONTROL INPUT * 

* * 
*************************************************** 


TITLE = EXAMPLE PROBLEM : STATISTICAL MATERIAL PARAMETER ESTIMATION 

3 = CONTROL INDEX FOR OUTPUT OF ANALYSIS (ID4A) 

1 : VOLUME FLAW ANALYSIS OUTPUT ONLY 

2 : SURFACE FLAW ANALYSIS OUTPUT ONLY 

3 : VOLUME FLAW AND SURFACE FLAW ANALYSIS OUTPUT 


1 = NUMBER OF MATERIALS FOR SURFACE FLAW ANALYSIS (NMATS) 


1 = NUMBER OF MATERIALS FOR VOLUME FLAW ANALYSIS (NMATV) 


30 = NUMBER OF GAUSSIAN QUADRATURE POINTS, EITHER 15 OR 30 (NGP) 


1 = CONTROL INDEX FOR STRESS OUTPUT (IPRINT) 

0 : DO NOT PRINT ELEMENT STRESSES AND/OR FRACTURE DATA 

1 : PRINT ELEMENT STRESSES AND/OR FRACTURE DATA 
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* * 


* ECHO OF MATERIAL CONTROL INPUT * 

* * 


TITLE = SI3N4 SPECIMEN DATA FROM ASEA CERAMA FOR VOLUME FLAW ANALYSIS 


1 = MATERIAL IDENTIFICATION NUMBER (MATID) 


1 = CONTROL INDEX FOR VOLUME OR SURFACE FLAW ANALYSIS (ID4) 

1 : VOLUME 

2 : SURFACE 


0 = CONTROL INDEX FOR METHOD OF CALCULATING WEIBULL PARAMETERS 
FROM THE EXPERIMENTAL FRACTURE DATA (MLORLE) 

0 : MAXIMUM LIKELIHOOD 

1 : LEAST-SQUARES LINEAR REGRESSION 


5 = CONTROL INDEX FOR EXPERIMENTAL DATA (ID1) 

1 : UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

2 : FOUR-POINT BEND TEST DATA 

3 : DIRECT INPUT OF THE WEIBULL PARAMETERS, M AND SP 

(SHAPE PARAMETER AND SCALE PARAMETER) 

4 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

5 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

FOUR-POINT BEND TEST DATA 


3 = CONTROL INDEX FOR VOLUME FRACTURE CRITERION (ID2V) 

1 : NORMAL STRESS FRACTURE CRITERION 

(SHEAR- INSENSITIVE CRACK) 

2 : MAXIMUM TENSILE STRESS CRITERION 

3 : COPLANAR STRAIN ENERGY RELEASE RATE CRITERION 

4 : WEIBULL PIA MODEL 

5 : SHETTY’S SEMI-EMPIRICAL CRITERION 


1 = CONTROL INDEX FOR SHAPE OF VOLUME CRACKS (ID3V) 

1 : GRIFFITH CRACK 

2 : PENNY-SHAPED CRACK 


1 = CONTROL INDEX FOR METHOD OF CALCULATING BATDORF CRACK DENSITY 
COEFFICIENT 

(K SUB B) FROM TEST DATA (IKBAT) 

0 : SHEAR-INSENSITIVE METHOD (MODE I FRACTURE ASSUMED) 

1 : SHEAR-SENSITIVE METHOD (FRACTURE ASSUMED TO OCCUR 

ACCORDING TO THE FRACTURE CRITERION AND CRACK 
SHAPE SELECTED BY THE ID2 AND ID3 INDICES) 


.2500 = POISSON’S RATIO (PR) 
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* * 

* STATISTICAL ANALYSIS OF FRACTURE SPECIMEN DATA * 

* * 
**************************************************** 


ECHO OF SPECIMEN INPUT DATA, IN ASCENDING ORDER OF FRACTURE STRESS 


.4000E-01 = OUTER LOAD SPAN OF FOUR-POINT BEND BAR 
. 2000E-01 = INNER LOAD SPAN 
. 3500E-02 = DEPTH OF SPECIMEN 
. 4500E-02 = WIDTH OF SPECIMEN 


79 = NUMBER OF SPECIMENS IN BATCH 70.000 = TEMPERATURE OF BATCH 

"S"URFACE OR "VOLUME OR "U"NKNOWN FLAW ORIGIN AND RESPECTIVE FAILURE STRESS 


vw 

.4159E+03 

.4584E+03 

. 5199E+03 

vw 

.5272E+03 

. 5458E+03 

. 5606E+03 

vw 

. 5725E+03 

. 5950E+03 

. 6035E+03 

sss 

. 6045E+03 

. 6093E+03 

. 6115E+03 

sss 

. 6144E+03 

. 6210E+03 

. 6215E+03 

sss 

•6216E+03 

. 6222E+03 

. 6224E+03 

sss 

. 6251E+03 

. 6258E+03 

. 6311E+03 

sss 

• 6402E+03 

. 6426E+03 

. 6490E+03 

sss 

. 6501E+03 

. 6516E+03 

. 6551E+03 

sss 

. 6569E+03 

. 6570E+03 

. 6597E+03 

sss 

. 6597E+03 

. 6617E+03 

. 6621E+03 

sss 

. 6621E+03 

. 6636E+03 

. 6641E+03 

sss 

.6645E+03 

. 6664E+03 

. 6688E+03 

sss 

. 6708E+03 

. 6712E+03 

. 6717E+03 

sss 

. 6721E+03 

. 6741E+03 

. 6768E+03 

sss 

. 6771E+03 

. 6784E+03 

. 6799E+03 

sss 

. 6831E+03 

. 6837E+03 

. 6855E+03 

sss 

. 6867E+03 

. 6870E+03 

. 6911E+03 

sss 

. 6940E+03 

. 6950E+03 

. 6998E+03 

sss 

.7029E+03 

.7029E+03 

. 7029E+03 

sss 

.7030E+03 

* 7045E+03 

. 7045E+03 

sss 

.7061E+03 

.7101E+03 

. 7126E+03 

sss 

.7156E+03 

. 7156E+03 

. 7158E+03 

sss 

♦7164E+03 

.7165E+03 

. 7173E+03 

sss 

.7252E+03 

.7254E+03 

. 7254E+03 

sss 

■7260E+03 

. 7271E+03 

. 7287E+03 

s 

. 7323E+03 
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— STEFANSKY OUTLIER TEST OF SPECIMEN FRACTURE STRESSES — 


RESULTS FROM THE STEFANSKY OUTLIER TEST FOR TEMP. = 70.0000 

FAILURE STRESS 


.4159E+03 
.4584E+03 
.5199E+03 
. 5272E+03 
. 5458E+03 


DEVIATES FROM THE MAIN 
SIGNIFICANCE LEVEL 
DEVIATES FROM THE MAIN 
SIGNIFICANCE LEVEL 
DEVIATES FROM THE MAIN 
SIGNIFICANCE LEVEL 
DEVIATES FROM THE MAIN 
SIGNIFICANCE LEVEL 
DEVIATES FROM THE MAIN 
SIGNIFICANCE LEVEL 


TREND 

OF 

THE 

DATA 

AT 

THE 

1% 

TREND 

OF 

THE 

DATA 

AT 

THE 

1 % 

TREND 

OF 

THE 

DATA 

AT 

THE 

1 % 

TREND 

OF 

THE 

DATA 

AT 

THE 

5% 

TREND 

OF 

THE 

DATA AT 

THE 

10% 


DEVIATION FROM THE MAIN TREND OF THE DATA MAY 
INDICATE BAD VALUES. MULTIPLE DEVIATIONS FROM THE SAME 
REGION OF THE DISTRIBUTION INDICATE THAT EITHER A CONCURRENT 
OR A PARTIALLY CONCURRENT FLAW POPULATION HAS BEEN DETECTED 
(NOTE THAT A CONCURRENT FLAW POPULATION MAY BE PRESENT BUT NOT 
BE DETECTED BY THE OUTLIER TEST) . DEVIATIONS OCCURRING 
IN THE SAME REGION OF THE DISTRIBUTION WITH ALL THREE 
SIGNIFICANCE LEVELS (1%, 5% AND 10%) PRESENT INDICATE A 
CONCURRENT FLAW POPULATION. DEVIATIONS SHOULD BE EXAMINED AND 
TREATED ACCORDINGLY (I.E. IGNORE, CENSOR, ADJUST OR ELIMINATE 
STRESS) . JUDGEMENT OF ACTION TAKEN CAN BE DETERMINED FROM THE 
GOODNESS- OF-FIT TESTS. 

THE OUTLIER TEST IS NO SUBSTITUTE FOR GRAPHICAL EXAMINATION!!! 


— TEMP. DEP. WEIBULL MODULUS AND CHARACTERISTIC STRENGTH WITH 90% — 
CONFIDENCE BOUNDS DETERMINED BY MAXIMUM LIKELIHOOD ANALYSIS 

NOTE: 90% CONFIDENCE BOUNDS ON PARAMETERS DETERMINED FROM COMPETING FAILURE 
MODES (CENSORED DATA) ARE APPROXIMATE. 

FOR CENSORED DATA THE UNBIASED VALUE OF THE PARAMETER "M" IS NOT 
GIVEN. 


FOR 

SAMPLE SIZES 

LESS THAN 

4, CONFIDENCE 

LIMITS ARE 

NOT GIVEN. 


TEMP. 

M BIASED 

M UNBIASED 

UP M 

LOW M 


CHAR. STR. 

UP C.S. 

LOW C.S. 

MEAN 

STD. DEV. 

VOLUME 

70.0000 

.4130E+01 


. 5666E+01 

. 2179E+01 


. 1128E+04 

. 1342E+04 

. 9542E+03 

. 1024E+04 

.2791E+03 

SURFACE 

70.0000 

. 2281E+02 


.2619E+02 

. 1920E+02 


. 6917E+03 

. 6981E+03 

. 6853E+03 

. 6755E+03 

.3684E+02 
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STATISTICS FROM THE GOODNESS-OF-FIT TESTS FOR TEMP. 


70.0000 


KOLMOGOROV-SMIRNOV TEST 
ORDER FRAC. STR. WBIB. 

1 . 4159E+03 

2 . 4584E+03 

3 . 5199E+03 

4 . 5272E+03 

5 . 5458E+03 

6 . 5606E+03 

7 . 5725E+03 

8 . 5950E+03 

9 . 6035E+03 

10 . 6045E+03 

11 . 6093E+03 

12 . 6115E+03 

13 . 6144E+03 

14 . 6210E+03 

15 . 6215E+03 

16 . 6216E+03 

17 . 6222E+03 

18 . 6224E+03 

19 . 6251E+03 

20 . 6258E+03 

21 . 6311E+03 

22 . 6402E+03 

23 . 6426E+03 

24 . 6490E+03 

25 . 6501E+03 

26 . 6516E+03 

27 . 6551E+03 

28 . 6569E+03 

29 . 6570E+03 

30 . 6597E+03 

31 . 6597E+03 

32 . 6617E+03 

33 . 6621E+03 

34 . 6621E+03 

35 .6636E+03 

36 . 6641E+03 

37 . 6645E+03 

38 . 6664E+03 

39 . 6688E+03 

40 . 6708E+03 

41 . 6712E+03 

42 . 6717E+03 

43 . 6721E+03 

44 . 6741E+03 

45 . 6768E+03 

46 . 6771E+03 

47 . 6784E+03 

48 . 6799E+03 

49 . 6831E+03 

50 . 6837E+03 

51 . 6855E+03 

52 . 6867E+03 


PROB. OF FAIL. D+ FACTOR 


.0161 

-.0035 

.0240 

.0013 

.0414 

-.0034 

.0443 

.0064 

.0529 

.0104 

.0620 

.0140 

.0714 

.0172 

.0983 

.0029 

.1130 

.0009 

.1151 

.0114 

.1254 • 

.0139 

.1305 

.0214 

.1379 

.0267 

.1567 

.0205 

.1581 

.0317 

.1585 

.0440 

.1603 

.0549 

.1610 

.0668 

.1704 

.0701 

.1726 

.0805 

.1930 

.0728 

.2346 

.0439 

.2476 

.0435 

.2853 

.0185 

.2924 

.0241 

.3020 

.0271 

.3262 

.0156 

3399 

.0146 

3407 

.0264 

3616 

.0181 

3616 

.0308 

3773 

.0277 

3804 

.0374 

3811 

.0493 

3935 

.0496 

3980 

.0577 

4015 

.0669 

4182 

.0628 

4398 

.0538 

4580 

.0484 

4619 

.0571 

4674 

.0643 

4709 

.0734 

4909 

.0661 

5182 

.0515 

5207 

.0616 

5338 

.0611 

5495 

.0581 

5842 

.0360 

5909 

.0420 

6100 

.0356 

6237 

.0345 


FACTOR 

SIGNIF. LEVEL 

.0161 

99 . 0000 

.0114 

99.0000 

.0161 

99.0000 

.0063 

99.0000 

.0023 

99.0000 

-.0013 

99.0000 

- . 0045 

99 . 0000 

.0097 

99.0000 

.0118 

99.0000 

.0012 

99.0000 

-.0012 

99.0000 

-.0088 

99.0000 

-.0140 

99.0000 

-.0079 

99.0000 

-.0191 

99.0000 

-.0313 

99.0000 

-.0422 

97.1339 

-.0541 

87 . 2422 

-.0575 

83.1735 

-.0679 

68.5012 

-.0602 

79.6191 

-.0312 

99.0000 

-.0309 

99 . 0000 

-.0058 

99.0000 

-.0114 

99.0000 

-.0145 

99.0000 

-.0030 

99.0000 

-.0019 

99.0000 

-.0137 

99.0000 

-.0055 

99.0000 

-.0181 

99.0000 

-.0151 

99.0000 

-.0247 

99.0000 

-.0366 

99.0000 

-.0369 

99.0000 

-.0450 

95.5201 

-.0542 

87.1359 

-.0502 

91.3909 

-.0412 

97.5998 

-.0357 

99.0000 

-.0445 

95.8855 

-.0516 

89.9719 

-.0607 

78.8750 

-.0534 

88.0606 

-.0388 

98.4919 

-.0489 

92.5627 

-.0484 

92.9632 

-.0454 

95.2717 

-.0234 

99.0000 

-.0293 

99.0000 

-.0229 

99.0000 

-.0219 

99.0000 
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53 

. 6870E+03 

.6264 

54 

. 6911E+03 

.6712 

55 

. 6940E+03 

.7027 

56 

. 6950E+03 

.7135 

57 

. 6998E+03 

.7635 

58 

. 7029E+03 

.7946 

59 

. 7029E+03 

.7946 

60 

. 7029E+03 

.7950 

61 

. 7030E+03 

.7955 

62 

. 7045E+03 

.8099 

63 

. 7045E+03 

.8105 

64 

.7061E+03 

.8247 

65 

.7101E+03 

.8598 

66 

. 7126E+03 

.8801 

67 

. 7156E+03 

.9022 

68 

. 7156E+03 

.9022 

69 

. 7158E+03 

.9031 

70 

.7164E+03 

.9076 

71 

.7165E+03 

.9081 

72 

. 7173E+03 

.9134 

73 

.7252E+03 

.9549 

74 

.7254E+03 

.9560 

75 

-7254E+03 

.9560 

76 

.7260E+03 

.9581 

77 

.7271E+03 

.9624 

78 

.7287E+03 

.9680 

79 

. 7323E+03 

.9786 


.0445 

-.0318 

99.0000 

.0124 

.0003 

99.0000 

- . 0065 

.0192 

99.0000 

- . 0046 

.0173 

99 . 0000 

-.0419 

.0546 

97.2613 

-.0604 

.0731 

79.2561 

-.0478 

.0604 

93.5064 

-.0355 

.0482 

99.0000 

-.0233 

.0360 

99.0000 

-.0251 

.0378 

99.0000 

-.0130 

.0257 

99.0000 

-.0146 

.0272 

99.0000 

-.0370 

.0496 

98.9973 

- . 0447 

.0574 

95.7363 

-.0541 

.0667 

87.3420 

-.0414 

.0541 

97 . 5075 

-.0297 

.0424 

99 . 0000 

-.0215 

.0341 

99.0000 

-.0094 

.0220 

99; 0000 

- . 0020 

.0146 

99.0000 

-.0309 

.0435 

99 . 0000 

-.0193 

.0320 

99 . 0000 

-.0067 

.0193 

99.0000 

.0039 

.0087 

99 . 0000 

.0123 

.0003 

99.0000 

.0193 

-.0067 

99.0000 

.0214 

- . 0088 

99 . 0000 


KOLMOGOROV-SMIRNOV TEST YIELDS STATISTIC D = MAX (D+.D-) = 0805 

WITH AN ASSOCIATED SIGNIFICANCE LEVEL OF 68.5% 


ANDERSON-DARLING TEST YIELDS STATISTIC A**2 = .6725 

WITH AN ASSOCIATED SIGNIFICANCE LEVEL OF 58.3% 


KANOFSKY-SRINIVASAN 90% CONFIDENCE BANDS ABOUT THE WEIBULL DISTRIBUTION FOR 

TEMP . = 70 . 0000 


THE KANOFSKY-SRINIVASAN FACTOR FOR THIS DISTRIBUTION IS .0836 FOR A SAMPLE 

SIZE OF 79 


ORD. 

FRAC. STR. WE IB. 

PROB. OF FAIL 

1 

.4159E+03 

.0161 

2 

.4584E+03 

.0240 

3 

.5199E+03 

.0414 

4 

.5272E+03 

.0443 

5 

.5458E+03 

.0529 

6 

.5606E+03 

.0620 

7 

.5725E+03 

.0714 

8 

. 5950E+03 

.0983 

9 

. 6035E+03 

.1130 

10 

. 6045E+03 

.1151 

11 

.6093E+03 

.1254 

12 

.6115E+03 

.1305 

13 

.6144E+03 

.1379 


UPP . CONF . BAND 

WED. RANK 

LOW. CONF. BAND 

.0997 

.0088 

.0000 

.1077 

.0214 

.0000 

.1250 

.0340 

.0000 

.1279 

.0466 

.0000 

.1365 

.0592 

.0000 

.1456 

.0718 

.0000 

.1551 

.0844 

.0000 

.1820 

.0970 

.0147 

.1967 

.1096 

.0294 

.1988 

.1222 

.0315 

.2090 

.1348 

.0418 

.2141 

.1474 

.0469 

.2215 

.1599 

.0542 
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14 

. 6210E+03 

.1567 

.2403 

.1725 

.0731 

15 

. 6215E+03 

.1581 

.2418 

.1851 

.0745 

16 

. 6216E+03 

.1585 

.2422 

.1977 

.0749 

17 

. 6222E+03 

.1603 

.2440 

.2103 

.0767 

18 

. 6224E+03 

.1610 

.2447 

.2229 

.0774 

19 

. 6251E+03 

.1704 

.2540 

.2355 

,0867 

20 

. 6258E+03 

.1726 

.2563 

.2481 

.0890 

21 

. 6311E+03 

.1930 

.2766 

.2607 

.1094 

22 

. 6402E+03 

.2346 

.3182 

.2733 

.1510 

23 

. 6426E+03 

.2476 

.3313 

.2859 

.1640 

24 

. 6490E+03 

.2853 

.3690 

.2985 

.2017 

25 

. 6501E+03 

.2924 

.3760 

.3111 

.2088 

26 

. 6516E+03 

.3020 

.3856 

.3237 

.2183 

27 

. 6551E+03 

.3262 

.4098 

.3363 

.2425 

28 

. 6569E+03 

.3399 

.4235 

.3489 

.2562 

29 

. 6570E+03 

.3407 

.4243 

.3615 

.2571 

30 

. 6597E+03 

.3616 

.4453 

.3741 

.2780 

31 

. 6597E+03 

.3616 

.4453 

.3866 

.2780 

32 

. 6617E+03 

.3773 

.4609 

.3992 

' . 2937 

33 

. 6621E+03 

.3804 

.4640 

.4118 

.2967 

34 

. 6621E+03 

.3811 

.4647 

.4244 

.2975 

35 

. 6636E+03 

.3935 

.4771 

.4370 

.3098 

36 

. 6641E+03 

.3980 

.4816 

.4496 

.3144 

37 

. 6645E+03 

.4015 

.4851 

.4622 

.3178 

38 

. 6664E+03 

.4182 

.5018 

.4748 

.3345 

39 

. 6688E+03 

.4398 

.5234 

.4874 

.3562 

40 

. 6708E+03 

.4580 

.5416 

.5000 

.3743 

41 

. 6712E+03 

.4619 

.5455 

.5126 

.3782 

42 

. 6717E+03 

.4674 

.5510 

.5252 

.3837 

43 

.6721E+03 

.4709 

.5546 

.5378 

.3873 

44 

.6741E+03 

.4909 

.5745 

.5504 

.4073 

45 

. 6768E+03 

.5182 

.6018 

.5630 

.4345 

46 

. 6771E+03 

.5207 

.6043 

.5756 

.4371 

47 

. 6784E+03 

.5338 

.6175 

.5882 

.4502 

48 

.6799E+03 

.5495 

.6332 

.6008 

.4659 

49 

. 6831E+03 

.5842 

.6679 

.6134 

.5006 

50 

. 6837E+03 

.5909 

.6745 

.6259 

.5073 

51 

. 6855E+03 

.6100 

.6936 

.6385 

.5263 

52 

.6867E+03 

.6237 

.7073 

.6511 

.5401 

53 

.6870E+03 

.6264 

.7100 

.6637 

.5428 

54 

.6911E+03 

.6712 

.7548 

.6763 

.5875 

55 

. 6940E+03 

.7027 

.7864 

.6889 

.6191 

56 

. 6950E+03 

.7135 

.7971 

.7015 

.6298 

57 

. 6998E+03 

.7635 

.8471 

.7141 

.6798 

58 

.7029E+03 

.7946 

.8782 

.7267 

.7110 

59 

.7029E+03 

.7946 

.8782 

.7393 

.7110 

60 

.7029E+03 

.7950 

.8786 

.7519 

.7114 

61 

.7030E+03 

.7955 

.8791 

.7645 

.7119 

62 

.7045E+03 

.8099 

.8936 

.7771 

.7263 

63 

.7045E+03 

.8105 

.8941 

.7897 

.7269 

64 

.7061E+03 

.8247 

.9083 

.8023 

.7411 

65 

.7101E+03 

.8598 

.9434 

.8149 

.7761 

66 

.7126E+03 

.8801 

.9638 

.8275 

.7965 

67 

.7156E+03 

.9022 

.9858 

.8401 

.8185 

68 

.7156E+03 

.9022 

.9858 

.8526 

.8185 

69 

.7158E+03 

.9031 

.9867 

.8652 

.8195 

70 

.7164E+03 

.9076 

.9912 

.8778 

.8239 

71 

.7165E+03 

.9081 

.9917 

.8904 

.8245 

72 

.7173E+03 

.9134 

.9970 

.9030 

.8297 

73 

.7252E+03 

.9549 

1.0000 

.9156 

.8713 
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74 

.7254E+03 

.9560 

75 

.7254E+03 

.9560 

76 

. 7260E+03 

.9581 

77 

.7271E+03 

.9624 

78 

.7287E+03 

.9680 

79 

. 7323E+03 

.9786 


1.0000 

.9282 

.8724 

1 . 0000 

.9408 

.8724 

1.0000 

.9534 

.8745 

1 . 0000 

.9660 

.8787 

1 . 0000 

.9786 

.8844 

1 . 0000 

.9912 

.8950 


* * 

* VOLUME FLAW PARAMETER ANALYSIS * 

* * 


**** BATDORF MODEL CRACK ORIENTATION, CRACK SHAPE, **»* 

AND FRACTURE CRITERION ARE CONSIDERED 


— TEMPERATURE DEPENDENT MATERIAL PARAMETERS FOR MATERIAL NUMBER 

WEIBULL MODULUS (SHAPE PARAMETER) , M (DIMENSIONLESS) 

NORMALIZED BATDORF CRACK DENSITY COEFFICIENT, K (DIMENSIONLESS) 
SCALE PARAMETER, SP (UNITS OF STRESS*VOLUME** (1/M)) 


TEMPERATURE 

M BIASED 

K 

SP 

70.0000 

. 4130E+01 

. 5130E+01 

. 1787E+02 


FRACTURE CRITERION = COPLANAR STRAIN ENERGY RELEASE RATE CRITERION 
CRACK SHAPE = GRIFFITH CRACK 


* * 

* ECHO OF MATERIAL CONTROL INPUT * 

* * 


TITLE = SI3N4 SPECIMEN DATA FROM ASEA CERAMA FOR SURFACE FLAW ANALYSIS 


2 = MATERIAL IDENTIFICATION NUMBER (MATID) 


2 = CONTROL INDEX FOR VOLUME OR SURFACE FLAW ANALYSIS (ID4) 

1 : VOLUME 

2 : SURFACE 
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0 = CONTROL INDEX FOR METHOD OF CALCULATING WEIBULL PARAMETERS 
FROM THE EXPERIMENTAL FRACTURE DATA (MLORLE) 

0 : MAXIMUM LIKELIHOOD 

1 : LEAST-SQUARES LINEAR REGRESSION 


5 = CONTROL INDEX FOR EXPERIMENTAL DATA (ID1) 

1 : UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

2 : FOUR-POINT BEND TEST DATA 

3 : DIRECT INPUT OF THE WEIBULL PARAMETERS, M AND SP 

(shape parameter and scale parameter) 

4 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

UNIFORM UNIAXIAL TENSILE SPECIMEN TEST DATA 

5 : CENSORED DATA FOR SUSPENDED ITEM ANALYSIS OF 

FOUR-POINT BEND TEST DATA 


3 = CONTROL INDEX FOR SURFACE FRACTURE CRITERION (ID2S) 
1 : NORMAL STRESS FRACTURE CRITERION 
(SHEAR-INSENSITIVE CRACK) 

3 : COPLANAR STRAIN ENERGY RELEASE RATE CRITERION 

4 : WEIBULL PIA MODEL 

5 : SHETTY’S SEMI-EMPIRICAL CRITERION 


1 = CONTROL INDEX FOR SHAPE OF SURFACE CRACKS (ID3S) 

1 : GRIFFITH CRACK 

(ASSOCIATED WITH STRAIN ENERGY RELEASE RATE CRITERION) 
(ASSOCIATED WITH SHETTY’S SEMI -EMPIRICAL CRITERION) 

3 : GRIFFITH NOTCH 

(ASSOCIATED WITH STRAIN ENERGY RELEASE RATE CRITERION) 
(ASSOCIATED WITH SHETTY’S SEMI-EMPIRICAL CRITERION) 

4 : SEMICIRCULAR CRACK 

(ASSOCIATED WITH SHETTY’S SEMI-EMPIRICAL CRITERION) 


1 = CONTROL INDEX FOR METHOD OF CALCULATING BATDORF CRACK DENSITY 
COEFFICIENT 

(K SUB B) FROM TEST DATA (IKBAT) 

0 : SHEAR-INSENSITIVE METHOD (MODE I FRACTURE ASSUMED) 

1 : SHEAR-SENSITIVE METHOD (FRACTURE ASSUMED TO OCCUR 

ACCORDING TO THE FRACTURE CRITERION AND CRACK 
SHAPE SELECTED BY THE ID2 AND ID3 INDICES) 


.2500 = POISSON’S RATIO (PR) 


* * 

* STATISTICAL ANALYSIS OF FRACTURE SPECIMEN DATA * 

* * 
**************************************************** 
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ECHO OF SPECIMEN INPUT DATA, IN ASCENDING ORDER OF FRACTURE STRESS 


. 4000E-01 = OUTER LOAD SPAN OF FOUR -POINT BEND BAR 
. 2000E-01 = INNER LOAD SPAN 
. 3500E-02 = DEPTH OF SPECIMEN 
. 4500E-02 = WIDTH OF SPECIMEN 


79 = NUMBER OF SPECIMENS IN BATCH 70.000 = TEMPERATURE OF BATCH 

"S"URFACE OR "VOLUME OR "U"NKNOWN FLAW ORIGIN AND RESPECTIVE FAILURE STRESS 


vvv 

.4159E+03 

. 4584E+03 

vw 

.5272E+03 

. 5458E+03 

vw 

.5725E+03 

. 5950E+03 

sss 

. 6045E+03 

. 6093E+03 

sss 

. 6144E+03 

. 6210E+03 

sss 

. 6216E+03 

. 6222E+03 

sss 

. 6251E+03 

. 6258E+03 

sss 

. 6402E+03 

. 6426E+03 

sss 

.6501E+03 

. 6516E+03 

sss 

. 6569E+03 

. 6570E+03 

sss 

. 6597E+03 

. 6617E+03 

sss 

-6621E+03 

. 6636E+03 

sss 

. 6645E+03 

. 6664E+03 

sss 

. 6708E+03 

. 6712E+03 

sss 

.6721E+03 

. 6741E+03 

sss 

. 6771E+03 

. 6784E+03 

sss 

. 6831E+03 

. 6837E+03 

sss 

. 6867E+03 

. 6870E+03 

sss 

.6940E+03 

. 6950E+03 

sss 

. 7029E+03 

. 7029E+03 

sss 

. 7030E+03 

. 7045E+03 

sss 

•7061E+03 

. 7101E+03 

sss 

.7156E+03 

-7156E+03 

sss 

.7164E+03 

. 7165E+03 

sss 

.7252E+03 

. 7254E+03 

sss 

. 7260E+03 

. 7271E+03 

s 

.7323E+03 



— STEFANSKY OUTLIER TEST 


. 5199E+03 
. 5606E+03 
. 6035E+03 
. 6115E+03 
. 6215E+03 
. 6224E+03 
. 6311E+03 
. 6490E+03 
. 6551E+03 
. 6597E+03 
. 6621E+03 
. 6641E+03 
. 6688E+03 
. 6717E+03 
. 6768E+03 
. 6799E+03 
. 6855E+03 
. 6911E+03 
. 6998E+03 
.7029E+03 
.7045E+03 
. 7126E+03 
. 7158E+03 
. 7173E+03 
. 7254E+03 
. 7287E+03 


SPECIMEN FRACTURE STRESSES — 


RESULTS FROM THE STEFANSKY OUTLIER TEST FOR TEMP. = 70.0000 

FAILURE STRESS 

.4159E+03 DEVIATES FROM THE MAIN TREND OF THE DATA AT THE 1 % 
SIGNIFICANCE LEVEL 

. 4584E+03 DEVIATES FROM THE MAIN TREND OF THE DATA AT THE 1 % 
SIGNIFICANCE LEVEL 

.5199E+03 DEVIATES FROM THE MAIN TREND OF THE DATA AT THE 1 % 
SIGNIFICANCE LEVEL 

. 5272E+03 DEVIATES FROM THE MAIN TREND OF THE DATA AT THE 5% 
SIGNIFICANCE LEVEL 

.5458E+03 DEVIATES FROM THE MAIN TREND OF THE DATA AT THE 10% 
SIGNIFICANCE, LEVEL 
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DEVIATION FROM THE MAIN TREND OF THE DATA MAY 
INDICATE BAD VALUES. MULTIPLE DEVIATIONS FROM THE SAME 
REGION OF THE DISTRIBUTION INDICATE THAT EITHER A CONCURRENT 
OR A PARTIALLY CONCURRENT FLAW POPULATION HAS BEEN DETECTED 
(NOTE THAT A CONCURRENT FLAW POPULATION MAY BE PRESENT BUT NOT 
BE DETECTED BY THE OUTLIER TEST) . DEVIATIONS OCCURRING 
IN THE SAME REGION OF THE DISTRIBUTION WITH ALL THREE 
SIGNIFICANCE LEVELS ( 1 %, 5% AND 10*) PRESENT INDICATE A 
CONCURRENT FLAW POPULATION. DEVIATIONS SHOULD BE EXAMINED AND 
TREATED ACCORDINGLY (I.E. IGNORE, CENSOR, ADJUST OR ELIMINATE 
STRESS) . JUDGEMENT OF ACTION TAKEN CAN BE DETERMINED FROM THE 
GOODNESS-OF-FIT TESTS. 

THE OUTLIER TEST IS NO SUBSTITUTE FOR GRAPHICAL EXAMINATION!!! 


— TEMP. DEP. WEIBULL MODULUS AND CHARACTERISTIC STRENGTH WITH 90% — 
CONFIDENCE BOUNDS DETERMINED BY MAXIMUM LIKELIHOOD ANALYSIS 


NOTE: 90% CONFIDENCE BOUNDS ON PARAMETERS DETERMINED FROM COMPETING FAILURE 
MODES (CENSORED DATA) ARE APPROXIMATE . 

FOR CENSORED DATA THE UNBIASED VALUE OF THE PARAMETER "M" IS NOT 
GIVEN. 


FOR 

SAMPLE SIZES 

LESS THAN 

4, CONFIDENCE 

LIMITS ARE 

NOT GIVEN. 


TEMP. 

M BIASED 

M UNBIASED 

UP M 

LOW M 


CHAR. STR. 

UP C.S. 

LOW C.S. 

MEAN 

STD. DEV. 

VOLUME 

70 . 0000 

. 4130E+01 


. 5666E+01 

. 2179E+01 


. 1128E+04 

. 1342E+04 

. 9542E+03 

. 1024E+04 

. 2791E+03 

SURFACE 

70.0000 

. 2281E+02 


. 2619E+02 

. 1920E+02 


. 6917E+03 

. 6981E+03 

. 6853E+03 

. 6755E+03 

. 3684E+02 


STATISTICS FROM THE GOODNESS-OF-FIT TESTS FOR TEMP. = 70.0000 


KOLMOGOROV-SMIRNOV TEST 


ORDER 

FRAC. STR. 

WEIB. PROB. OF FAIL. 

D+ FACTOR 

D- FACTOR 

SIGNIF. LEVEL 

1 

.4159E+03 

.0161 

-.0035 

.0161 

99.0000 

2 

.4584E+03 

.0240 

.0013 

.0114 

99 . 0000 

3 

.5199E+03 

.0414 

-.0034 

.0161 

99.0000 

4 

. 5272E+03 

.0443 

.0064 

.0063 

99.0000 

5 

. 5458E+03 

.0529 

.0104 

.0023 

99.0000 

6 

.5606E+03 

.0620 

.0140 

-.0013 

99.0000 

7 

.5725E+03 

.0714 

.0172 

-.0045 

99.0000 

8 

. 5950E+03 

.0983 

.0029 

.0097 

99.0000 
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9 

. 6035E+03 

.1130 

.0009 

.0118 

99.0000 

10 

. 6045E+03 

.1151 

.0114 

.0012 

99.0000 

11 

. 6093E+03 

.1254 

.0139 

-.0012 

99.0000 

12 

. 6115E+03 

.1305 

.0214 

-.0088 

99.0000 

13 

. 6144E+03 

.1379 

.0267 

-.0140 

99.0000 

14 

. 6210E+03 

.1567 

.0205 

-.0079 

99.0000 

15 

. 6215E+03 

.1581 

.0317 

-.0191 

99.0000 

16 

. 6216E+03 

.1585 

.0440 

-.0313 

99.0000 

17 

. 6222E+03 

.1603 

.0549 

-.0422 

97.1339 

18 

. 6224E+03 

.1610 

.0668 

-.0541 

87.2422 

19 

. 6251E+03 

.1704 

.0701 

-.0575 

83.1735 

20 

. 6258E+03 

.1726 

.0805 

-.0679 

68.5012 

21 

. 6311E+03 

.1930 

.0728 

-.0602 

79.6191 

22 

. 6402E+03 

.2346 

.0439 

-.0312 

99.0000 

23 

. 6426E+03 

.2476 

.0435 

-.0309 

99.0000 

24 

. 6490E+03 

.2853 

.0185 

- . 0058 

99 . 0000 

25 

. 6501E+03 

.2924 

.0241 

-.0114 

99.0000 

26 

. 6516E+03 

.3020 

.0271 

-.0145 

99.0000 

27 

. 6551E+03 

.3262 

.0156 

-.0030 

99.0000 

28 

. 6569E+03 

.3399 

.0146 

-.0019 

99 . 0000 

29 

. 6570E+03 

.3407 

.0264 

-.0137 

99 . 0000 

30 

. 6597E+03 

.3616 

.0181 

-.0055 

99 . 0000 

31 

. 6597E+03 

.3616 

.0308 

-.0181 

99.0000 

32 

. 6617E+03 

.3773 

.0277 

-.0151 

99.0000 

33 

. 6621E+03 

.3804 

.0374 

-.0247 

99.0000 

34 

. 6621E+03 

.3811 

.0493 

-.0366 

99.0000 

35 

. 6636E+03 

.3935 

.0496 

-.0369 

99.0000 

36 

. 6641E+03 

.3980 

.0577 

-.0450 

95.5201 

37 

. 6645E+03 

.4015 

.0669 

-.0542 

87.1359 

38 

. 6664E+03 

.4182 

.0628 

-.0502 

91.3909 

39 

. 6688E+03 

.4398 

.0538 

-.0412 

97.5998 

40 

. 6708E+03 

.4580 

.0484 

-.0357 

99.0000 

41 

. 6712E+03 

.4619 

.0571 

-.0445 

95.8855 

42 

. 6717E+03 

.4674 

.0643 

-.0516 

89.9719 

43 

. 6721E+03 

.4709 

.0734 

-.0607 

78.8750 

44 

. 6741E+03 

.4909 

.0661 

-.0534 

88.0606 

45 

. 6768E+03 

.5182 

.0515 

-.0388 

98.4919 

46 

. 6771E+03 

.5207 

.0616 

-.0489 

92.5627 

47 

. 6784E+03 

.5338 

.0611 

- . 0484 

92.9632 

48 

. 6799E+03 

.5495 

.0581 

-.0454 

95.2717 

49 

. 6831E+03 

.5842 

.0360 

-.0234 

99.0000 

50 

. 6837E+03 

.5909 

.0420 

-.0293 

99.0000 

51 

. 6855E+03 

.6100 

.0356 

-.0229 

99.0000 

52 

. 6867E+03 

.6237 

.0345 

-.0219 

99 . 0000 

53 

- 6870E+03 

.6264 

.0445 

-.0318 

99.0000 

54 

. 6911E+03 

.6712 

.0124 

.0003 

99.0000 

55 

. 6940E+03 

.7027 

-.0065 

.0192 

99.0000 

56 

. 6950E+03 

.7135 

-.0046 

.0173 

99.0000 

57 

. 6998E+03 

.7635 

-.0419 

.0546 

97.2613 

58 

.7029E+03 

.7946 

-.0604 

.0731 

79.2561 

59 

.7029E+03 

.7946 

-.0478 

.0604 

93.5064 

60 

.7029E+03 

.7950 

-.0355 

.0482 

99.0000 

61 

. 7030E+03 

.7955 

-.0233 

.0360 

99.0000 

62 

. 7045E+03 

.8099 

-.0251 

.0378 

99.0000 

63 

. 7045E+03 

.8105 

-.0130 

.0257 

99.0000 

64 

.7061E+03 

.8247 

-.0146 

.0272 

99.0000 

65 

.7101E+03 

.8598 

-.0370 

.0496 

98.9973 

66 

.7126E+03 

.8801 

-.0447 

.0574 

95.7363 

67 

.7156E+03 

.9022 

-.0541 

.0667 

87 . 3420 

68 

. 7156E+03 

.9022 

-.0414 

.0541 

97.5075 
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69 

. 7158E+03 

.9031 

70 

. 7164E+03 

.9076 

71 

.7165E+03 

.9081 

72 

.7173E+03 

.9134 

73 

. 7252E+03 

.9549 

74 

. 7254E+03 

.9560 

75 

. 7254E+03 

.9560 

76 

.7260E+03 

.9581 

77 

. 7271E+03 

.9624 

78 

. 7287E+03 

.9680 

79 

. 7323E+03 

.9786 


- . 0297 

.0424 

99.0000 

- . 0215 

.0341 

99.0000 

-.0094 

.0220 

99.0000 

-.0020 

.0146 

99.0000 

-.0309 

.0435 

99.0000 

-.0193 

.0320 

99.0000 

- . 0067 

.0193 

99.0000 

.0039 

.0087 

99 . 0000 

.0123 

.0003 

99.0000 

.0193 

- . 0067 

99.0000 

.0214 

-.0088 

99.0000 


KOLMOGOROV-SMIRNOV TEST YIELDS STATISTIC D = MAX (D+.D-) = 0805 

WITH AN ASSOCIATED SIGNIFICANCE LEVEL OF 68.5% 


ANDERSON-DARLING TEST YIELDS STATISTIC A* *2 = .6725 

WITH AN ASSOCIATED SIGNIFICANCE LEVEL OF 58.3% 


KANOFSKY-SRINTVASAN 90% CONFIDENCE BANDS ABOUT THE WEIBULL DISTRIBUTION FOR 

TEMP. = 70.0000 


THE KANOFSKY-SRINTVASAN FACTOR FOR THIS DISTRIBUTION IS .0836 FOR A SAMPLE 

SIZE OF 79 


ORD. FRAC. STR. WEIB. PROB. OF FAIL. UPP . CONF . BAND MED. RANK LOW . CONF . BAND 


1 

.4159E+03 

.0161 

2 

.4584E+03 

.0240 

3 

. 5199E+03 

.0414 

4 

. 5272E+03 

.0443 

5 

, 5458E+03 

.0529 

6 

.5606E+03 

.0620 

7 

.5725E+03 

.0714 

8 

. 5950E+03 

.0983 

9 

. 6035E+03 

.1130 

10 

. 6045E+03 

.1151 

11 

. 6093E+03 

.1254 

12 

. 6115E+03 

.1305 

13 

. 6144E+03 

.1379 

14 

.6210E+03 

.1567 

15 

. 6215E+03 

.1581 

16 

. 6216E+03 

.1585 

17 

. 6222E+03 

.1603 

18 

. 6224E+03 

.1610 

19 

. 6251E+03 

.1704 

20 

. 6258E+03 

.1726 

21 

. 6311E+03 

.1930 

22 

. 6402E+03 

.2346 

23 

. 6426E+03 

.2476 

24 

. 6490E+03 

.2853 

25 

. 6501E+03 

.2924 

26 

. 6516E+03 

.3020 

27 

. 6551E+03 

.3262 

28 

. 6569E+03 

.3399 

29 

.6570E+03 

.3407 


.0997 

.0088 

.0000 

.1077 

.0214 

.0000 

.1250 

.0340 

.0000 

.1279 

.0466 

.0000 

.1365 

.0592 

.0000 

.1456 

.0718 

.0000 

.1551 

.0844 

.0000 

.1820 

.0970 

.0147 

.1967 

.1096 

.0294 

.1988 

.1222 

.0315 

.2090 

.1348 

.0418 

.2141 

.1474 

.0469 

.2215 

.1599 

.0542 

.2403 

.1725 

.0731 

.2418 

.1851 

.0745 

.2422 

.1977 

.0749 

.2440 

.2103 

.0767 

.2447 

.2229 

.0774 

,2540 

.2355 

.0867 

,2563 

.2481 

.0890 

,2766 

.2607 

.1094 

3182 

.2733 

.1510 

3313 

.2859 

.1640 

3690 

.2985 

.2017 

3760 

.3111 

.2088 

3856 

.3237 

.2183 

4098 

.3363 

.2425 

4235 

.3489 

.2562 

4243 

.3615 

.2571 
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30 

. 6597E+03 

.3616 

31 

. 6597E+03 

.3616 

32 

. 6617E+03 

.3773 

33 

. 6621E+03 

.3804 

34 

. 6621E+03 

.3811 

35 

. 6636E+03 

.3935 

36 

. 6641E+03 

.3980 

37 

. 6645E+03 

.4015 

38 

. 6664E+03 

.4182 

39 

. 6688E+03 

.4398 

40 

. 6708E+03 

.4580 

41 

. 6712E+03 

.4619 

42 

. 6717E+03 

.4674 

43 

. 6721E+03 

.4709 

44 

. 6741E+03 

.4909 

45 

.6768E+03 

.5182 

46 

. 6771E+03 

.5207 

47 

. 6784E+03 

.5338 

48 

. 6799E+03 

.5495 

49 

. 6831E+03 

.5842 

50 

. 6837E+03 

.5909 

51 

. 6855E+03 

.6100 

52 

. 6867E+03 

.6237 

53 

. 6870E+03 

.6264 

54 

. 6911E+03 

.6712 

55 

. 6940E+03 

.7027 

56 

. 6950E+03 

.7135 

57 

. 6998E+03 

.7635 

58 

.7029E+03 

.7946 

59 

. 7029E+03 

.7946 

60 

. 7029E+03 

.7950 

61 

. 7030E+03 

.7955 

62 

.7045E+03 

.8099 

63 

.7045E+03 

.8105 

64 

, 7061E+03 

.8247 

65 

.7101E+03 

.8598 

66 

. 7126E+03 

.8801 

67 

. 7156E+03 

.9022 

68 

. 7156E+03 

.9022 

69 

. 7158E+03 

.9031 

70 

, 7164E+03 

.9076 

71 

.7165E+03 

.9081 

72 

.7173E+03 

.9134 

73 

.7252E+03 

.9549 

74 

.7254E+03 

.9560 

75 

. 7254E+03 

.9560 

76 

. 7260E+03 

.9581 

77 

. 7271E+03 

.9624 

78 

. 7287E+03 

.9680 

79 

. 7323E+03 

.9786 


.4453 

.3741 

.2780 

.4453 

.3866 

.2780 

.4609 

.3992 

.2937 

.4640 

.4118 

.2967 

.4647 

.4244 

.2975 

.4771 

.4370 

.3098 

.4816 

.4496 

.3144 

.4851 

.4622 

.3178 

.5018 

.4748 

.3345 

.5234 

.4874 

.3562 

.5416 

.5000 

.3743 

.5455 

.5126 

. .3782 

.5510 

.5252 

.3837 

. 5546 

.5378 

.3873 

.5745 

.5504 

.4073 

.6018 

.5630 

.4345 

.6043 

.5756 

.4371 

.6175 

.5882 

.4502 

.6332 

.6008 

.4659 

.6679 

.6134 

.5006 

.6745 

.6259 

.5073 

.6936 

.6385 

.5263 

.7073 

.6511 

.5401 

.7100 

.6637 

.5428 

.7548 

.6763 

.5875 

.7864 

.6889 

.6191 

.7971 

.7015 

.6298 

.8471 

.7141 

.6798 

.8782 

.7267 

.7110 

.8782 

.7393 

.7110 

.8786 

.7519 

.7114 

.8791 

.7645 

.7119 

.8936 

.7771 

.7263 

.8941 

.7897 

.7269 

.9083 

.8023 

.7411 

.9434 

.8149 

.7761 

.9638 

.8275 

.7965 

.9858 

.8401 

.8185 

.9858 

.8526 

.8185 

.9867 

.8652 

.8195 

.9912 

.8778 

.8239 

.9917 

.8904 

.8245 

.9970 

.9030 

.8297 

1.0000 

.9156 

.8713 

1 . 0000 

.9282 

.8724 

1 . 0000 

.9408 

.8724 

1 . 0000 

.9534 

.8745 

1 . 0000 

.9660 

.8787 

1 . 0000 

.9786 

.8844 

1 . 0000 

.9912 

.8950 
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* SURFACE FLAW PARAMETER ANALYSIS * 

* * 


**** BATDORF MODEL — CRACK ORIENTATION, CRACK SHAPE, **** 
AND FRACTURE CRITERION ARE CONSIDERED 


TEMPERATURE DEPENDENT MATERIAL PARAMETERS FOR MATERIAL NUMBER 2 

WEIBULL MODULUS (SHAPE PARAMETER) , M (DIMENSIONLESS) 

NORMALIZED BATDORF CRACK DENSITY COEFFICIENT, K (DIMENSIONLESS) 
SCALE PARAMETER, SP (UNITS OF STRESS*AREA**(1/M)) 

TEMPERATURE M BIASED K SP 

70.0000 . 2281E+02 .6050E+01 .4613E+03 


FRACTURE CRITERION = COPLANAR STRAIN ENERGY RELEASE RATE CRITERION 
CRACK SHAPE = GRIFFITH CRACK 
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7.0 CODE DESCRIPTION 


Included on the distribution disk is a file called PCCARES.FOR which con- 
tains the FORTRAN source code for the PC-CARES program. The code was written 
with VAX and Microsoft extensions and was compiled and linked with Microsoft 
FORTRAN 5.0. The specific extensions used include the D0...ENDD0, the DO 
WHILE ... ENDDO, the ALLOCATE and DEALLOCATE routines and the EOF function. 

Those wishing to modify the code and compile it with a compiler other than 
Microsoft FORTRAN 5.0 may have to alter the sections of the code which utilize 
these extensions. For the rest of this description, you may want to have a 
copy of the code in front of you. 

The ALLOCATE and DEALLOCATE routines are called from the following 
PC-CARES routines: MAIN, MATL, CRACKV, CRACKS, and NORMAL. These calls should 

be removed. However note how they are dimensioned in the ALLOCATE call so that 
when you dimension them normally you will know to what size. You will have to 
remove the ALLOCATABLE attribute from the DIMENSION statements as well, when 
you insert the dimension numbers. For those arrays dimensioned using IMAXF and 
IMAXT parameters, it is suggested that you simply dimension the arrays using 
these labels and then set them to interger constants in the MAIN routine using 
the FORTRAN PARAMETER statement. The default PC-CARES values are IMAXF = 150 
and IMAXT = 10. In addition the READINI subroutine should be altered to remove 
the assignment of the IMAXF and IMAXT parameters. 

The EOF function is used in only two routines: the MAIN routine and the 

READINI subroutine. The code can be modified to support whatever end of file 
support your particular compiler supplies or the calls may simply be removed. 
Note that the $END keywords normally stop the FORTRAN file reading, however, if 
those keywords are missing then your code may attempt to read beyond the end of 
file and abort abnormally generating a run-time error. 

The DO... ENDDO and DO WHILE. .. ENDDO extensions are used throughout the 
code, so if your compiler does not support these extensions then you will have 
to change these structures to standard FORTRAN 77 by using the DO label ... label 
CONTINUE structure for the DO... END DO and by substituting an appropriately 
placed IF and GOTO statement combination for the DO WHILE. .. ENDDO. An example 
1 s given as fol lows : 

Microsoft/VAX Extension Standard FORTRAN 

DO I = 1, 1000 DO 100 I = 1, 1000 


END DO 


100 CONTINUE 
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DO WHILE (CONDITION) 


105 IF (.NOT. (CONDITION)) GOTO 110 


END DO GOTO 105 

110 CONTINUE 

Finally here Is a list of descriptions of the routines used In PC-CARES. 
The MAIN Routine 

The PC-CARES main routine controls the logical flow of the PC-CARES pro- 
gram as diagrammed In figure 2.1. Specifically the MAIN routine Initializes 
and dimensions the code's variables and arrays, after which It calls the 
READINI subroutine to read the Initialization file. If the initialization file 
Is not present then the Initialization parameters are set to their default 
values, the routine then allocates the array space needed by the program by 
calling the Microsoft FORTRAN extension ALLOCATE, after which the PC-CARES 
input file Is read to obtain the Master Control Input and the Material Control 
Input, echoing the Input along the way by calling the subroutines PRINTA and 
PRINTB. The routine continues by reading the fracture stress data and their 
respective fracture origins If present. If the material parameters (specifi- 
cally the shape and scale parameters) are supplied, they are read from the 
Input file. Finally the subroutine MATL Is called which performs all the sta- 
tistical analysis of the material. Following the MATL subroutine call, the 
MAIN routine deallocates all of the allocated array space using the DEALLOCATE 
command from Microsoft FORTRAN, closes all the files and exits. 

In addition to the main program, the following subroutines appear in the 
PC-CARES listings: 

Subroutine ANGLE 

This subroutine evaluates Q(E,o cr ) or u(r,a cr ) for volume and/or sur- 
face flaw analysis, respectively, when the Batdorf method Is selected. ANGLE 
employs the quadratic solution procedure described In the appendix A.l Volume 
Flaw Reliability and the appendix A. 2 Surface Flaw Reliability sections of 
this manual. It determines the critical Intervals where o e > o cr for various 
angles of a and values of a cr about the unit sphere for volume flaw analy- 
sis. For some specific stress states (02 = 03 ), this evaluation Is Independent 
of 0. For surface flaw analysis these Intervals are determined about the unit 
circle. For volume flaw analysis the critical Intervals correspond to the 3 
Integral In equation (A. 27). For surface flaw analysis the critical Intervals 
correspond to the Integral In equation (A. 55). Figure 2.2 shows the fracture 
criteria and flaw geometries for which the coding has been developed. These 
correspond to the effective stress equations (A. 18), (A. 19), (A. 22), (A. 23), 

(A. 25), (A. 26), and (A. 40) for volume flaw analysis and to equations (A. 51) to 
(A. 54) for surface flaw analysis. The equations listed In tables A.l to A. Ill 
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are used to find when o e - o cr , and the procedure outlined in equations (A. 36) 
and (A. 37) is used to find the intervals where c e > o cr . 

ANGLE is called from the MATBAT subroutine In order to aid in the calcu- 
lation of the Batdorf crack density coefficient when IKBAT = 1 (in other words 
when the user specifies the fracture criterion and the crack geometry and wants 
shear sensitivity taken into account). Arguments R2 to R4 correspond to 02 , 
< 73 , and (02 - < 73 ), divided by cti for the given stress state. Arguments P and 
Q represent squared trigonometric functions of the angles a or ( 3 , whereas 
argument H represents the values of a cr at locations of the Gaussian quadra- 
ture points. These arguments are required within ANGLE to calculate Intermedi- 
ate variables ai to a 3 that are coefficients in the quadratic equation for 
cos 2 a or cos 2 |3 listed in tables A. I to A. III. 

Subroutine ANGLES 

This subroutine is used with the Batdorf volume flaw model to integrate 
over the surface area of a quadrant of the unit sphere when the Shetty failure 
criterion Is used and <j p „, u > 01 . ANGLES determines the intervals where 

a e 2 a cr f° r constant angles of a about the unit sphere and stores the 
limits of these Intervals in the INTVAL array. The critical Intervals corre- 
spond to the Integral of |3 described In equation (A. 27). The limits of these 
Intervals are determined for each transformed Gaussian value of o cr and a. 
Each consecutive pair of integers in the third index of array INTVAL represents 
an Interval where o e > a cr for an angle of a denoted by the second Index. 
The first index corresponds to values of a cr at locations of the Gaussian 
quadrature points. The limits of Integration stored In INTVAL are integers 
representing 1° Increments of angle (3 counted from -ir/2 to ir/2. 

Function CONLIC 

This function performs a table lookup of the factors for obtaining 90 per- 
cent upper and lower confidence bounds of the MLE of oq . These factors have 
been taken from reference 10. They are obtained from a Monte-Carlo simulation 
by using maximum likelihood analysis and uncensored data. The confidence bound 
calculations are performed in subroutine MATL. 

Function CONLIM 

This function performs a table lookup of the factors for obtaining 90 per- 
cent upper and lower confidence bounds of the MLE of m. These factors have 
been taken from reference 10. They are obtained from a Monte-Carlo simulation 
by using maximum likelihood analysis and uncensored data. The confidence bound 
calculations are performed in subroutine MATL. 

Subroutine CRACKS 

Subroutine CRACKS is called from MATBAT and serves as an interface with 
the SORMAL, SNGLES, SVALP3, and FINDP subroutines for the calculation of the 
shear-sensitive (IKBAT = 1) normalized Batdorf surface crack density coeffi- 
cient kgs for the Shetty criterion when °e max > ° 1 • 
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Subroutine CRACKV 


Subroutine CRACKV Is called from MATBAT and serves as an Interface with 
the NORMAL, ANGLES AND EVALP3 subroutines for the calculation of the shear-_ 
sensitive (IKBAT = 1) normalized Batdorf volume crack density coefficient kgv 
for the Shetty criterion when o emax > • 

Subroutine EVALP3 

Subroutine EVALP3 Is used with the Batdorf model for volume flaw analysis 
with the Shetty failure criterion when <7e max > • It performs the Integra- 

tion 



which Is used In equations (A. 27) and (A. 38a). Legendre-Gauss quadrature is 
used for the numerical Integrations of da and da cr . The stored values In 
the INTVAL array previously calculated in the ANGLES subroutine are used to 
perform the integration. 

Function F 

This function computes the polynomial approximation to the Gamma function 
as per the "Handbook of Mathematical Functions." 

Subroutine FINDP 

Subroutine FINDP Is used with the Batdorf model for surface flaw analysis 
with the Shetty failure criterion when o 8max > a] . It calculates P 2 S as 
defined by equation (A. 55). The interval Is determined from transforming val- 
ues stored in the INTVAL array Into real numbers. FINDP is called from the 
SVALP3 subroutine. 

Subroutine GAUSS 

Subroutine GAUSS contains roots of the Legendre polynomials and the weight 
factors for the Gauss quadrature. It Is employed In the calculation of the 
Batdorf crack density coefficient, when a closed-form solution Is not avail- 
able. The number of Gauss points (NGP) Is specified by the user in the program 
Input. Data are available In GAUSS for NGP ■ 2 to 10, 15, and 30 although only 
15 and 30 are recommended. The weights and locations are contained In the W 
and H arrays, respectively. 

Subroutine LEAST2 

This subroutine calculates the Welbull strength parameters m and C 
(see eq. (A. 59)) by using the least-squares analysis method for complete or 
censored samples. The slope m and the Intercept (In (O) of the line of 
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best fit are obtained by solving two simultaneous equations (ref. 17). For 
uncensored data, median rank regression analysis (eq. (A. 61)) is used to 
calculate the failure probability, Pf. However, in case of censored data, the 
median rank regression analysis cannot be used directly because of the effect 
of competing failure modes. Instead, the rank increment technique (eq. 

(A. 62)) is used to adjust rank values. These adjusted rank values are then 
used with median rank regression analysis to calculate the failure 
probability. The LEAST2 subroutine is called from the MATL subroutine. 

Subroutine MATBAT 

Subroutine MATBAT calculates the surface and/or volume scale parameters, 
oqS and a 0 y, and the normalized Batdorf crack density coefficients, kgs and 
Rbv* respectively. For a given material, parameters are found for each temper- 
ature level that is input by the user. If m$ and o 0 s or my and o 0 y are 
directly input, then only Rg$ or Rgy is calculated, respectively. If 
experimental fracture stresses are input for either four-point bend or uniaxial 
tensile specimens, then all required parameters are calculated. The scale 
parameter for volume flaws is calculated from equation (A. 72), and for surface 
flaws, equation (A. 87) is used. The scale parameter is determined from the 
specimen geometry and from the values of m and C estimated in the LEAST2 
or MAXL subroutines. The coefficient Rg$ calculated from equation (A. 91) 
or (A. 94), and kgy is calculated from equation (A. 78) or (A. 84). The ANGLE 
subroutine Is called from MATBAT to evaluate co(E,o cr )/2ir or Q(E,o cr )/4ir for 
a uniaxial stress state to find Rg$ or Rgy. If the Shetty criterion is 
selected by the user with IKBAT = 1 option when a emax > oi , then the CRACKS 

subroutine finds w(E,a cr )/2ir and the CRACKV subroutine finds Q(E,o cr )/4ir to 
calculate Rg$ and kgy, respectively. 

Subroutine MATL 

Subroutine MATL controls the program logic flow for the determination of 
the statistical material parameters and other useful statistical quantities as 
shown In the flowchart of figure 2.1. Kolmogorov-Smlrnov and Anderson-Darling 
goodness-of-f 1 t tests, Kanofsky-Srinivasan 90-percent confidence bands, Weibull 
mean, Weibull variance, and 90-percent confidence bounds on the parameters are 
all calculated in this subroutine. Ancillary subroutines to detect outliers 
(OUTLIE), to perform least-squares (LEAST2) or maximum likelihood (MAXL) analy- 
sis, to calculate Weibull scale parameters and the Batdorf crack density coef- 
ficient (MATBAT), and to print out results of the analysis (PRINTP) are called 
from MATL. 

Subroutine MAXL 

This subroutine determines the MLE's of the Weibull strength Parameters m 
and o@ by using the maximum likelihood method for both uncensored and cen- 
sored data. The logarithm of the likelihood function Is differentiated with 
respect to m and C, and the resulting expressions are set equal to zero 
(eqs. (A. 59), and (A. 63) to (A. 65)). The Newton-Raphson Iterative technique is 
used to obtain the parameter MLE's by solving these nonlinear equations. The 
estimate for m for the first iteration is obtained from least-squares analy- 
sis via the LEAST2 subroutine. If the convergence criterion it not met after 
50 Iterations, the maximum likelihood method is terminated, a warning message 
Is printed, and the results from the least-squares analysis are subsequently 
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used in the program. If the fracture data are a complete sample, then the 
unbiased estimate of the shape parameter is calculated with the factors stored 
in the UNBIAS subroutine. The unbiased estimate is passed to the MATL subrou- 
tine and is later printed out in the PRINTP subroutine. This parameter is not 
employed in any subsequent calculations. Reference 17 contains a detailed 
description of the method of calculation of these statistical quantities. The 
MAXL subroutine is called from the MATL subroutine. 

Subroutine NORMAL 

This subroutine is called when the Batdorf volume flaw model is used with 
the Shetty failure criterion and a e max > a l • NORMAL calculates the normalized 

effective stress about the unit sphere as a function of the angles a and 3, 
and stores these values in the SEANGL array. The effective stress is deter- 
mined for P increments of 3 and i s stored in the second index of the array. 
The first index denotes angles of a with values corresponding to the trans- 
formed Gaussian points. The array values are normalized by the maximum effec- 
tive stress o e found for the stress state being evaluated. The effective 

stress is calculated from equations (A. 25) and (A. 26). 

Subroutine OUTLIE 

In this subroutine, the available specimen fracture stress data at each 
temperature level are examined for outliers or inconsistent data. At the start 
of the subroutine, the sample mean and sample standard deviation are calcula- 
ted. From these values, the normed residual for each specimen is obtained 
(ref. 6). The normed residuals are normalized deviations of the data about the 
sample mean. The Weibull distribution is not symmetrical about its mean, and 
therefore this technique is only approximate. The absolute maximum of the 
normed residual (MNR) statistic is compared with the critical value (CV) at 1-, 
5-, and 10-percent significance levels. If the MNR statistic is smaller than 
the three critical values, then no outliers are detected. However, if the MNR 
is larger than at least one of the three critical values, the corresponding 
data value with the MNR statistic is detected as an outlier with the appro-, 
priate si gnlf ic&nce level. The outlier test can only flag one point per trial 
as an outlier. If an outlier is detected at the 10-percent or less signif- 
icance level, it is removed from the sample and the remaining points are then 
retested. This process is repeated until the sample is reduced such that no 
more outliers are detected. Once all such points are detected, each of these 
potential outliers is retested against the remaining "good" data, and those 
points that maintain significance levels at or below 10 percent are flagged 
with the appropriate significance level in the ISKIP array. The results of 
the outlier test are output in the PRINTP subroutine via the MATL subroutine. 

A discussion of the equations used with this method is given in reference 17. 

Subroutine PRINTA 

Subroutine PRINTA echoes the user input or default values from the Master 
Control Input. 

Subroutine PRINTS 

Subroutine PRINTB echoes the user input or default values from the Mate- 
rial Control Input. 
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Subroutine PRINTP 


This subroutine prints the values of all of the statistical quantities 
calculated in MATL, OUTLIE, LEAST2 , MAXL, and MATBAT subroutines at each dis- 
crete temperature. A detailed description of the information printed can be 
found in section 5.0 PC-CARES OUTPUT INFORMATION, of this manual. 

Subroutine READINI 

This routine reads the initialization file, PCCARES.INI, from the direc- 
tory from which the PC-CARES program was executed and sets the INFILE, OUTFILE, 
IMAXF, and IMAXT parameters according to commands in the initialization file, 
if present. If the initialization file is not present then the default values 
of these parameters are used. 

Subroutine SNGLES 

This subroutine is used with the Batdorf surface flaw model to integrate 
over the contour of the unit circle when Shetty's failure criterion is used 
and a P > 01 . SNGLES determines the intervals where o 9 > o cr about the 

unit circle and stores the limits of these intervals as a function of a cr in 
the INTVAL array. The critical intervals correspond to the integral described 
by equation (A. 55) to determine to(Z,a cr ). The limits of these intervals are 
determined for each transformed Gaussian value of o cr . Each consecutive pair 
of integers in the second index of array INTVAL represents an interval where 
CT e 2 CT cr f° r a va 1 ue °f a cr denoted by the first index. The limits of 
integration stored in INTVAL are integers representing 1° increments of angle 
counted from -ir/2 to ir/2. 

Subroutine SORMAL 

This subroutine is called when the Batdorf surface flaw model is used 
with Shetty's failure criterion when a p _ > 01 . SORMAL calculates the nor- 

malized effective stress about the unit circle as a function of the angle a 
and stores it in the SEANGL array. The effective stress is determined for 1° 
increments of a. The array values are normalized by the maximum effective 
stress 0 P found for the stress state being evaluated. The effective 

stress is calculated from equations (A. 52) to (A. 54). 

Subroutine SORTRA 

This subroutine sorts the experimental fracture stresses at a given tem- 
perature level into ascending order (IASEND = 1) along with the corresponding 
fracture origins. It is invoked at the beginning of subroutine MATL. Array D 
contains the fracture stresses to be sorted, and the parameter NSORT equals 
the number of fracture stresses. The alphanumeric AINDEX array contains frac- 
ture origins (S, V, or U) that correspond in position to the sorted stresses. 

Subroutine SVALP3 

Subroutine SVALP3 is used with the Batdorf model for surface flaw analysis 
with Shetty's criterion when o e(nax > 0 ]. It performs the integration 
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which is used in equation (A. 55). Gauss-legendre quadrature is used for the 
numerical integration of doer- The stored values in the INTVAL array, which 
were previously calculated in the SNGLES subroutine, are used to perform the 
integration. The FINDP subroutine is called from SVALP3 to perform the evalu- 
ation of P 2 S for each transformed Gaussian value of o cr . 

Function UNBFTR 

This function performs a table lookup of the unbiasing factors for the 

estimated We i bu 11 modulus m. These factors are a function of sample size and 
are taken from reference 10. They are obtained from a Monte Carlo simulation 
of unimodal fracture data by using maximum likelihood analysis. The factors 
are based on the sample mean. The unbiased estimate of m is obtained by 
multiplying the biased estimate of m by the unbiasing factor. 
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APPENDIX 


SYMBOL LIST 

A surface area 

A-D Anderson-Darling 

A^ Anderson-Darling goodness-of-f i t test statistic 

a crack half length or penny-shaped crack radius or radius of semi- 

circular surface crack 

a-j coefficients of quadratic equation used for calculating <t>i where 

j = 1,2,3 

B risk of rupture in We 1 bu 11 1 s cumulative failure distribution 

C modified Neibull parameter (C = l/<70) m , or centigrade measure of 

temperature 

c contour of a unit radius circle in two-dimensional principal stress 

space 

C Shetty's constant in mixed-mode fracture criterion 

cm centimeter 

CV critical value 

Dj constants used in calculating ?2 for j = 1,2,3 

D + ,D - Kolmogorov-Smi rnov goodness-of-f i t test statistic 

D Kolmogorov-Smi rnov goodness-of-f i t test statistic defined as D + 

or D~ whichever is the largest 

E Young's modulus of elasticity 

EDF empirical distribution function 

F Fahrenheit 

F(x) cumulative distribution function of a random variable 

F(a) Weibull cumulative distribution of material strength 

Fnj(x) empirical distribution function 

G strain energy release rate, or crack extension force 

Gc critical value of strain energy release rate 


65 



Gl 

Gil 

Gill 

Gmax 

Gt 

GC 

GN 

GPa 

9 

h 

i 

i n . 

K I 

K IC 

KlI 

K III 

<& 

K(N) 

k B 

k B 

K-S 

k w 

k wp 

Ll 

L2 

lb 


strain energy release rate for crack opening mode crack extension 

strain energy release rate for crack sliding mode crack extension 

strain energy release rate for crack tearing mode crack extension 

maximum strain energy release rate 

total strain energy release rate 

Griffith crack (flattened elliptical cylinder) 

Griffith notch 

gigapascal 

gram 

total height of MOR bar with rectangular cross section 

ranking of ordered fracture data in statistical analysis or any 
counter 

inch 

opening mode stress intensity factor 
critical opening mode stress intensity factor 
sliding mode stress intensity factor 
tearing mode stress intensity factor 
<H or Km 

Kanofsky-Srini vasan confidence band factors 

Batdorf crack density coefficient or flaw distribution parameter 
normalized Batdorf crack density coefficient 
Kolmogorov-Smi rnov 

Weibull crack density coefficient, (l/o 0 ) m 

polyaxial Weibull crack density coefficient 

length between outer loads in four-point bending 

length between symmetrically applied inner loads in four-point 
bending 

pound 
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1 ,m,n 


direction cosines of oblique plane normal in principal stress space 
for the Cauchy infinitesimal tetrahedron 

In natural logarithm 

MLE maximum likelihood estimate 

MNR maximum normed residual 

MOR modulus of rupture or extreme fiber fracture stress 

MOR 0 characteristic modulus of rupture or extreme fiber fracture stress 

at which 63.2 percent of MOR bars will fail 

MPa megapascal 

m We i bu 11 modulus or shape parameter; also Batdorf crack density 

function exponent or flaw distribution parameter 

mm millimeter 

N number of MOR specimens at a given temperature 

NGP number of Gauss base points used in numerical integration 

N(ct) We i bu 1 1 crack density function or number of flaws per unit 

volume/area with strength <o in uniaxial stress state 

N(o cr ) Batdorf crack density function which is a material property 

independent of stress state and is the number of cracks per unit 
volume/area with strength £o cr 

n number of links in a structure 

n unit vector along oblique plane normal determined by angles a 

and |3 in principal stress space 

P load applied to MOR bar specimen 

Pf cumulative failure probability 

P$ cumulative survival probability 

PSC penny-shaped crack (flattened oblate spheroid) 

P^ probability of existence in incremental volume or area of a crack 

with strength between <a cr 

P2 probability of crack with strength <yj cr being so oriented that 

°e > °cr 

PIA principle of independent action 
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ps 1 pounds per square inch 

r number of remaining specimens in censored data analysis 

r^ Inside radius 

r 0 outside radius 

SC semicircular crack 

sec second 

t thickness 

V volume 

WLT weakest link theory 

w total width of MOR bar with rectangular cross section 

X ordered statistics 

x any variable 

x,y,z Cartesian coordinate directions 

Zi predicted failure probability at the fracture strength of the I** 1 

specimen 

a angle between o n and the maximum principal stress, a] 

(figs. A.l and A. 2 ); also significance level 

a defined as cos - V 5 , when oi > 02 =* <73 for volume flaws and also 

for surface flaws when o] > 02 

3 angle between c n projection and the Intermediate principal 

stress 02 in plane perpendicular to oj (fig. A.l) 

3 defined as root of cos - V 5 , when 02 * 03 

A Increment 

AP] probability of existence In Incremental volume or area of a crack 

with strength between o cr and a cr + A o cr 

r gamma function; tabulated In mathematical handbooks 

v material Poisson's ratio 

n usual product notation 

ir 3.1416 
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E 

a 


CT cr 


a e 


ae max 




°n 


CT n 


applied multidimensional stress state or summation notation 

applied stress distribution; also the traction or stress vector on 
oblique plane of Cauchy infinitesimal tetrahedron 

remote, macroscopic, uniaxial, normal fracture stress of a crack 

effective stress acting on a crack plane, o e = f(a n ,x) 

maximum effective stress for the particular stress state 

extreme fiber fracture stress in MOR bar test 

normal stress acting on oblique plane whose normal is determined by 
angles a and (5 (figs. A.l and A. 2)) 

normal stress averaged about a unit radius sphere or unit radius 
circle 


°o 


Weibull scale parameter or normalizing stress 


ct u Weibull location parameter or threshold stress 


a l , CT 2 > CT 3 principal stresses (o; > 03 > 03 ) 

(70 volume or area characteristic strength or characteristic modulus of 

rupture, MOR 0 . This is the stress or extreme fiber stress at which 
63.21 percent of the specimens will fail 


T 

<t> 

Q 


(0 


shear stress acting on oblique plane whose normal is determined by 
angles a and |3 (figs. A.l and A. 2)) 

defined as cos^j! or cos^a, depending on stress state, for which 
a e - a cr = 0 

solid angle in three-dimensional principal stress space for which 
a e 2 a cr 

angle in two-dimensional principal stress space for which 
a e 2 a cr 


Subscripts : 


B Batdorf 

cr critical 

e effective 

f failure, fracture 

g gage 
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I 

crack opening 

mode 

II 

crack sliding 

mode 

III 

crack tearing 

mode 

n 

normal 


P 

polyaxial 


S 

surface 


s 

survival 


V 

vol ume 


w 

Wei bul 1 



Superscript: 


estimated parameter 



APPENDIX 


THEORY 

The use of advanced ceramic materials in structural applications requiring 
high component integrity has led to the development of a probabilistic design 
methodology. This method combines three major elements: (1) linear elastic 

fracture mechanics theory which relates the strength of ceramics to the size, 
shape and orientation of critical flaws, (2) extreme value statistics to obtain 
the characteristic flaw size distribution function, which is a material prop- 
erty, and (3) material microstructure. Inherent to this design procedure is 
that the requirement of total safety must be relaxed and an acceptable failure 
probability must be specified. 

The statistical nature of fracture in engineering materials can be viewed 
from two distinct models (ref. 18). The first was presented by Weibull and 
used the "weakest link" theory as originally proposed by Pierce (ref. 19). The 
second model was also analyzed by Pierce (ref. 19) and, in addition, by Daniels 
(ref. 20) and is referred to as the "bundle" or "parallel" model. In this 
model, a structure is viewed as a bundle of parallel fibers. Each fiber will 
support a load less than its breaking strength indefinitely but will break 
immediately under any load equal to or greater than its breaking strength. 

When a fiber fractures, a redistribution of load occurs and the structure may 
survive. Failure occurs when all of the fibers have fractured. The weakest 
link model assumes that the structure is analogous to a chain with "n" links. 
Each link may have a different limiting strength. When a load is applied to 
the structure such that the weakest link fails, then the structure fails. 
Observations show that advanced monolithic ceramics closely follow the weakest 
link hypothesis. A component fails when an equivalent stress at a flaw 
reaches a critical value which depends on a fracture mechanics criterion, 
crack configuration, crack orientation, and the crack density function of the 
material. In comparison with the bundle model, WLT is in most cases more 
conservative . 

One of the important features of WLT is that it predicts a size effect. 

The number and severity of flaws present in a structure depends on the mate- 
rial volume and surface area. The largest flaw in a big specimen is expected 
to be more severe than the worst flaw in a smaller specimen. Another conse- 
quence of WLT is that component failure may not be initiated at the point of 
highest nominal stress (ref. 21), as would be true for ductile materials. A 
large flaw may be located in a region far removed from the most highly stressed 
zone. Therefore, the complete stress solution of the component must be 
obtai ned . 

Classical WLT does not predict behavior in a multiaxial stress state. A 
number of concepts such as the PIA, Weibull 's normal stress averaging method, 
and Batdorf's model have been applied to account for polyaxial stress state 
response. Batdorf's model (ref. 22) assumes the following: (1) microcracks in 

the material are the cause of fracture, (2) cracks do not interact, (3) each 
crack has a critical stress, o cr , which is defined as the stress normal to the 
crack plane which will cause fracture, and (4) fracture occurs under combined 
stresses when an effective stress, o e , acting on the crack is equal to o cr . 

For an assumed crack shape o e may be obtained through the application of a 
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fracture criterion. These concepts are applied in the PC-CARES code to obtain 
the normalized Batdorf crack density coefficient. 


A. 1 Volume Flaw Reliability 


Consider a stressed component containing many flaws and assume that 
failure is due to any number of independent and mutually exclusive mechanisms 
(links). Each link involves an infinitesimal probability of failure APfv- 
Discretize the component into n incremental links. The probability of sur- 
vival, P s v> the i th link is 


<P sV > i = 1 - <AP fV > i 


(A. 1 ) 


where the subscript V denotes volume dependent terms. The resultant proba- 
bility of survival of the whole structure is the product of the individual 
probabilities of survival 


n n 

p sv -n<vi -nt’- < 4 p fv>i] 

i =1 i»1 


n 


n 

n 

] = exp 

- E <ap f V >i 

i-i 


. 1-1 


(A. 2) 

Assume the existence of a power function N\/(ct), referred to as the crack 
density function, representing the number of flaws per unit volume having a 
strength equal to or less than a. In uniform tension of magnitude a, the 
probability of failure of the i th link, represented the incremental volume 
AVi, is 

( AP fv >. = N v <o)AV. (A. 3) 

and substituting into equation (A. 2), the resultant probability of survival is 

P sV = exp[-N v <o)V] (A. 4) 

and the probability of failure is 

P fv - 1 - exp[-N v <o)V] (A. 5) 

where V is the total volume. If the stress is a function of location then 
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fV ' 


1 - exp 



Ny(o)dV 


1 - expt-By] 


(A. 6) 


A term called the "risk of rupture" by Weibull and denoted by the symbol 
By is commonly used in reliability analysis. Equations similar to (A. 5) and 
(A. 6) are applicable to surface distributed flaws where surface area replaces 
volume and the flaw density function is surface area dependent. 
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Welbull Introduced a three-parameter power function for the crack density 
function Ny(cr), 


N y (a) 



(A. 7) 


where a u \/ is the threshold stress (location parameter), usually taken as 
zero for ceramics. The location parameter is the value of applied stress 
below which the failure probability is zero. When the location parameter is 
zero, the two-parameter Weibull model is obtained. The scale parameter, a 0 y, 
then corresponds to the stress level where 63.2 percent of specimens with 
unit volumes woulcj^fracture. The scale parameter has dimensions of 

stress x(volume) mv . my is the shape parameter (Weibull modulus) which 
measures the degreeof strength variability, my is a dimensionless 
quantity. As my increases, the dispersion is reduced. For large values of 
my (my > 40), such as those obtained for ductile metals, the magnitude of the 
scale parameter corresponds to the material ultimate strength. These three 
statistical parameters are material properties, which are temperature and 
processing dependent. 


Three-parameter behavior is rarely observed in as-processed monolithic 
ceramics and statistical estimation of the three material parameters is very 
Involved. Therefore, the PC-CARES program uses the two-parameter model. The 
subsequent reliability predictions are more conservative than for the three- 
parameter model since we have taken the minimum strength of the material as 
zero. 


The two-parameter crack density function is expressed as 


N v (o) 



(A. 8 ) 


and substituting equation (A. 8 ) into equation (A. 6 ), the failure probability 
becomes 
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expj -k 


wV 


J V 




(A. 9) 


where k w y = o 0 y mv is the uniaxial Weibull crack density coefficient. Var- 
ious methods have been developed to calculate <?qV and my for a given 
material using fracture strength data from simple uniaxial specimen tests 
(ref. 17). 


The two most common techniques for using uniaxial data to calculate Pfy 
In polyaxlal stress states are the PIA (refs. 23 and 24) and the Weibull normal 
tensile stress averaging method (refs. 25 and 26). In the PIA model, the prin- 
cipal stresses ai > a£ > 03 are assumed to act Independently. If all princi- 
pal stresses are tensile, the probability of failure according to this approach 
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(A. 10) 


Compressive principal stresses are assumed not to contribute to the failure 
probability. It has been shown that this equation yields nonconservative esti- 
mates of Pfv in comparison with the Weibull normal stress method (ref. 27). 

The failure probability using the Weibull normal tensi le. stress averaging 
method, which has been described through an integral formulation (ref. 28), can 
be calculated from 


where 




[ "V 

A a n 


dA 


dA 


(A. 11) 


The area integration is performed in principal stress space over the surface, 

A, of a sphere of unit radius for regions where a n , the projected normal 
stress on the surface, is tensile. The polyaxial Weibull crack density coeffi- 
cient is k WD v- The relationship between k wp v and k w v is found by equating 
the fai 1 ure probability for uniaxial loading to that obtained for the polyaxial 
stress state when the latter is reduced to a uniaxial condition. The result is 


k wpV = <2m V + ,)k wV 


(A. 12) 


Batdorf and Crose (ref. 22) proposed a statistical theory in which atten- 
tion is focused on cracks and their failure under stress. Flaws are taken to 
be uniformly distributed and randomly oriented in the material bulk. Fracture 
is assumed to depend only on the tensile stress acting normal to the crack 
plane hence, shear-insensitivity is inherent to the model. Subsequently, 
Batdorf and Heinisch (ref. 29) included the detrimental effects of shear trac- 
tion on a flaw plane. Their method applies fracture mechanics concepts by com- 
bining a crack geometry and a mixed-mode fracture criterion to describe the 
condition for crack growth. Adopting this approach, the PC-CARES program con- 
tains several fracture criteria and flaw shapes for volume and surface analyses 
(fig. 2.2). 

Consider a small uniformly stressed material element of volume AV. The 
incremental probability of failure under the applied state of stress, E, can be 
written as the product of two probabilities, 

aP fv (r,« cr .4V) - 4 P )v p 2V <A.i3> 
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where APiv is the probability of existence in AV of a crack having a criti- 
cal stress between o cr and o cr + Ao cr . As previously noted, critical stress 
is defined as the remote, uniaxial, fracture strength of a given crack in 
mode I loading. p£v denotes the probability that a crack of critical stress 
o cr will be oriented in a direction such that an effective stress, o e (func- 
tion of fracture criterion, stress state, and crack configuration) satisfies 
the condition o e 2 a cr- The effective stress is defined as the equivalent 
mode I stress a flaw would experience when subjected to a multiaxial stress 
state which results in modes I, II, and III crack surface displacements. 


Crack dimensions are related to crack strength, and crack size is never 
explicitly used in statistical fracture theories. Batdorf and Crose (ref. 22) 
describe APiv as 


4P, V - AV 


d W> 


da 


cr 


cr 


(A. 14) 


and P 2 V Is expressed as 
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2V “ 


Q<E -V > 

4-rr 


(A. 15) 


Where Nv<o cr ) is the Batdorf crack density function. Q(E,o cr ) is the area of 
the solid angle projected on the unit radius sphere in principal stress space 
containing all the crack orientations for which o e > o cr . The constant 4ir 
is the surface area of a unit radius sphere and corresponds to a solid angle 
containing all possible flaw orientations. 


The probability of survival in a volume element AVi is 
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(P sV ) = exp 2 -AV 1 


n a e 
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Q( E ,o ) dN.,(o ) 
cr V cr 


4ir 
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do. 


cr 


cr 


(A. 16a) 


where <*e max i s the maximum effective stress a randomly oriented flaw could 

experience from the given stress state. Hence, the component failure proba- 
bility is 
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(A. 16b) 


The Batdorf crack density function Ny(or r ) is a material property, 
independent of stress state, and is usually approximated by a power function 
(ref. 29). This leads to the Batdorf crack density function of the form 


Vv’ * k BV°cr 


(A. 17) 
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where the material Batdorf crack density coefficient kgy and Weibull modulus 
my can be evaluated from experimental fracture data. Batdorf and Crose 
(ref. 22) initially proposed a Taylor series expansion for N \/<a cr ), but this 
method has computational difficulties. A more convenient integral equation 
approach was recently formulated and extended to the use of data from four- 
point MOR bar tests (ref. 30). Note that N v (a cr ) has units of inverse volume. 

Although the Weibull (eq. (A. 8)) and Batdorf (eq. (A. 17)) crack density 
functions are similar in form, they are not the same. The Weibull function 
simply depends on the applied stress, o, and is the only term other than the 
volume necessary to calculate Pfy. The Batdorf function depends on the mode I 
strength of the crack, o cr , which is probabilistic and must be integrated over 
a range of values for a given stress state. Furthermore, to obtain Pfy, a 
crack orientation function, P2V» mus t be considered in addition to the density 
function and the volume. Finally, the Batdorf coefficient, k B y, cannot be cal- 
culated from uniaxial data until a fracture criterion and crack shape are 
chosen, in contrast to the Weibull coefficient, k w y, which depends only on the 
data itself. 

Assuming a shear-insensitive condition, fracture occurs when 
0 = a e > a cr , where a n is the normal tensile stress on the flaw plane. 
However, for a flat crack it is known from fracture mechanics analysis that a 
shear stress, x, applied parallel to the crack plane (mode II or III), also 
contributes to fracture. Therefore, for polyaxial stress states expressing 
the effective stress, a e , as a function of both on and t is more accurate 
than assuming shear-insensitivity. Batdorf and Heinisch (ref. 29) give effec- 
tive stress expressions for two flaw shapes using both Griffith s maximum ten- 
sile stress criterion and Griffith's critical coplanar strain energy release 
rate (Gj> criterion. Arranged in order of increasing shear-sensitivity, for 
the maximum tensile stress criterion the effective stress equations are 


a 


e 




(Griffith flaw) 


(A. 18) 



(Penny-shaped flaw) (A. 19) 


where v is Poisson's ratio. 

The total coplanar strain energy release rate criterion is calculated from 


G t = G i + G n + G ni 


(A. 20) 


where G is the energy release rate for various crack extension modes. In 
terms of stress intensity factors, the effective stress equation can be derived 
from (plane strain condition assumed) enforcing the condition Gj = Gq, where 
Gq Is the critical strain energy release rate. Thus, 
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(A. 21) 
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For a Griffith crack, assuming that modes I and II dominate the response 
with Kj = a n Vira and Kjj = x-s/ira, where 2a Is the crack length, we have 
from equation (A. 21 ) 


a 


e 



(A. 22) 


For a penn y-s haped crack at the critical po int on the crack periphery, we have 
Ki = 2 o n Va/ir and Kjj = [4x/(2 - v) ] v/a/tr (ref. 31) where a now is the 
crack radius. The resulting effective stress equation is 
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(A. 23) 


The equations given by Batdorf and Heinlsch consider only self-similar 
(coplanar ) crack extension. However, a flaw experiencing a mu 1 1 1 a x 1 a 1 stress 
state usually undergoes crack propagation initiated at some angle to the flaw 
plane (noncoplanar crack growth). Shetty (ref. 14) performed experiments on 
polycrystalline ceramics and glass considering crack propagation as a function 
of an applied far field mu ltiaxial stress state. He modified an equation pro- 
posed by Palani swamy and Knauss (ref. 13) to empirically fit experimental data. 
This multimodal interaction equation takes the form 



(A. 24) 


where K$ is either Kjj or Kjjj, whichever is dominant, and C is a 
constant adjusted to best fit the data. Shetty (ref. 14) found a range of val- 
ues of 0.80 < C < 2.0 for the materials he tested which contained large 
Induced flaws. As C Increases, the response becomes progressively more 
shear-insensitive. 

Using this relationship with assumed modes I and II dominance for the 
Griffith crack yields 
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e 




(A. 25) 


and for a penny-shaped crack, we get 
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To determine a component probability of failure from equation (A. 16), P 2 V 
has to be evaluated for each elemental volume AV within which a uniform 
stress state 1 ( 01 , 02 , 03 ) is assumed. The solid angle Q(E,o cr ) depends on 
the selected fracture criterion, crack configuration and on the applied stress 
state. For multiaxial stress states, with few exceptions, Q(E,o cr ) must be 
determined numerically. For a sphere of unit radius (fig. A.l), an elemental 
surface area of the sphere is dA = sin a d3 da. Project onto the spherical 
surface the equivalent stress o e (E,a,(3). The solid angle Q(E,o cr ) is the 
area of the sphere containing all the projected equivalent stresses where 
CT e > °cr- Noting the symmetry of o e , and addressing the first octant of the 
unit sphere, then 


Q(E,o cr ) = 4 ttP 2V = 8 


p ir/2 

J 0 



sin a da 


(A. 27) 


where (5 is evaluated between 0 and ir/ 2 . 

To obtain the limits of integration, 3i and 32 , f° r the interval where 
a > a cr , the principal stresses must first be transformed to normal and shear 
stresses. Selecting an arbitrary plane and imposing equilibrium of forces 
(fig. A.l), the following equations are obtained: 


where a is 
tion cosines 
functions of 
figuration 


2 2 n 2 22 22 
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(A. 28) 
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(A. 29) 

2 2 2 
t = a - a n 



(A. 30) 

total traction vector acting on the 
m, and n are given in figure A.l 
and 3. From the selected fracture 
s obtained as a function of E, a, 

crack plane and the direc- 
in terms of trigonometric 
criterion and crack con- 
and 3 - 


By defining $ = cos 2 3 and enforcing the failure condition of a e = o cr , 
we obtain a quadratic equation in 4* satisfying either 
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cr 


(°n * Ih2 ) ' 0 


(A. 31) 
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where D Is some constant defined by the specific fracture criterion and 
crack geometry. Equation (A. 31) Is used with the effective stress equa- 
tions (A. 22) and (A. 23). Equation (A. 32) Is used when the effective stress 
equations (A. 18), (A. 19), (A. 25), and (A. 26) are selected. The quadratic 
equation takes the form 



+ a 2 <t> + a 3 = 0 


(A. 33) 


and the roots and 4> 2 are 





- 4a l a 3 


2a 


1 


(A. 34) 


where <J>i < $ 2 - The expressions for coefficients ai to a 3 are given In 
tables A. I and A. II. The values for 3 are then found as 


3, = cos' 1 ^ 


^ = 0 
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0 < $ 2 < ! 




*2 < 0 or 4>2 > 1 or 4> 2 Isa complex number 


32 = cos -1 y'fFJ" 0 < < 1 
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(A. 35) 


< 0 or 4>i > 1 or 4>^ Isa complex number 


J 


After obtaining 3] and 32 for a given stress state E, a, and a cr , 
care must be taken In evaluating the Integral. The solution of the integral 
In equation (A. 27) Is either 


p3i 


„ir/2 


d3 + 


d3 = 3, - 3 2 + f 


(A. 36) 


or 


d3 - 3 2 - 3, (A. 37) 

^1 

The correct solution Is determined by checking If o e - o cr > 0 at some 
angle 3 between 0 and ir/2. In the PC-CARES program extensive logic has been 
devised to examine all possible permutations the roots may have. Including 
Imaginary roots. 
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An alternative approach to calculate is to Increment the angles a 

and (3 over the surface of a unit radius sphere. By symmetry only one octant 
needs to be considered. At each discrete point on the surface, the effective 
stress Is evaluated and the associated area element Is summed depending on 
whether a e > a cr . This procedure Is computationally Intensive, and, whenever 
possible PC-CARES employs the more efficient approach described previously. 


For a given stress state and value of a cr j , a is varied from 0 to ir/2 
and £KZ,a cri ) Is evaluated. The values of a cri vary from 0 to o emax , for 
the Gauss-Legendre Integration used by PC-CARES. The probability of survival 
In volume AV^ Is obtained by substituting equation (A. 17) Into equation 
(A. 16a) to get 


r 




^sV^I = ex P % ~ AV i m v k BV 


n°e 


max 


Q(E,cr ^.) m^-1 

: a da 

4ir cr cr 


(A. 38a) 


and the component failure probability Is 
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(A. 38b) 




Consider the simple stress state o] > 02 - 03 . For this case o e and 
Q(E,a cr ) are Independent of (3 and equation (A. 27) reduces to 


Q(E,a cr > = 8 


17 r W2 A r 

d|3jsin a da =» 4ir J sin a 


da 


(A.39) 


where a is Integrated between 0 and ir /2 In a manner similar to the Integra- 
tion of |3 in equation (A. 27). The quadratic equation (A. 33) Is reformulated 
as a function of a, where now 4> = cos^a. Table A. 11(a) contains the coeffi- 
cients a] to a 3 for calculating and thence Q(E,a cr ) for various frac- 
ture criteria and crack shapes for this stress state. The logic for evalua- 
ting the a Integral Is the same as that for the |3 Integral, as described 
In equations (A. 35) to (A. 37). With the possible exception of the Shetty cri- 
terion or when 03 < 0 , the quadratic equation can have only one root between 
0 and 1 and Q<Z,a cr ) is simply 4ir<l - cos a) where a corresponds to the 
single root. If both roots lie outside the range 0 to 1 , then a sample point 
Is required to determine whether fXE.acr) * 0 or £KE,a cr ) = 4ir. Additional 
equations for calculating P 2 V are also listed In table A. 11(b) for special 
stress states, such as the uniaxial, equl biaxial and equltrlaxlal loading 
condl tlons . 
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For certain stress states and crack plane orientations, the normal stress 
on the crack plane can be compressive. When this situation occurs in the 
PC-CARES program the normal stress is set to zero and only the shear stress is 
assumed to contribute to crack growth. This is generally a conservative 
assumption since friction between the crack faces is ignored. If friction were 
considered, the effective applied shear would be reduced. 

For most fracture criteria, c em;iv = a], that is the maximum effective 

stress is equal to the maximum tensile principal stress. For_noncopl anar crack 
extension using equations (A. 25) and (A. 26), if 1/C and 2/C(2 - v) are <1.0, 
respectively, then o emax = oi . If these terms are greater than 1, then 

° e max > °1 P oss ^b^ e - Also a e(liax > a] is possible when 03 < 0. For these 

conditions, the values of 3 of equation (A. 35) are found by a surface element 
sampling scheme. 

For the special case of shear insensitivity, the projected equivalent 
stress on a unit radius sphere is equal to the normal stress, that is, 
o e = a n . Substituting for o n , we obtain 

2 2 P 

a e = °3 + (o l ~ a 3 )cos a + (a 2 - a 3 )cos |3 sin a (A. 40) 

The value of 3 satisfying o e - a cr = 0 is obtained by defining 

* = cos^3 and calculating the coefficients ai for equation (A. 33). For 
this shear-insensitive case, we get 
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(A. 41 ) 


and 


a, = (a, - ct-,)C 0S a + o~. - a 
413 3 cr 


We can now solve for <t> to obtain 


(A. 42) 
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( o ^ - a^) S i n a 


(A. 43) 


It is obvious that only one value of 4> satisfies equation (A. 43), from which 
the limits of integration become 
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= o 

i! 2 = cos -1 sfi if 

3 2 = 0 if 

h-T if 

and equation (A. 37) Is used for all cases. 
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0 < $ < 1 

<t> > 1 

<t> < 0 
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(A. 44) 


A. 2 Surface Flaw Reliability 

For surface flaw analysis (ref. 2), many of the equations from section 
A 1 Volume Flaw Reliability remain the same, except that the statistical mate- 
rial parameters are a function of surface area instead of volume and the equiv- 
alent stress projections are onto the contour of a circle of unit radius rather 
than onto the surface of a unit radius sphere. The cracks are assumed to be 
randomly oriented in the plane of the external boundary with their planes nor- 
mal to the surface. 

For surface flaw induced failure in ceramic structures the probability of 
failure for the two-parameter Weibull distribution, which is analogous in form 
to equation (A. 9) is 
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fS ‘ 



(A. 45) 


where k w c = (l/aoc)" 15 , is the Weibull surface crack density coeff 
subscript S denotes the terms that are surface area dependent. 

the surface scale parameter with units of stress x(area) s and 
stressed surface area. For biaxial stress states, the Weibull dis 
combination with the PIA hypothesis yields 


i dent. The 
Here a 0 $ is 
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(A. 46) 


where ai and a? are the principal tensile in-plane stresses acting on the 
surface of the structure. The failure probability using the Weibull normal 
stress averaging method can be calculated from 
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(A. 47) 
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where 



Here k W nS Is the polyaxlal Weibull crack density coefficient for surface 
flaws. The line Integration is performed over the contour, c, of a unit radius 
circle where the projected normal stress, c n , is tensile. The relationship of 
£wpS To k w $ is obtained by carrying out the Integration In equation (A. 47) 
for a uniaxial stress and equating the resultant failure probability to that of 
equation (A. 45) (ref. 17). This results in 



m^T(m^) Vir 

r ( m s ♦ I) 



(A. 48) 


where r is the gamma function. Equation (A. 47) is the shear-insensitive case 
of the more general Batdorf polyaxlal model. 

For mixed-mode fracture due to surface flaws the Batdorf polyaxial failure 
probability equation (analogous to eq. (A. 16b)) is 
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(A. 49) 


where 


AP 


IS 


d N s (a ) 

AA - ^ Cr da 
da cr cr 


and 


For randomly oriented cracks w(Z,o cr > Is the total arc length on a unit radius 
circle in principal stress space on which the projection of the equivalent 
stress satisfies a e 2 <*cr ar| d 2ir is the total arc length of the circle. 

The same as for volume flaws, the Batdorf surface crack density function is 
approximated by the power function, 
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(A. 50) 
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m$ 

k BS a cr 


where kgs is the Batdorf surface crack density coefficient. 

Fracture occurs when the equivalent stress ag >. ^cr- For the shear- 
insensitive case fracture depends only on the value of the normal tensile 
stress such that c e = <*n- For shear-sensitive cracks and col 1 near crack 
extension (Gj criterion), assuming a Griffith crack with Kj = a n Vwa and 
Kji = x Vira we obtain as before 


V* 


+ X 


(A. 51a) 


while for a Griffith notch subjected to plane strain conditions with 
= 1.121 5o n ^/xra and Kjn = x Vrra (ref. 31) we get 


a e = 
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2 0.7951 2 

a n + (1 - v) x 


(A. 51b) 


Note that the equivalent stress for the Griffith crack is dependent on modes I 
and II, while for the Griffith notch, the equivalent stress is dependent on 
modes I and III (ref. 2) . 

For noncoplanar crack growth, from equation (A. 24) the effective stress 
equations for the Griffith crack and Griffith notch, respectively, are 




(A. 52) 


and 


°e = 2 




+ 3.1803 


(A. 53) 


For a semicircular surface crack Kj = 1.366o n Va, <n - l-241xVa, and 
Kttt = 0.133xVa (refs. 32 and 33). Since the contribution of K nI is smal 
it is neglected, and thus, the effective stress for this case is 


a e = 2 
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o 2 + 3.301 — 


(A. 54) 
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For the same stress state and Identical C, the Griffith crack is the most 
shear-sensitive, while the Griffith notch and the semi-circular crack give 
almost identical predictions. 

The solution procedure for c o(E,a cr ) is similar to the methods outlined 
In section A.l Volume Flaw Reliability. The probability that the crack orien 
tatlon is such that a e > a cr can be calculated from 
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2S " 


to(E ’ g cr ) 

2-rr 



(A. 55) 


where over the unit radius circle, 0 < a, < ir/2. The limits of integration 
«1 and «2 are obtained through the enforcement of the failure condition 
°e = a cr- The required normal and shear stresses are calculated from force 
equilibrium on a crack plane. As shown in figure A. 2, the stress vector a, 
the normal stress a n and the shear stress -z can be expressed as 


2 

X 


2 

a 



- Or 


“ O n 


- a. 


1 2 2 

fcos a + a 2 

(A. 56) 
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(A. 57) 

\2 2 /. 2 \ 

J COS all - COS a) 

(A. 58) 


^ P °" S ^ t1tut1on of a n* x and satisfaction of o e = o cr , equations (A. 51) 
to (A. 54) are reduced to a quadratic expression of the same form as equa- 
tlon (A. 33) with $ = cos 2 a. However, since 3=0° (fig. A. 1(a)) In the 
oj - 02 plane, the constants ai to a 3 are dependent only on the two prln- 
C ' Stresses, o cr , and In some cases on Poisson's ratio. Using the solution 
methods outlined in the previous section we obtain the roots of the quadratic 
equation. These values are in table A. Ill along with the coefficients ai . 

For cases where the roots, 4^, of the quadratic equation are not between 0 
and 1, the calculation of P 2S in equation (A. 55) follows the same logic as 
has been given in equations (A. 35) to (A. 37), with a replacing 0. Specific 

examples for this situation have been given in reference 2. For the equibi- 

axlal surface stress state, we always have w/ 2ir = 1 for o cr < OQ , since 

the In-plane shear stress is zero and hence o e = oi for all values^of a and 

any effective stress equation. 


When the normal stress is compressive (a n < 0) It Is equated to zero 
and the shear stress alone contributes to crack, growth. The maximum egulv- 
alent stress a e(n ax Tor most cases is equal to a], the maximum principal 
stress. However, for noncoplanar crack extension, using equation (A. 24), 

° e max 1s dependent on the value of C, and may exceed a]. Also when 02 < 0, 
then cFg max > a] is possible. Again a sampling scheme Is used to evaluate 
w(E,a cr ) when this condition occurs. 
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A. 3 Estimation of Statistical Material Strength Parameters 

Selected statistical theories and equations for parameter estimation are 
explained in detai 1 in reference 17. The following Is a brief description of 
these methods and how they are used in the PC-CARES code. Typically for bn t- 
Hp materials the Weibull parameters are determined from simple specimen geom 
r/ IB' loldin conations, such as beams under flexure and etercy in real 
or flat sDecimens under uniform uniaxial tension. The flexural test failure 
probabmty can be expressed in terms of the extreme fiber fracture stress, of, 
or modulus of rupture, MOR, using the two-parameter Wei bull form as 


P f = 1 - exp(- 


Ca f) = 1 


exp 




C(MOR) m ] 


= 1 - exp 


MOR 

V a 0 y 


(A. 59) 


where m i s the volume or area Wei bul 1 modul us , C is the modified Wei bull 
narampter (C - (l/ae) m ) and oa is the volume or area specimen characteristic 
o cSa ritJc modulus of rupture, M0R o For uniform unaxiate- 
sion tests of in equation (A. 59) would just be replaced by 01 . The We bull 
- r;i 1 e parameter o n . as defined in equations (A. 7) and (A. 45) for volume and 
surface cracks respectively, is determined from o 0 , m, the specimen geometry, 
and the loading configuration. The scale parameter, a 0 , is based on a unit 
volume or area whereas a 0 includes the effects of the specimen dimensions. 
Th^charactertstlc strength c 6 .1; defined as the uniform stress or extreme 
fiber stress at which the probability of failure is 0.632. 

Before computing the estimates of the statistical material parameters, it 
is essential to P carefully examine the available specimen data to scree [J.^ e '" 
for oStlle?s. Very often, a data set may contain one or more values which may 

outl1ers n a t t different r signtf?cance°levels e are a explalried P in C references 6 andU. 

tion and censored statistics is only approximate. 

Various methods are available to estimate the statistical material param- 
eters from experimental data for the two-parameter Wei ^ \ distnbut on^ lhe 
success of the statistical approach depends upon how well the den 

sity function fits the data. Two popular techniques used to es 

rharacteristic strength and shape parameter (o 0 and m) are the least scares 
analysis and the maximum likelihood method. Least-sguares analysis is a spe- 
?U1 clle of the maximum likelihood method where the error is norm y 
distributed and has a zero mean and constant variance. The leas 9 
method i s not suitable for calculating confidence intervals and ^biasing 
factors, which quantify the statistical uncertainties in the available data. 

Equation (A. 59) can be linearized by taking the natural logarithm twice 
yielding 


an 



= an 


an 


1 - Pd 


= an C + m an a. 


(A. 60) 


86 



For the least-squares analysis, it is necessary to obtain the line of best fit 
with slope m and an intercept b which, as seen in equation (A. 60), is 
equal to the natural log of C. The failure probability Pf is determined by 
conducting fracture tests on N specimens. The fracture stresses are ranked 
such that of] < of 2 < ■ ..< of] < ... < ofN- For median rank regression 
analysis, the probability of failure of a specimen with rank i is 


i - 0.3 
N + 0.4 


(A. 61 ) 


By taking the partial derivative of the sum of the squared residuals with 
respect to m and C, and by equating the derivates to zero, values of m 
and C can be estimated. 


With censored data, one cannot directly use the median rank regression 
analysis as given in equation (A. 61) because of the competing failure modes. 
To take into account the influence of the suspended items, Johnson (ref. 9) 
developed the rank increment technique. For this technique, all observed 
fracture stresses are arranged in ascending order, and rank increment values 
are calculated for each failure stress from the following equation: 


Rank increment = s t + 1) - (previous adju sted rank) 

1 + (number of items beyond present suspended item) 

In the PC-CARES program for volume flaw analysis, all fracture stresses desig- 
nated as V's are considered as failure data; for surface flaw analysis the 
u £ r . e cons1dered . a s failure data. The new adjusted rank values are obtained 
by adding the rank increment value to the previously adjusted rank. These 
adjusted rank values and the median rank regression analysis (i.e., eq. (A. 61)) 
are then used to calculate the failure probability Pf. Finally, the Weibull 
parameters m and C are obtained. 


i-i d ' s tn bution of errors from the data is not normal, the maximum 

likelihood method is often preferred in Weibull analysis. This method has 

inherent properties. The likelihood equation from which the maximum 
likelihood estimates (MLE's) are obtained will have a unique solution. In 
addition, as the sample size increases the solution converges to the true val- 
ues of the parameters. Another _ feature of the maximum likelihood method is 
that there are no ranking functions or linear regression analysis when complete 
or censored samples are analyzed. The likelihood equation for a complete sam- 
ple is given by 



The values of m and ctq which maximize the likelihood function L, are 
determined by taking the partial derivative of the logarithm of the likelihood 

function with respect to m and with respect to oq. The values of m and 

OQ are obtained by equating the resulting expressions to zero and solving the 
simultaneous equations using the Newton-Raphson iterative technique. The MLE 

of m and oq are designated by m v and S ev and by me and S 0 e for 
volume flaw analysis and surface flaw analysis, respectively. For censored 
statistics we have 
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where r is the number of remaining specimens failed by the flaw mode for 
which parameters are being calculated. For a complete (uncensored) sample, r 
is replaced by N which is the total size of the sample. 

The MLE of the shape parameter, m, is always a biased estimate that 
depends on the number of specimens in the sample. Unbiasing of the shape 
parameter estimate is desired to minimize the deviation between the sample and 
the true population. The unbiased estimate of m is obtained by multiplying 
the biased estimate with an unbiasing factor (ref. 10). The confidence inter- 
vals for complete samples can also be obtained (ref. 10). For censored sam- 
ples a rigorous method for obtaining confidence intervals has not yet been 
developed due to the complexity of competing failure modes. Confidence bounds 
for censored statistics are instead estimated in the PC-CARES code from the 
factors obtained from complete samples (ref. 17). Confidence bounds enable 
the user to estimate the uncertainty in the parameters as a function of the 
number of specimens. Bounds at 90-percent confidence level and therefore, 

5 and 95 percentage points of distribution of the MLE's of the parameters, have 
been incorporated into the PC-CARES program, with data taken from reference 10. 

Subjective judgement is needed to test the goodness-of-f i t of the data to 
the assumed distribution. When graphical techniques are used, it can be very 
difficult to decide if the hypothesized distribution is valid, especially for 
small sample sizes. Therefore, many statistical tests have been developed to 
quantify the degree of correlation of the experimental data to the proposed 
distribution. 

In general, a statistic is a numerical value computed from a random sample 
of the total population. The difference between an empirical distribution 
function (EDF) and a hypothesized distribution function is called an EDF sta- 
tistlc There are two major classes of EDF statistics and they differ in the 
manner' in which the functional (vertical) difference between the EDF and the 
proposed distribution function F(x) is considered. The Kolmogorov-Smi rnov 
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(K-S) goodness-of-f i t statistic D belongs to the supremum class and is very 
effective for small samples. It uses the largest vertical difference between 
the two distribution functions to determine the goodness-of-f i t . For the K-S 
test, the sample is arranged in ascending order, and the empirical distribution 
function Ffg < x ) is a step-function obtained from the following expressions: 


F N (x) = 0 

x < X 1 


> 




X, < x < X. . 
1 - i + l 

i = 1,2,. . 

, . , N — 1 


(A. 66) 

F n (x) = 1 

X 

I'N 

X 






where X] < X 2 < . . . X} . . . < are the ordered fracture stresses from a 
sample of size N . The statistic D is obtained by initially evaluating two 
other statistics D + and D _ , the largest vertical differences when Fm(x) is 
greater than F(x) and the largest vertical differences when F|^(x) is smaller 
than F(x) , respectively. All three statistics are calculated by using the 
following expressions: 


D + 


h - F(x). 




1 


D = 


F«>t -Sr 1 


1 = 1,2 N 


D = max(D + ,D~) 


J 


(A. 67) 


For ceramics design, the F(x)j's are equal to Pf's and are calculated bv 
using equation (A. 59) . 


On the other hand, the Anderson-Darling statistic, A 2 , belongs to the 
quadratic class and is a more powerful goodness-of-f i t statistic. It evalu- 
ates the discrepancy between the two distributions through squared differences 
and the use of an appropriate weighting function. The statistic A 2 is qiven 
by 


A 2 = 


-N - 


N 

BE 

i + l 


(21 - 1 ) 


{“" Si )] ♦ tn [' - 


(N+l-i ) 


(A. 68) 


In this case, Z}'s are the predicted failure probabilities obtained from equa- 
tion (A. 59). Corresponding significance levels a are calculated from the D 
and A^ statistics. From previous surveys (ref. 17) there is no specific men- 
tion of an absolute accepted significance level. Therefore, the user has to be 
subjective, using his own judgement in either accepting or rejecting the 
hypothesi s that the data fit a Weibull distribution. However, a higher value 
of a indicates that the data fit the proposed distribution to a greater 
extent. 


For complete samples, the 90-percent Kanofsky-Srini vasan confidence band 
values about the proposed distribution are also calculated to ascertain the fit 
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of the data. These values are similar to the K-S statistic D centered around 
the EDF . The bands are generated by 

Confidence bands = [ F ( x ) - K(N), F(x) + K(N)] (A. 69) 

where F(x) Is the failure probability obtained by substituting the Weibull 
parameters in equation (A. 59). The Kanofsky functions, denoted by K(N), are 
described in reference 34. 

Some limitations are intrinsic to a purely statistical approach to design. 
One problem occurs when the design stress is well below the range of experimen- 
tal data as shown in figure A. 3. Extrapolation of the Weibull distribution 
into this regime may yield erroneous results if other phenomena are present. 
When two flaw populations exist concurrently, but only one (population A) is 
active in the strength regime tested, the predicted failure probability may be 
incorrect. Furthermore, if the threshold strength is not zero, the strength 
may be underestimated. Finally, an approach based only on statistics can allow 
for stress state effects only in an empirical fashion. 


A. 4 Material Strength Characterization 

Ceramic strength is an ambiguous entity since, for brittle materials, ten- 
sile strength, compressive strength, shear strength, flexural strength, and 
theoretical strength all have unique meanings and different values. The theo- 
retical strength Is defined as the tensile stress required to break atomic 
bonds which typically ranges from one-tenth to one-fifth of the elastic modu- 
lus for ceramic materials. Because of processing flaws, this strength is never 
obtained. A much more meaningful strength measurement Is the tensile strength 
In uniaxial tension or through flexural testing. In flexural strength testing 
the bend strength of of a ceramic Is defined as the maximum tensile stress in 
the extreme fiber of a beam specimen (modulus of rupture, MOR). The main 
objective of the PC-CARES program Is to characterize ceramic strength in terms 
of the MOR or pure uniaxial strength and to use this data with appropriate 
analysis to predict component response under complex multlaxial stress states. 
The PC-CARES program calculates required polyaxlal statistical material 
strength parameters from uniaxial tensile specimen or four-point bend specimen 
fracture data. After evaluating the Initial parameters as described in sec- 
tion A. 3 Estimation of Statistical Material Strength Parameters, additional 
calculations are performed to determine the material scale parameter and 
Batdorf crack density coefficient for use In the reliability calculations. 

For volume flaw analysis using four-point MOR bar data with known geometry 
(fig. A. 4), the value of Cy in equation (A. 59) are obtained from the least- 
squares or maximum likelihood analysis. The tensile stress distribution In a 
four-point bend specimen Is 
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4 xyo. 


a x " (L, - U)h 


2ya- 


CT x = h 


4 ( L ^ - x)ycy 


a x ~ <L 1 - l_ 2 >h 


0 < x < 
L, - L 


L 1 - 4 




2 

L, + L 


2 L 1 + L 2 
< x < = — - 


> 


< x < L 1 


(A. 70) 


By equating the risk of ruptures of equations (A. 9) and (A. 59), we obtain 

m, 


r !x_y v 

, CT ovJ 


m v 

dV = CyOy 


(A. 71) 


V 


and, after Integrating over the tensile portion of the bar, the scale param- 
eter Is 

l/n V 


oV 


wh (L 1 + m V L 2 > 


2 C v (m y + l) 2 


1 /m. 


(A. 72) 


where V e Is the effective volume. For uniaxial tensile loading, the effec- 
tive volume is equal to the gage volume V g , which is the uniformly stressed 
region where fracture Is expected to occur. Note when L 2 Is zero the solution 
for the three-point MOR bar Is obtained. 

For the Batdorf model, using the shear-insensitive case from 
table A.II(b), we have 


n ( r,a cr ) 

4ir 


= 1 - 



(A. 73) 


From equations (A. 16b) and (A. 17), after performing the do cr Integration, the 
risk of rupture for the four-point-bend specimen Is 


o°e_ 


'max 


a<E,a cr ) 

4tt 


' „ V 1 

m V k BV a cr 


da 


cr 


V L 


dV = 


BV 


2m^ + 1 


m v 

a x dV 


(A. 74) 


Equating risk of ruptures from equations (A. 74) and (A. 59) gives 
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(A. 75) 


m V k BV 
C V°f = 2m v + 1 



d V 


from which, after Integrating the stress over the tensile loaded volume with 
o x defined In equation (A. 70), we get 


k 


BV “ 


<2m v + 


1) 


2C v <m v + l) 2 
wh(L^ + myl^) 



(A. 76) 


Using equation (A. 72) and the previously defined Welbull crack density coeffl 
clent, we have 


r £ k V 

C„ = m e 


(A. 77) 


M v 


Substituting equation (A. 77) Into equation (A. 76) and rearranging gives the 
normalized Batdorf crack density coefficient for the shear-insensitive ca , 


BV 0m ^ , 
BV = k wV = 2m V 


(A. 78) 


For the Batdorf shear-sensitive case, assuming a Griffith an J rom 
coplanar strain energy release rate criterion, we obta n Q<S.ocr from 
table A.II(b). For uniaxial loading o, after performing the Indicated t 

gratlon we get 


J V 


,a Q(Z.o^) dN x/ (o rr ) 
0 


4tr 


cr' -V-c r_ d<) 


da 


cr 


cr 


d V 


'BV 


"V 


"V 

a V dV 


'V 


(A. 79) 


Again equating the risk of ruptures from equations (A. 59) and (A. 79) In terms 
of the effective volume gives 



'BV 


m v + 1 



Using equation (A. 77), we substitute for C v to get 

k BV 

k wV = m v + 1 


(A. 80) 


(A. 81) 


from which the normalized crack density coefficient for the selected shear 
sensitive case Is 
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(A. 82) 


r _BV , 

k BV " k .. = m V + 1 
wV 

In the PC-CARES program, kgy is computed numerically for a shear- 
sensitive material for the general case where no closed-form solution exists. 
By using equations (A. 16b) and (A. 17) and equating the appropriate risk of 
ruptures we obtain 



then substituting for a x and performing the integration in equation (A. 85), 
the scale parameter is obtained as 


93 



r r 


a oS " 1 



+ D‘ 


1 /m. 


' m c w \ 

— H + 1 <w 1- h)l, > 
w + h I 1 


J 


(A. 87a) 


or 


a oS 



(A. 87b) 


where A e Is the effective area. For uniaxial tensile loading, the effective 
area Is equal to the specimen gage area Ag, which Is the total specimen sur- 
face area of Interest. Note when i _2 Is zero the solution for the three-point 
MOR bar Is obtained. 

For surface flaw reliability analysis with the Welbull normal stress aver- 
aging method, we calculate the polyaxlal crack density coefficient k^pj from 
the following equation (refs. 17 and 28): 



m<. si* r(m s )k wS 

r (”s * ?) 


(A. 88) 


where k w $ has been previously defined In equation (A. 45). 

By combining equations (A. 49) and (A. 50) for the Batdorf surface flaw 
model , we can express Pfs as 


r 

> 

p 

< -m s k BS 




A 



“max 


<o(E,o ) m c -l 
— lT~ V d °cr 


dA 


J 


(A. 89) 


For uniaxial tension with a shear-insensitive fracture criterion, substituting 
for &)(E,a ct -)/2ir from table A. Ill (02 = 0), we obtain 


P 


fS 


1 - exp < - 


2m S k BS 


i>0i 


COS 



dA > 


(A. 90) 





J 


Equating the risk of rupture In equation (A. 90) with that of equation (A. 45) 
results in 
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BS - k 


kgs m^ r(m^) 


ws r 


( m s * ?) 


Wo^hni! S S f eci f 1 1 case , the Batdorf crack density coefficient is identica' 
to the Welbuil polyaxlal crack density coefficient: that is 

k BS = k wpS <A.921 

Similar results were obtained for volume flaw based analysis as well, 
nn general shear-sensitive case, kgs is computed numerically since 

lUMM'X .l9;% t Wes ThUS - equat,n9 the r,Sk 0f rUBtUres * 


C S a f = m S k BS 


) m^-1 

U °cr do cr dfl 


from which we obtain 



(A. 94) 
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AREA OF PLANE-A 


°3 




DIRECTION COSINES: 


n 


cos a 

sin a cos & 
sin a sin (S 


(a> 



(b) 


FIGURE A. 1 . - STRESSES ON CAUCHY INFINITESIMAL 
TETRAHEDRON IN PRINCIPAL STRESS SPACE (8>, PRO- 
JECTED ONTO A TANGENT PLANE TO THE UNIT RADIUS 
SPHERE AS SHOWN IN <b). 





DESIGN 

STRESS, 

°d 


! A 


DATA 
f RANGE 


) 


i / 
I / 

/ 


/ w EXTRAPOLATION RANGE 


/ 


14 


BASIC DATA 



NONZERO THRESHOLD STRENGTH 



UNDETECTED FLAW POPULATIONS 



FIGURE A. 3. - LIMITATIONS OF EXPERIMENTAL DATA EXTRAPOLATION AND 
STATISTICAL APPROACH TO DESIGN. 
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TABLE A. I. - FORMS OF P 2V FOR VARIOUS FRACTURE CRITERIA AND SELECTED CRACK CONFIGURATIONS 

02 ) o 3 , o 3 >. 0 and Og - o^ r .] 



where 4^ i *2 


5 (a,o ) = cos */*? 

1 cr f 2 -1 r 

jj(a,a ) = COS Y# 

cr * 

0 2 <a ’°cr ) = C ° S_1 


After obtaining roots for a given stress state E, o cr and varying a, care must be taken in 
evaluating j d|3. The relation of o e to o cr in the neighborhood of S must be known to 
obtain the proper limits of the integral. 




TABLE A. I. - Concluded. 


criterion configuration 


Strain 
energy 
rel ease 
rate, 

Gt 


Quadratic equation coefficients 
for a n (E,a,|3) > 0 


Gri f f i th a = 0 

crack 


= “(°2 ” a 3 ) Sin2 “ 


a 2 2 2.2 2 
= -<*. COS a - a sin a + o 
J I 3 cr 


Penny-shaped 

crack 

(PS C) 


a, = D 


:(°2 ' “5) 


2 . 4 

sin a 


2 . 2 \2 2 
y cos a + a si n a ) + a 

1 3 / cr 


Gri f f i th 

crack 

(GC) 


2 = D l (°2 * a + ^° 2^°2 ~ a 3^ S ^ n ^ a (^3 sin ^ a + ° 1 cos 2 a ^ 

_ n / 2 2 2 . 2 \ / 2 2 \2 2 

3 D l\l C ° S a + a 3 Sina y + ^2^1 C0S a + s i n + a 


J_ ( \2 . 4 

1 " -2 \°2 " ° 3 / Sin a 

2 = ^°2 ” a 3 ^ S ^ n a |“ 0 cr + -2 (^ a i C ° S a + ^°3 s ^ n2<3t “ a 2 “ a ^) 
i = °cr “ °cr( a i C ° s2a + °3 sin2<l ) ’ ^ [(°, " ° 3 ) 2 si " 2 « «* 2 “ 


a 2 = (°2 " ° 3 ^ sin 


2 o • 2 

4a COS a + 2o^ sin a - a - a 
\ 1 3 2 3 


a = a 
3 cr 


Penny-shaped 

crack 

(PSC) 




2 . 4 

I si n a 


2 

* 2 = ( a 2 “ a 3 ^ sin c 


D 

ot -a + ~~ [2a 
cr -2 \ 


2 o • 2 

, COS a + 2a„ sin a - a - a 

1 3 2 3 


a , = o 
3 cr 


a cr(°l COs2 “ + ° 3 sin2 “) “ i [(«*, - ° 3 ) 


\2 . 2 2 
1 sin a cos a 
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TABLE A. 11(a) . - FORMS OF P 2 v FOR VARIOUS FRACTURE CRITERIA, 
CRACK CONFIGURATIONS, AND STRESS STATES 
) 02 — 03, 03 2 0 and Oq — 


Q(E,o ) f 

= sin a da 


2V 4ir 


* 2,1 mcM S.Z 

where 


V Vi - 4a i a 3 

2 *. 


a, (a ) = COS 
1 cr 


-1 


5 fa ) = cos 
2 cr 


-1 




s 

or $ - - when a = 0 

d 2 ’ 


a(a ) = cos 

cr T 


D , = 


_L 


(1 - 0 . 5 v) 


.2 ’ 


D , vll - 0,25 »1 , „ , 

2 (1 - 0 . 5 v) 2 


Fracture 
cri teri on 


Maximum 
tensi 1 e 
stress 


Crack 

conf i gurati on 


Gri f f i th 

crack 

(GC) 


Penny-shaped 

crack 

(PSC) 


Quadratic equation coefficients 
for o n (E,a) > 0 


■, ■ (-, - » 2 ) 2 

*2 = “(°1 " ° 2) 2 “ 4( °1 " °2 


)o 


cr 


= 4o (o - a ) 
3 cr cr 2 


■, = d i(°i - “ 2 ) 

>2 = - D ,[(°. - ° 2 ) 


2 ±_ , , 

+ Dl ( °1 " °2 °cr 


a„ = 4o (a - o ) 
3 cr cr 2 
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TABLE A.II(a). - Concluded. 


Fracture Crack Quadratic equation coefficients 

criterion configuration f or 0 ^( Ea j > q 




a o = (°i 

2 cr 1 


"± ( V 
-2 (°1 V 


o ( a 
cr cr 


) 


Penny-shaped 

crack 

(PSC) 














TABLE A. 11(b) . - FORMS OF P 2 v FOR VARIOUS FRACTURE CRITERIA, CRACK 
CONFIGURATIONS, AND STRESS STATES 
[°3 > 0 .] 






















TABLE A. III. - FORMS OF_ P 2S FOR VARIOUS FRACTURE CRITERIA AND SELECTED 

CRACK CONFIGURATIONS 


[oj > o 2 , o 2 > 0 , 03 = 0 and o e = o cr .] 


“(£,0 ) „ r 

p = — = 2 U, 

2S 2 it 11 


-a 0 ± Wa A - 4a 


or = 


when a = 0 


where 0 < <& 

1 ~ 2 


a (a ) = cos 
1 cr 


a 0 (a ) = cos 
2 cr 


a(a ) = cos 1 \f® 
cr V 


After obtaining roots for a given stress state Z and o cr , care must 
be taken in evaluating da. The relation of o e to a Qr in the 
neighborhood of a must be known to obtain the proper limits of the 


i ntegral . 


3 " 1 - 


Fracture Crack 

criterion conf i gurati on 


Quadratic equation coefficients 
for o n (£,a) > 0 


Normal 
stress 
(shear- 
i nsensi ti ve 

cracks ) 3 


Independent a = 0 
of crack 

shape o 


a 2 = °1 - °2 


a = o - o 
3 2 cr 


Strain 

energy 

release 

rate, 



a i = 0 


2 2 

a 2 = °2 - °1 


2 2 
a, = a - o 
3 cr 2 


"l = ( °1 ‘ °2 ) 2 

(1.1215) 


a 2 = -2o z (o 1 - « 2 ) - 


(1.1215) 




2 2 

a = a - a 
3 cr 2 


a For cases where neither nor 4> 2 is between 0 and 1 , see ref. 2. 




TABLE A. III. - Concluded. 


Fracture 
cri ten on 

Crack 

configuration 

Quadratic equation coefficients 
for o n (E,a) 2 0 

Shetty 

Gri f f i th 

crack 

(GC) 

*1 °zf 

* 2 * "V (a i ' V ‘ ^ (”l ' V 



a 3 = °cr (o cr " °2 ) 


Gri f f i th 

notch 

(GN) 

^ (•, - ° 2 ) 2 
(1.1215C) v 7 

a 2 = - ° C r ( °l - °2 ) ’ - 2 (°1 " °z) 

2 cr (1.12150 v 7 

a„ = o (a - o ) 

3 cr cr 2 


Semi ci rcul ar 
crack 

■, ■ (“f 5 ) 2 (•, - 

H -- - V - (T) 2 (•' - ^ 



a_ = o (o - o ) 

3 cr cr 2 
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